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FOREWORD 


The  Earth's  magnetic  field  continues  to  play  an  essential  role  in  global  navigation.  All 
navigational  aids  or  Attitude/Heading  Reference  Systems  (AHRS),  regardless  of  their  operating 
principles,  must  speak  a  common  language.  That  common  language  is  in  terms  of  the  Earth's 
directional-field  components,  magnetic  declination  and  magnetic  inclination.  Magnetic-related 
navigational  aids  are  integrated,  in  the  forms  of  computer  hardware  and  software,  into  virtually 
every  major  weapons  system  of  the  Army,  Air  Force,  Navy  and  Marines,  as  well  as  many  North 
Atlantic  Treaty  Organization  (NATO)  systems.  Most  particularly,  all  Department  of  Defense 
sponsored  Global  Positioning  System  (GPS)  receivers  have  the  World  Magnetic  Model 
imbedded  into  them.  In  order  to  maintain  optimum  performance,  these  AHRS  must  be 
periodically  updated  with  respect  to  the  Earth's  magnetic  field,  which  is  a  dynamic  entity  that 
changes  slowly  but  erratically  with  time. 

For  well  over  a  century  it  has  been  the  responsibility  of  the  Naval  Oceanographic  Office  to 
monitor  the  Earth's  changing  magnetic  field  and  periodically  report  on  these  changes  in  the 
forms  of  magnetic  charts  and  mathematical  models.  During  the  past  forty-five  years,  this  task 
has  involved  an  intensive  data  collection  effort  through  the  Navy's  Project  MAGNET  program 
and  more  recently  through  the  Navy's  Polar  Orbiting  Geomagnetic  Survey  (POGS)  satellite 
program.  The  Navy's  participation  in  this  effort,  however,  will  be  discontinued  as  of  1  October 
1995,  at  which  time  the  World  Magnetic  Modeling  program  will  be  returned  to  its  sponsor,  the 
Defense  Mapping  Agency,  for  further  disposition,  while  the  rest  of  the  Navy's  magnetic  data 
collection  and  analyses  efforts  will  be  discontinued.  Extreme  gratitude  is  expressed  to  all  those 
who  have  participated  in  the  Navy's  magnetics  program  in  general  and  its  Project  MAGNET 
program,  which  began  in  1951,  in  particular. 

This  report  is  a  comprehensive  summary  of  the  cooperative  effort  between  the  Naval 
Oceanographic  Office  and  the  British  Geological  Survey  in  producing  the  1995  Epoch  World 
Magnetic  Model,  WMM-95. 
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ABSTRACT 


The  1995  Epoch  World  Magnetic  Model  (WMM-95)  is  the  product  of  a  cooperative  modeling 
effort  between  the  United  States  Naval  Oceanographic  Office  (NAVOCEANO)  and  the  United 
Kingdom's  British  Geological  Survey  (BGS).  The  model  is  based  on  data  from 
NAVOCEANO's  Project  MAGNET  high-level  (>  15,000  ft.)  vector-aeromagnetic  surveys  flown 
from  1988  to  1994,  scalar  total  intensity  data  collected  from  NAVOCEANO's  Polar  Orbiting 
Geomagnetic  Survey  (POGS)  satellite  during  the  period  1991  through  1993,  and  all  available 
magnetic  observatory  vector  annual-means  and  repeat  station  data  supplied  by  countries  around 
the  world.  The  WMM-95  model  describes  that  portion  of  the  geomagnetic  field  generated 
within  the  Earth's  liquid  core.  It  is  presented  in  terms  of  one  set  of  spherical  harmonic 
coefficients  which  characterize  the  Earth's  Main  magnetic  field  at  the  fixed  epoch  1995.0  and  a 
second  set  of  spherical  harmonic  coefficients  which  predict  the  Secular  Variation  (i.e.,  slow 
temporal  change)  of  the  Earth's  Main  magnetic  field  between  1995.0  and  2000.0.  These 
coefficients  (Main  and  Secular  Variation)  are  collectively  called  Gauss  coefficients  in  honor  of 
the  German  scientist  Carl  Frederick  Gauss,  who  in  1834  created  the  first  spherical-harmonic, 
geomagnetic  field  model  and  used  it  to  show  that  the  bulk  of  the  geomagnetic  field  originated 
below  the  Earth's  surface. 

This  report  contains  a  detailed  summary  of  the  data  used,  analyses  performed,  modeling 
techniques  employed,  and  results  obtained  during  the  course  of  the  1995  Epoch  World  Magnetic 
Modeling  effort.  This  report  also  contains  the  GEOMAG  algorithm  and  describes  its  uses  and 
limitations.  Charts  derived  from  the  WMM-95  model  and  the  GEOMAG  algorithm  for  both  the 
Main  geomagnetic  field  components  and  their  Secular  Variations  are  presented  on  Mercator  and 
polar  stereographic  projections.  Additionally,  the  numerical  values  of  the  Main  geomagnetic 
field  components  and  their  Secular  Variations  are  tabulated  on  a  5-degree  worldwide  grid. 

The  Defense  Mapping  Agency's  Hydrographic/Topographic  Center  (DMA/HTC)  publishes 
wall-sized  charts  of  the  Total  Intensity,  Declination,  Inclination,  and  the  Horizontal  and  Vertical 
geomagnetic  field  components  and  their  respective  secular  variations  on  Mercator  and  polar 
stereographic  projections. 
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THE  1995  EPOCH  GEOMAGNETIC  MODEL  AND  THE  GEOMAG  ALGORITHM 


1 


1.0  Introduction 


The  Earth's  magnetic  field,  as  measured  by  a  magnetic  sensor  on  or  above  the  Earth's  surface,  is 
actually  a  composite  of  several  magnetic  fields  generated  by  a  variety  of  sources.  These  fields 
are  superimposed  onto  each  other  and,  through  inductive  processes,  interact  with  each  other. 
The  most  important  of  these  geomagnetic  sources  are: 

a.  the  Earth's  conducting,  fluid  outer  core; 

b.  the  Earth's  crust/upper  mantle; 

c.  the  ionosphere;  and 

d.  the  magnetosphere. 

More  than  90  percent  of  the  geomagnetic  field  is  generated  by  the  Earth's  outer  core.  It  is  this 
portion  of  the  geomagnetic  field  that  is  represented  by  the  1995  Epoch  World  Magnetic  Model 
(WMM-95).  Those  portions  of  the  geomagnetic  field  not  represented  by  the  model  are 
collectively  referred  to  as  the  anomalous  geomagnetic  field,  which  varies  both  spatially  and 
temporally  with  respect  to  the  model. 

The  model  itself  consists  of  a  degree  and  order  12  spherical-harmonic  Main  (i.e., 
core-generated)  Field  (MF)  model  comprised  of  168  spherical-harmonic  Gauss  coefficients  and  a 
degree  and  order  12  spherical-harmonic  Secular-Variation  (SV)  (core-generated,  slow  temporal 
variation)  field  model  comprised  of  an  additional  168  spherical-harmonic  Gauss  coefficients. 
As  in  previous  World  Magnetic  Models,  the  1995.0  Predictive  SV  Gauss  coefficients 
corresponding  to  degree  and  order  9,  or  larger,  were  set  to  zero  due  to  a  lack  of  data. 

The  primary  geomagnetic  data  set  for  the  1995  Epoch  MF  model  is  the  Polar  Orbiting 
Geomagnetic  Survey  (POGS)  satellite  data  set,  which  provided  long-term  Total  Intensity  data 
from  1991  through  1993.  This  data  set  was  supplemented  by  Project  MAGNET 
vector-aeromagnetic  data,  which  spanned  the  period  from  1988  through  1993.  These  two  data 
sets  and  the  geomagnetic  observatory  annual  means  data  set  also  provided  sufficient  spatial  and 
temporal  coverage  to  permit,  for  the  first  time,  the  computation  of  the  Definitive  SV  Gauss 
coefficients  to  the  full  degree  and  order  12  for  the  definitive  1992.5  Epoch  SV  model,  which  in 
turn  played  a  key  role  in  determining  the  1995  Epoch  MF  model  coefficients. 

The  MF  Gauss  coefficients  characterize  the  geomagnetic  field  at  one  instant  of  time  called  the 
base  epoch,  which  for  WMM-95  is  1995.0.  The  predictive  SV  coefficients  characterize  the  slow 
rate  of  change  of  the  geomagnetic  field  for  the  5-year  period  from  the  base  epoch  to  the 
termination  epoch,  which  for  WMM-95  is  2000.0,  at  which  time,  WMM-95  will  be  replaced  by 
WMM-2000.  The  magnetic  field  components  are  computed  via  the  magnetic  variation 
algorithm  (GEOMAG),  which  is  a  FORTRAN  (also  translated  into  C,  ADA,  and  other 
high-level  languages)  subroutine  that  uses  the  WMM  Gauss  coefficients  in  conjunction  with  the 
spherical-harmonic  expansions  associated  with  each  field  component. 

The  WMM  coefficients  are  produced  jointly  by  the  British  Geological  Survey  (BGS)  in 
Edinburgh,  Scotland,  and  the  Naval  Oceanographic  Office  (NAVOCEANO)  at  Steimis  Space 
Center,  Mississippi,  on  behalf  of  the  British  Hydrographic  Office  in  Taunton,  England,  and  the 
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Defense  Mapping  Agency  (DMA),  Washington,  D.C.  The  model,  associated  software,  and 
documentation  are  distributed  by  NAVOCEANO  on  behalf  of  DMA  in  accordance  with  DMA 
Instructions  8000.1  and  8000.2.  These  models  are  produced  at  5-year  intervals,  as  are  DMA's 
Declination/Grid-Variation  charts,  while  charts  of  the  other  geomagnetic  components  are 
published  by  DMA  every  10  years.  The  1995.0  Epoch  corresponds  to  a  10-year  interval  when 
all  of  the  geomagnetic  components  will  be  published  in  chart  form  by  DMA.  The  military 
specifications  for  the  WMM  are  contained  in  MIL-W-89500  (DMA  [1993]).  Magnetic  model 
requirements  that  are  more  stringent  than  those  set  forth  in  this  military  specification  (e.g.,  those 
which  must  include  magnetic  effects  of  the  Earth's  crust,  ionosphere,  or  magnetosphere  and/or 
require  greater  spatial  or  temporal  resolution  on  a  global,  regional,  or  local  basis)  should  be 
addressed  to: 

Director,  Defense  Mapping  Agency 

8613  Lee  Highway 

Fairfax,  VA  22031-2137 

ATTN:  PR,  ST  A- 13 

It  is  extremely  important  to  recognize  that  the  WMM  series  of  geomagnetic  models  and  the 
charts  produced  from  these  models  characterize  only  that  portion  of  the  Earth's  magnetic  field 
which  is  generated  by  the  Earth's  fluid  outer  core.  The  portions  of  the  geomagnetic  field 
generated  by  the  Earth's  crust,  upper  mantle,  ionosphere,  and  magnetosphere  are  not  represented 
in  these  models.  Consequently,  a  magnetic  sensor  such  as  a  compass  or  magnetometer  may 
observe  spatial  and  temporal  magnetic  anomalies  when  referenced  to  the  appropriate  WMM.  In 
particular,  certain  local,  regional,  and  temporal  magnetic  declination  anomalies  can  exceed  10 
degrees.  Anomalies  of  this  magnitude  are  not  common  but  they  do  exist.  Declination 
anomalies  on  the  order  of  3  or  4  degrees  are  not  uncommon  but  are  of  small,  spatial  extent  and 
are  relatively  isolated.  On  land,  spatial  anomalies  are  produced  by  mountain  ranges;  ore 
deposits;  ground  struck  by  lightning;  geological  faults;  and  cultural  features  such  as  trains, 
planes,  tanks,  railroad  tracks,  power  lines,  etc.  In  ocean  areas  these  anomalies  occur  most 
frequently  along  continental  margins;  near  seamounts;  and  near  ocean  ridges,  trenches,  and 
fault  zones,  particularly  those  of  volcanic  origin.  Ships  and  submarines  are  also  sources  of 
magnetic  anomalies  in  the  ocean. 

Temporal  anomalies  over  either  ocean  or  land  areas  can  last  from  a  few  minutes  to  several  days 
and  are  produced  by  ionospheric  and  magnetospheric  processes  which  are  driven  by  the  solar 
wind.  In  particular,  magnetic  storms  generated  by  solar  flares  and  other  solar  activity  can, 
through  modulation  of  the  solar  wind,  cause  severe  and  persistent  magnetic  anomalies  in  the 
Earth's  environment.  Even  during  periods  of  quiet  solar  activity,  significant  spatial  and 
temporal  magnetic  anomalies  are  found  in  the  polar  and  equatorial  regions  of  the  Earth,  where 
magnetic  fields  produced  by  ionospheric  current  systems,  such  as  the  auroral  electrojets  and  the 
equatorial  electrojet,  are  always  present.  Most  sources  of  magnetic  anomalies  are 
comparatively  isolated  in  either  space  or  time.  Therefore,  from  a  global  perspective,  the 
WMM-95  root-mean-square  (RMS)  Declination  (D),  Inclination  (I),  and  Grid  Variation  (GV) 
errors  of  the  WMM  are  estimated  to  be  less  than  0.5  degrees  in  ocean  areas  and  less  than  1.0 
degree  over  land  areas  at  the  Earth's  surface  over  the  entire  5-year  life  of  a  particular  model. 
Also,  the  RMS  errors  at  sea  level  for  the  full  5-year  life  of  the  WMM-95  model  for  the 


4 


Horizontal  Intensity  (H)  and  the  Vertical  component  (Z),  over  the  oceans,  are  estimated  to  be 
less  than  200  nanoTeslas  (nT),  while  for  the  Total  Intensity  (F)  the  RMS  error  is  estimated  to  be 
less  than  280  nanoTeslas.  Over  land  areas  the  H,  Z,  and  F  components  may  be  somewhat  larger 
and  are  more  difficult  to  determine.  So,  estimates  for  these  are  not  given. 


1.1  The  Mathematical  Model 

The  Earth's  core-generated  magnetic  field  has  associated  with  it  a  geomagnetic  potential 
V(r,6,(p,t),  which  can  be  expressed  in  spherical  coordinates  in  terms  of  a  spherical-harmonic 
expansion  of  the  following  form: 

V(r,Q,(p,t)  =  Re  J  X  iSnmit)  cosim(p)  +  sin(w(p)}P“(0)  (1) 

where  the  spherical  coordinates  (r,6,(p)  correspond  to  the  radius  from  the  center  of  the  Earth,  the 
colatitude  (i.e.,  90®  -  latitude),  and  the  longitude.  is  the  mean  radius  of  the  Earth  (6371.2 
km);  g„„(t)  and  h„Jt)  are  referred  to  as  the  Gauss  coefficients  at  time  t,  where  t  is  the  time  in 
years  (e.g.,  1997.312).  P^(d)  represents  a  particular  Schmidt-normalized  associated  Legendre 
polynomial  of  spherical-harmonic  degree  n  and  order  m.  These  are  polynomials  in  terms  of  the 
cosine  of  the  colatitude  6.  The  Gauss  coefficients  are  slowly  varying  functions  of  time  and  are 
expressed  in  the  form  of  a  Taylor  series  expansion,  where  only  terms  up  to  first  order  in  time 
are  retained  so  that: 


gnmif)  gnmiP Epoch)  Snm  T Epoch) 


■  Epoch 


^  t  ^  T Epoch  5 


(2a) 


^nm(^^)  ^nm(^PEpoch)  PEpoch) 


T Epoch 


^  t  <  Tppoch  +  5 


(2b) 


where  is  the  base  epoch  of  the  model,  which  for  WMM-95  is  1995.0.  Thus, 

and  are  the  Schmidt-normalized  Gauss  coefficients  of  the  WMM  at  the  model's  base 

epoch,  while  the  Schmidt-normalized  SV  Gauss  coefficients,  and  (pronoimced  dot 
and  dot,  where  the  dot  represents  differentiation  with  respect  to  time:  jj),  are  the  annual 

rates  of  change  of  the  MF  Gauss  coefficients  g„„  and  h„„  and  are  evaluated  at  the  middle  of  the 
model's  lifespan  (i.e.,  at  +  2.5).  The  MF  Gauss  coefficients  and  SV  field  Gauss 

coefficients  are  collectively  referred  to  as  spherical-harmonic  coefficients. 

Taking  the  time  derivative  of  eq.  (1)  yields  the  spherical-harmonic  expression  for  the 
Secular-Variation  V(r,d,(p,t)  of  the  geomagnetic  potential: 
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which,  due  to  the  assumption  that  the  MF  Gauss  coefficients  vary  linearly  with  time  during  the 
course  of  a  5-year  interval,  as  expressed  by  eqs.  (2a)  and  (2b),  is  approximately  time 
independent  over  this  short  time  span.  The  maximum  degree  N  of  the  spherical-harmonic 
expansions  in  eqs.  (1)  and  (3)  is  equal  to  12.  This  value  is  determined  by  noting  that  when  the 
spectral  density  of  the  MF  Gauss  coefficients  is  plotted  as  a  function  of  harmonic  degree,  a 
distinct  break  in  this  density  function  occurs  between  degree  12  and  degree  15.  This  is 
interpreted  to  mean  that  the  low-degree  harmonics  corresponding  to  A/^<  12  are  dominated  by 
core-generated  magnetic  fields,  while  those  high-degree  harmonics  for  which  N  >  15  are 
dominated  by  fields  generated  within  the  crust  and  upper  mantle.  These  fields  are  primarily 
associated  with  permanent  and  induced  magnetization.  This  magnetization  is  limited  to  depths 
for  which  the  ambient  temperature  does  not  exceed  the  Curie  temperature.  Spherical-harmonic 
degrees  13  and  14  correspond  to  a  transition  region  where  neither  magnetic  fields  generated 
within  the  Earth's  fluid  core  nor  those  generated  within  the  Earth's  crust  dominate.  We  therefore 
use  iV  =  12  as  the  spherical-harmonic  cutoff  which  permits  the  best  description  of  the 
core-generated  magnetic  field  and  its  slow  temporal  change.  This  means  that  the  shortest 
wavelength  contained  in  the  model  is: 

Kin  =  ^  =  3336  (4) 

Thus,  the  WMM  is  a  low-resolution  model.  High-resolution  (short  wavelength)  descriptions  of 
that  part  of  the  magnetic  field  generated  by  the  Earth's  upper  crust  are  better  characterized  via 
rectangular  harmonic  modeling  of  small  local  areas  (Quinn  and  Shiel  [1993]),  while 
intermediate  wavelength  descriptions  of  the  Earth's  magnetic  field  generated  by  the  lower  crust 
and  upper  mantle  are  best  characterized  via  spherical-cap  harmonic  models  of  large  regional 
areas  (Haines  [1985a,  1985b,  1985c,  and  1990]).  Global  geomagnetic  data  sets  currently 
available  do  not  support  high-resolution  (i.e.,  X  <  500  km)  models  and  only  marginally  support 
intermediate-resolution  (i.e.,  500  km  <  X  <  3336  km)  magnetic  field  modeling.  Special  local 
and/or  regional  magnetic  surveys  are  required  to  generate  intermediate-resolution  and 
high-resolution  geomagnetic  models.  Consequently,  there  are  some  applications  for  which  the 
use  of  the  WMM  will  be  entirely  inadequate. 

The  Earth's  magnetic  field  'B(r,d,<p,t)  is  a  vector  quantity  having  three  components  which 
correspond  to  the  projection  of  the  magnetic  field  vector  onto  the  three  coordinate  axes.  Thus, 
B/r,d,(p,t)  is  that  portion  of  the  field  pointing  radially  outward  from  the  Earth's  center  (i.e., 
perpendicular  to  the  surface  of  the  Earth);  Bg(r,d,(p,t)  is  that  portion  of  the  field  pointing  locally 
due  south;  and  B/r,e,(p,t)  is  that  portion  of  the  field  pointing  locally  due  east.  The  magnetic 
field  vector  can  be  computed  from  the  geomagnetic  potential  by  taking  its  negative  gradient, 
thus: 


B(r,  0,  (p,  t)  =  -  V  V(r,  0,  (p,  /)  (5) 

Consequently,  the  magnetic  field  components  are  related  to  the  geomagnetic  potential  as 
follows: 
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which  yield  the  following  spherical-harmonic  expansions: 


Br{r,e,<p,t)  =  2  {gnm(0cos(»i(p) -I- A„m(Osin(m9)}P”(0)  (7a) 


5e(r,e,(p,0  =  -  J  X  {gnmit)cos(m(p)  +  h„m(t)sin(m(p)} 


dP^m 


5(p(r,e,(p,0  =  iX  f'^1  m{gnmit)sm(m(f)  -  h„ni(t)cos{m(p)}P';i(Q)  (7c) 


These  expressions  are  solutions  of  the  Laplace  equation,  which  in  turn  is  derived  from 
Maxwell's  famous  electromagnetic  field  equations  under  the  assumptions  that  the  magnetic  field 
exists  in  a  source-free  region  (i.e.,  no  charges  or  currents  are  present),  and  that  the  fields  are 
slowly  varying. 

Similarly,  for  the  geomagnetic  S  V  field  we  have: 


B(r,  0, 9, 0  =  -  V  F(r,  0, 9,  t) 


so  that: 


5^(r,0,9,t)  - Yr — 
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which  yield  the  following  spherical-harmonic  expressions: 


Br{r,%,%f)  =  1)(^]  S  |^„m(0cos(/n(p)  +  /i„^(/)sin(w9)jp;r(e)  (10a) 
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1.2  Spherical-Harmonic  Normalization 

• 

The  Gauss  coefficients  g^(t),  hji),  g^(t),  and  hj^x),  as  well  as  the  associated  Legendre 
polynomials  and  their  derivatives,  are  Schmidt  normalized  by  international  agreement  (circa 
1930)  of  the  International  Union  of  Geodesy  and  Geophysics  (lUGG).  This  particular 
normalization  allows  one  to  determine  which  spherical-harmonic  terms  of  a  particular  model  are 
the  most  significant  simply  by  a  cursory  inspection  of  the  model  coefficients'  relative 
magnitudes.  The  Schmidt-normalized  associated  Legendre  polynomials  P„”'(Q)  are  related  to  the 
unnormalized  associated  Legendre  polynomials  F''”(6)  (note  position  of  indices)  by  the 
following  relation: 

PJ'(0)  =  5""”P"'”(0)  (11) 

The  Schmidt  normalization  factors  5"'”  and  the  unnormalized  associated  Legendre  polynomials 
F"”(d)  are  computed  via  recurrence  relations  as  follows  (Cain  et  al,  1967): 


pOO(0)  ^  1 

(12a) 

P'""(Q)  =  sin0P”-i''”-H0) 

m  =  n  ^  0 

(12b) 

P'""(0)  =  cos0P"“^’'”(0)  -  K'”"P”'^’'”(0) 

m  ^  n,  n  >  \ 

(12c) 
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1.3  Coordinate  Transformations 

Although  the  magnetic  field  model  is  defined  in  terms  of  spherical  coordinates,  the  intended 
application  is  in  geodetic  coordinates.  So,  a  coordinate  transformation  is  necessary  (Cain  et  al., 
1967).  The  1984  World  Geodetic  System  (WGS-84)  (DMA  [1991])  and  its  corresponding 
ellipsoid  are  used  as  the  reference  datum  for  this  purpose.  Computing  the  magnetic  field 
components  at  a  given  location  expressed  in  geodetic  coordinates  using  the  WMM-95  model  is  a 
three-step  procedure: 
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a.  Convert  the  geodetic  latitude,  longitude,  and  altitude  (X,(p,h)  to  spherical  coordinates 

(r,0,(p). 

b.  Compute  the  magnetic  field  components  BXr,6,(p,t),  Bg(r,0,(p,t),  and  B^(r,0,(p,t). 

c.  Rotate  the  magnetic  field  components  from  spherical  coordinates  back  to  geodetic 
coordinates,  thus  yielding  the  magnetic  field  components  Bx(X,(p,h,t),  By(A,,(p,h,t),  and 
B2(?t,<p,h,t),  which  are  projections  of  the  magnetic  field  vector  B(X,(p,h,t)  onto  the  X-north, 
Y-east,  and  Z-vertically  down  coordinate  axes  of  the  local  rectangular  coordinate  system 
defined  by  the  tangent  plane  to  the  ellipsoid,  which  is  concentric  about  the  WGS-84  ellipsoid 
and  which  encompasses  the  point  (X.,(p,h). 

The  transformations  corresponding  to  step  a  are  as  follows: 


COS0  = 


cos2X,+  sin^X 


(16a) 


sin0  =  yi  -cos^0  (16b) 

where,  if  a  and  b  are  respectively  the  semi-major  and  semi-minor  axes  of  the  WGS-84 
ellipsoid,  then: 


Furthermore: 


-  (a^  -  6^)  sin^X 


a*  -{a* -b*)  sin^X 
g2  _  (^2-^2)  sin^X 


(18) 


The  transformations  corresponding  to  step  c  depend  on  the  angle  a  through  which  the  magnetic 
field  vector  must  be  rotated  while  transforming  from  spherical  to  geodetic  coordinates.  This 
rotation  angle  is  defined  by  the  following  transformation  equations. 


h  +  -Ja^  cos^X+b^  sin^^ 


(19a) 


sin  a 


(a^-b^)  cos  X  sin  X. 


(19b) 


a  =  +  0  -  f 


(19c) 


Consequently,  the  components  of  the  magnetic  field  vector  in  geodetic  coordinates  may  be 
computed  as  follows: 


10 


BxO^,<^,h,t)  =  -cosa^eCr,  0,(p,/)  -  sina5r(r,6,  (p,^) 


(20a) 


5y(X,,  9,^,0  =  5(p(r,0,(p,O 


(20b) 


5z(X.,  9, /z,  ?)  =  sina5e(^,0, 9,0  ~  cos  a5r(r,  0,9,0 


(20c) 


From  these  three  rectangular  geomagnetic  field  components,  it  is  possible  to  compute  all  others. 
In  particular,  the  following  magnetic  components  can  be  computed: 

^i/(^, 9,^,0  =  jB^X,(p,h,t)  +  (Horizontal  Intensity)  (21a) 


B fQ^,  9,  h,  0  =  9,  h,  t)  +  B^(K,  9,  h,  t)  (Total  Intensity) 


(21b) 


BD(^,(?,h,t)  =  tan' 


J  By(X,(?,h,t) 

1  Bx(X,(p,h,t) 


(Declination) 


Bi(k,(^,h,t)  =  tan 


-1  f 

1 


(Inclination) 


(21d) 


Bd  ~  9  X  >  0 
Bd  9  X  <  0 


(Grid  Variation) 


Inclination  is  often  referred  to  as  the  Dip  angle,  while  the  magnetic  declination  is  sometimes 
referred  to  as  the  magnetic  variation.  Frequently,  the  magnetic  field  components  in  eqs.  (21a) 
through  (21e)  are  simply  referred  to  in  terms  of  their  subscripts:  X,  Y,  Z,  H,  F,  D,  I,  and  G.  The 
Total  Magnetic  Intensity  is  sometimes  referred  to  as  77,  while  the  Grid  Variation  is  sometimes 
referred  to  as  GV.  Some  additional  and  quite  useful  relationships  among  these  magnetic  field 
components  are: 

H(k,  9,  h,  t)  =  F(X,  9,  h,  t)  cos[7(X,  9,  h,  t)]  (22a) 


X(K,  9,  h,  t)  =  H{%,  9,  h,  t)  cos[D(X,  9,  h,  t)] 


(22b) 


7(X,  9,  h,  t)  =  H(k,  9,  h,  t)  sin[D(X,  9,  h,  /)] 
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Z(X,  (p,  h,  t)  =  ¥(k,  cp,  h,  t)  sin  (p,  h,  t)]  (22d) 

1.4  The  GEOMAG  Algorithm 

The  Main  Field  Gauss  coefficients  at  the  base  epoch,  are  stored  in  array  C  of  the 

GEOMAG  algorithm  (sometimes  referred  to  as  the  MAGYAR  algorithm),  which  is  listed  in  the 
appendix,  such  that  the  lower  half  of  array  C  is  occupied  by  those  Gauss  coefficients  SnmC^Epoch) 
corresponding  to  the  cosine  terms  in  the  potential  function  of  eq.  (1),  while  the  upper  half  of 
array  C  is  occupied  by  those  Gauss  coefficients  h^(T£p^i,)  corresponding  to  the  sine  terms  in 
eq.  (1).  Table  1  illustrates  the  details  of  this  storage  scheme,  which  is  equivalent  to  the 
following  mathematical  assignments: 

^  _  J 

^nm  ""  I  1 

[  nm,n+l 

(23) 

which  implies  that: 

g„m  =  C„„.  m  <  n  (24a) 


hnm  =  Cm-i.n  m  <  n,  m  0  (24b) 

The  Secular- Variation  Gauss  coefficients  which  describe  the  Main  Field's  slow  annual  change 
are  stored  in  array  CD  (which  stands  for  C  [pronounced  C  dot])  such  that  the  lower  half  of  array 
CD  is  occupied  by  the  Gauss  coefficients  g.^ ,  which  correspond  to  the  cosine  terms  in  eq.  (3), 

while  the  upper  half  of  the  array  is  occupied  by  the  Gauss  coefficients  h,„„,  corresponding  to 
the  sine  terms  in  eq.  (3).  Table  2  illustrates  the  details  of  this  storage  scheme  for  array  CD.  It 
takes  essentially  the  same  form  as  table  1  for  array  C  and  corresponds  to  the  following 
mathematical  assignments: 

c.„  =  1  f "  “  -  "  I  (25) 

[  hm.n+1  m  >  n 
which  implies  that: 

=  C„m  m  <  n  (26a) 


m-l,n 


m  <  n,  m  9^  0 


(26b) 
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Table  1.  Arrangement  of  Main  Field  Coefficients  in  Array  Cnm 


10  11  12 


h\D,\  All, I  Ai2,i 

A  10,2  A  11,2  A  12,2 

Aio,3  All, 3  Ai2,3 

A  10,4  A  11,4  A  12,4 
Aio,5  An,5  Ai2,5 

Aio,6  All, 6  Ai2,6 

A  10,7  A  11,7  A  12,7 

Aio,8  An, 8  Ai2,8 

A  10,9  A  11,9  A  12,9 
Aio,io  All, 10  Ai2,io 
gl0,10  All, 11  Ai2,11 
g\\,\0  gu,\\  Ai2,12 


g'la.io  gn,n  gn,i2\ 


Table  2.  Arrangement  of  Secular  Variation  Coefficients  in  Array 


n\m 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

0 

^00 

• 

hu 

• 

A21 

• 

hsi 

• 

h4i 

• 

hsi 

• 

h6i 

• 

hi] 

• 

^81 

• 

A91 

• 

hio,i 

• 

hii,i 

• 

^12,1 

1 

S'lo 

Su 

• 

^22 

• 

hii 

• 

^42 

• 

^52 

• 

hoi 

m 

hi2 

• 

ho2 

• 

^92 

• 

hio,2 

• 

^11,2 

• 

h\2,2 

2 

S20 

k} 

^22 

• 

A33 

• 

h43 

• 

/?53 

m 

h63 

• 

hi3 

• 

ho3 

• 

A  93 

• 

hl0,3 

• 

^11,3 

• 

hl2,3 

3 

S30 

hi 

^32 

^33 

• 

/Z44 

• 

hs4 

• 

A  64 

• 

hi  4 

00 

4^ 

• 

^94 

• 

hl0,4 

• 

hii,4 

• 

^12,4 

4 

ho 

§4] 

S42 

S43 

5^44 

• 

hss 

• 

^65 

• 

hi  5 

• 

ho5 

• 

^95 

• 

^10,5 

• 

^11,5 

• 

h]2,5 

5 

S50 

4 

hi 

hi 

§53 

^54 

§55 

• 

^66 

• 

hi  6 

• 

ho6 

• 

^96 

• 

hl0,6 

• 

hii,6 

• 

hl2,6 

6 

Seo 

§61 

Sei 

§63 

§64 

§65 

S66 

• 

hii 

• 

hoi 

• 

A  97 

• 

hio,i 

• 

hii,i 

• 

h]2,i 

7 

Sio 

hi 

S12 

§13 

§14 

§15 

§16 

§11 

• 

^88 

• 

h90 

• 

^10,8 

• 

hii,o 

• 

hl2,0 

8 

^80 

Ssi 

s%2 

§%3 

§04 

§05 

§06 

§01 

Sgg 

• 

A  99 

• 

hl0,9 

• 

^11,9 

• 

hl2,9 

9 

^90 

S91 

S92 

§93 

§94 

§95 

§96 

§91 

S98 

S99 

• 

^10,10 

• 

^11,10 

• 

^12,10 

10 

• 

^10,0 

Sio,i 

Sl0,2 

§10,3 

• 

^10,4 

S’ 10,5 

• 

Si  0,6 

§10,1 

• 

Si  0,8 

• 

Si  0,9 

• 

Sio,io 

• 

hii,ii 

• 

^12,11 

11 

S’n,o 

Sn,i 

Sll,2 

§11,3 

§11,4 

Si  1,5 

Si  1,6 

§11,1 

Sll,8 

Sll,9 

• 

Sii,io 

Sii,ii 

• 

h\2,l2 

12 

^12,0 

Sl2,l 

Sl2,2 

§12,3 

§12,4 

S  12,5 

Sl2,6 

§12,1 

Sl2,8 

Si  2,9 

• 

Si2,10 

Sl2,ll 

• 

Sl2,12 
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The  numerical  values  of  the  Gauss  coefficients  at  the  base  epoch  and  their  corresponding 
predictive  annual  rates  of  change  for  the  WMM-95  geomagnetic  model  are  listed  in  table  3. 
These  numerical  values  are  inserted  into  arrays  C  and  CD  through  data  statements  in  the 
GEOMAG  algorithm.  Replacing  the  Gauss  coefficients  in  these  data  statements  and  the  date  of 
their  base  epoch  are  the  only  changes  that  need  to  be  made  to  update  the  algorithm  from  an  older 
model  to  the  new  model.  In  all  other  respects  the  GEOMAG  routine  remains  unaltered.  Other 
versions  of  the  GEOMAG  routine  exist  for  which  the  coefficients  can  be  read  in  from  an 
external  file.  Then,  only  the  external  coefficient  data  file  needs  to  be  updated,  while  the 
algorithm  remains  unchanged,  except  for  comment  statements  that  may  be  revised. 

Important  parameters  in  the  GEOMAG  routine  and  their  mathematical  correspondences  are. 


A 

a  =  6378.137  Am 

B 

- 

b  =  6356.7523142 

RE 

- 

7?£  =  6371.2  km 

TIME 

t 

EPOCH 

-- 

T 

^  Epoch 

DT 

Epoch 

ALT 

h 

SNORM(N,M) 

-- 

Snm 

K(N,M) 

K"” 

GLAT 

X 

GLON 

- 

<P 

SP(M) 

- 

sin(m(p) 

CP(M) 

- 

cos(m(p) 

ST 

-- 

sin  (9) 

CT 

- 

cos(9) 

CA 

cos(oc) 

SA 

sin  (a) 

BR 

- 

Br 

BT 

Be 

BP 

BX 

Bx 

BY 

By 

BZ 

Bz 

DEC 

-- 

Bo 

DIP 

B, 

TI 

- 

B, 

MAXDEG 

N 

MAXORD 

M  =  N 

P(N,M) 

DP(N,M) 

TC 

- 

C  +  (t  - 

CD 

- 

C 

Q2 

- 
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Table  3.  WMM-95  Model  Coefficients 


n 

m 

gn 

urn 

•  m 

gn 

•  m 

bn 

1 

0 

-29,682.1 

0.0 

17.6 

0.0 

1 

1 

-1,782.2 

5,315.6 

13.2 

-18.0 

2 

0 

-2,194.7 

0.0 

-13.7 

0.0 

2 

1 

3,078.6 

-2,359.1 

4.0 

-14.6 

2 

2 

1,685.7 

-418.6 

-0.3 

-7.2 

3 

0 

1,318.8 

0.0 

0.8 

0.0 

3 

1 

-2,273.6 

-261.1 

-6.6 

4.0 

3 

2 

1,246.9 

301.0 

-0.5 

2.2 

3 

3 

766.3 

-416.5 

-8.5 

-12.6 

4 

0 

940.0 

0.0 

1.2 

0.0 

4 

1 

782.9 

259.4 

1.1 

1.3 

4 

2 

290.9 

-230.9 

-6.8 

1.0 

4 

3 

-418.9 

99.8 

0.3 

2.5 

4 

4 

113.8 

-306.1 

-4.5 

-1.2 

5 

0 

-209.5 

0.0 

0.9 

0.0 

5 

1 

354.0 

43.7 

0.5 

0.5 

5 

2 

238.2 

157.6 

-1.4 

1.5 

5 

3 

-122.1 

-150.1 

-1.7 

0.6 

5 

4 

-162.8 

-59.2 

0.0 

1.7 

5 

5 

-23.3 

104.4 

2.1 

0.6 

6 

0 

68.5 

0.0 

0.4 

0.0 

6 

1 

65.6 

-15.2 

-0.3 

0.7 

6 

2 

64.1 

74.3 

0.3 

-1.5 

6 

3 

-169.1 

69.4 

2.1 

-0.5 

6 

4 

-0.5 

-55.3 

0.0 

-0.7 

6 

5 

16.5 

3.0 

-0.4 

1.1 

6 

6 

-91.0 

33.3 

-0.4 

2.6 

7 

0 

78.0 

0.0 

-0.3 

0.0 

7 

1 

-68.1 

-76.1 

-1.1 

0.3 

7 

2 

0.1 

-24.5 

-0.5 

0.0 
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Table  3.  WMM-95  Model  Coefficients  (Con.) 


n 

m 

pm 

fen 

K 

•  m 

gn 

•  m 

h„ 

7 

3 

29.6 

1.6 

0.5 

0.7 

7 

4 

6.0 

20.0 

1.3 

-0.6 

7 

5 

8.7 

16.5 

0.1 

0.1 

7 

6 

9.2 

-23.6 

0.0 

-0.6 

7 

7 

-2.4 

-6.8 

-0.9 

-0.4 

8 

0 

24.7 

0.0 

0.1 

0.0 

8 

1 

3.4 

14.9 

0.0 

0.4 

8 

2 

-1.5 

-19.5 

0.4 

-0.3 

8 

3 

-9.6 

6.3 

0.3 

0.1 

8 

4 

-16.5 

-20.4 

-1.3 

0.8 

8 

5 

2.6 

12.2 

0.5 

-0.1 

8 

6 

3.6 

7.0 

0.4 

-1.3 

8 

7 

-4.9 

-19.0 

-0.9 

-0.9 

8 

8 

-8.5 

-8.8 

0.1 

-1.1 

9 

0 

2.9 

0.0 

0.0 

0.0 

9 

1 

7.5 

-19.8 

0.0 

0.0 

9 

2 

0.4 

14.6 

0.0 

0.0 

9 

3 

-10.3 

10.9 

0.0 

0.0 

9 

4 

9.7 

-7.5 

0.0 

0.0 

9 

5 

-2.3 

-6.8 

0.0 

0.0 

9 

6 

-2.4 

9.3 

0.0 

0.0 

9 

7 

6.8 

7.7 

0.0 

0.0 

9 

8 

-0.5 

-8.1 

0.0 

0.0 

9 

9 

-6.5 

2.6 

0.0 

0.0 

10 

0 

-2.9 

0.0 

0.0 

0.0 

10 

1 

-3.3 

3.2 

0.0 

0.0 

10 

2 

2.8 

1.7 

0.0 

0.0 

10 

3 

-4.3 

2.9 

0.0 

0.0 

10 

4 

-3.1 

5.6 

0.0 

0.0 

10 

5 

2.4 

-3.4 

. 

0.0 

0.0 
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Table  3.  WMM-95  Model  Coefficients  (Con.) 


n 

m 

pm 

fen 

K 

•  m 

gn 

•  m 

h„ 

10 

6 

2.8 

-0.7 

0.0 

0.0 

10 

7 

0.7 

-2.9 

0.0 

0.0 

10 

8 

4.1 

2.3 

0.0 

0.0 

10 

9 

3.6 

-1.6 

0.0 

0.0 

10 

10 

0.6 

-6.6 

0.0 

0.0 

11 

0 

1.7 

0.0 

0.0 

0.0 

11 

1 

-1.6 

0.3 

0.0 

0.0 

11 

2 

-3.6 

1.0 

0.0 

0.0 

11 

3 

1.2 

-3.6 

0.0 

0.0 

11 

4 

-0.6 

-1.4 

0.0 

0.0 

11 

5 

0.1 

1.9 

0.0 

0.0 

11 

6 

-0.7 

0.2 

0.0 

0.0 

11 

7 

-0.8 

-1.3 

0.0 

0.0 

11 

8 

1.3 

-2.4 

0.0 

0.0 

11 

9 

-0.3 

-0.6 

0.0 

0.0 

11 

10 

2.2 

-2.2 

0.0 

0.0 

11 

11 

4.2 

1.3 

0.0 

0.0 

12 

0 

-1.8 

0.0 

0.0 

0.0 

12 

1 

0.9 

0.3 

0.0 

0.0 

12 

2 

-0.1 

1.4 

0.0 

0.0 

12 

3 

-0.5 

0.8 

0.0 

0.0 

12 

4 

0.8 

-3.0 

0.0 

0.0 

12 

5 

0.2 

0.7 

0.0 

0.0 

12 

6 

0.5 

0.5 

0.0 

0.0 

12 

7 

0.4 

-0.8 

0.0 

0.0 

12 

8 

-0.4 

0.6 

0.0 

0.0 

12 

9 

0.3 

0.1 

0.0 

0.0 

12 

10 

0.2 

-1.3 

0.0 

0.0 

12 

11 

0.4 

-0.4 

0.0 

0.0 

12 

12 

0.6 

0.9 

0.0 

0.0 
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Note  that  is  not  intended  to  be  the  mean  radius  of  the  WGS-84  ellipsoid.  By  international 
convention  established  by  the  International  Association  of  Geomagnetism  and  Aeronomy 
(I  AG  A)  circa  1968,  it  is  the  mean  radius  of  a  modified  ellipsoid  established  by  the  International 
Astronomical  Union  (lAU)  in  1966.  This  ellipsoid  is  referred  to  as  the  modified  IAU-66 
ellipsoid  (Zmuda  [1971]). 

The  GEOMAG  algorithm  is  organized  into  two  modules,  each  with  its  own  entry  point.  The 
first  is  an  Initialization  Module.  Its  purpose  is  to  compute  all  constants  such  as  the  recursion 
relation  factors  for  the  associated  Legendre  polynomials  k™” ,  the  Schmidt  normalization  factors 
6™*,  and  any  other  parameters  that  do  not  depend  on  position  or  time.  The  entry  point  for  this 
module  is: 

GEOMAG(MAXDEG) 

The  parameter  MAXDEG  determines  the  maximum  degree  and  order  of  the  magnetic  model  to 
be  used  in  the  computations.  Normally,  MAXDEG=12,  which  is  the  maximum  degree  and  order 
of  the  WMM  series  of  geomagnetic  models.  In  order  to  reduce  computation  time,  MAXDEG 
may  be  set  to  a  number  less  than  12  (e.g.,  8  or  10).  However,  the  accuracy  of  the  computed 
magnetic  parameters  is  correspondingly  reduced.  MAXDEG  must  be  set  in  the  calling  program. 

The  second  module  is  the  Processing  Module,  which  has  the  following  entry  point: 

GEOMGl(ALT,GLAT,GLON,TIME,DEC,DIP,TI,GV) 

The  purpose  of  this  module  is  to  compute  the  magnetic  Declination,  Inclination,  Total  Intensity, 
and  the  Grid  Variation  at  each  geodetic  position  and  time  supplied  to  it.  The  units  of  the 
parameters  in  the  argument  list  of  the  GEOMGl  entry  point  are  as  follows: 


ALT 

kilometers 

(e.g.,  5.314) 

(IN) 

GLAT 

degrees 

(e.g.,  33.716) 

(IN) 

GLON 

degrees 

(e.g., -163.315) 

(IN) 

TIME 

years 

(e.g.,  1997.427) 

(IN) 

DEC 

degrees 

(e.g., -121.734) 

(OUT) 

DIP 

degrees 

(e.g.,  48.387) 

(OUT) 

TI 

nanoTeslas 

(e.g.,  35781.7) 

(OUT) 

GV 

degrees 

(e.g.,  51.768) 

(OUT) 

computed  magnetic  field  parameters 

are  referenced  to  the 

WGS-84  ellipsoid. 

The  last 
•1 

parameter,  GV,  is  the  Grid  Variation  which  is  computed  only  for  the  polar  regions  (i.e.,  above 
+55®  latitude  or  below  -55®  latitude).  Outside  of  these  regions,  a  default  value  of  -999.0  is 
dummied  in.  The  Grid  Variation  is  referenced  to  Grid  North  of  a  polar  stereographic 
projection.  The  model  is  considered  to  be  a  valid  representation  of  the  Earth's  core  magnetic 
field  at  geodetic  altitudes  ranging  from  the  ocean  bottom  to  +1000  km  for  all  geodetic  latitudes 
and  longitudes. 
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The  SV  computation  of  a  geomagnetic  component  at  a  fixed  time  t  =  x  is  accomplished  by 
making  two  calls  to  the  entry  point  GEOMGl,  one  at  time  t,  =  x  -  0.5  and  one  at  time 
tj  =  X  +  0.5,  where  t  is  expressed  in  years.  This  yields  the  Declination,  Inclination,  and  the 
Total  Intensity  at  two  different  times  spaced  one  year  apart.  Using  these  three  magnetic 
components,  any  other  magnetic  component  can  be  calculated  at  these  same  two  times  via  eqs. 
(21a)  through  (21e)  and  eqs.  (22a)  through  (22d).  The  SV  is  then  determined  by  differencing 
the  two  MF  values  of  a  particular  component.  For  example,  the  Horizontal  component's  SV 
is  computed  by  inserting  eq.  (22a)  into  the  following: 

/f(X,(p,/i,x)  =  -  iy(X,(p,/j,ii)]/At  (27) 


where 


At  =  t2  -  ti  =  lyear  (28) 

The  FORTRAN  code  for  the  internal  version  (i.e.,  coefficients  embedded  in  the  algorithm  as 
data  statements)  of  the  GEOMAG  algorithm  is  listed  in  the  appendix. 
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SECTION  2 

THE  1995  EPOCH  GEOMAGNETIC  MODEL  DERIVATION 
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2.0  Overview 


Three  major  data  sets  were  available  for  the  1995  epoch  modeling  effort.  These  were  scalar 
magnetic  data  from  NAVOCEANO's  POGS  satellite  collected  from  January  1991  through  July 
1993;  NAVOCEANO's  Project  MAGNET  high-altitude  (>  15,000  ft.,  or  >  4.57  km) 
vector-aeromagnetic  survey  data  collected  from  October  1988  through  December  1993;  and 
Geomagnetic  Observatory  vector  annual-magnetic-means  data  collected  between  1985  and  1994 
and  provided  by  many  cooperating  nations  throughout  the  world.  Considerable  effort  was 
expended  toward  both  the  acquisition  and  the  reduction  of  these  data  sets  in  order  to  make  them 
suitable  for  modeling  purposes.  The  data  acquisition/reduction  procedures  for  POGS  data  are 
described  by  Quinn  et  al.  (1993b),  while  those  for  Project  MAGNET  data  are  described  by 
Coleman  (1992).  Those  for  observatory  vector-annual-magnetic-means  data  are  described  by 
Macmillan  (1994). 

In  addition  to  the  usual  concern  of  isolating  the  core  field  from  the  magnetic  field 
measurements,  which  also  contain  ionospheric,  magnetospheric,  and  crustal  magnetic  field 
contaminations,  the  quality  of  the  world  magnetic  model  generated  from  these  data  are 
influenced  by  four  major  factors: 

a.  Data  age  relative  to  the  model  epoch; 

b.  Data  temporal  coherence; 

c.  Data  spatial  uniformity;  and 

d.  Data  spatial  density. 

None  of  the  three  available  data  sets  are  ideal  in  all  of  these  respects.  However,  the  bulk  of 
these  data  sets  do  overlap  in  time  and  are  fairly  recent  data,  being  less  than  5  years  old. 

Small  POGS  data  subsets  (i.e.,  having  a  20-  to  60-day  time  span)  may  be  considered  temporally 
coherent,  uniformly  distributed  in  space,  and  of  sufficient  spatial  density  for  degree  12 
spherical-harmonic  modeling.  The  POGS  satellite  data  distribution  used  in  the  1995  epoch 
modeling  sequence  is  illustrated  in  the  Aitoff  equal-area  charts  1  through  42.  The  data  shown 
are  specially  selected  for  geomagnetically  quiet  times.  The  data  selection  criteria  consists  of  the 
following; 


a.  Kp  <  2" 

b.  lDst(0)l  <  50  nT 

c.  Local  Time:  night-side,  7  PM  to  5  AM 

The  planetary  K  (Kp)  index  is  a  measure  of  solar  activity  and  ranges  from  0  to  9  with  0  being 
the  quietest  time.  The  Disturbance  storm-time  (Dst)  index  is,  at  least  classically,  a  measure  of 
the  strength  of  the  magnetospheric  Ring  current.  It  is  based  on  geomagnetic  measurements 
made  by  a  select  group  of  geomagnetic  observatories  located  near  the  geomagnetic  equator.  The 
Dst  index  is  referenced  to  the  geomagnetic  equator  [i.e.,  Dst(O)].  Its  value  at  higher 
geomagnetic  latitudes,  0,  has  been  taken  to  be  the  equatorial  value  provided  by  M.  Sugiura 
through  the  National  Geophysical  Data  Center  (NGDC),  multiplied  by  the  cosine  of  the 
geomagnetic  latitude,  and  is  denoted  as  Dst(0).  Note  that  inadvertently,  the  POGS  data  fdes 
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Chart  1.  POGS  Distribution:  1991,  Days  Oil  -  020,  F  Component  Residuals 
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Chart  2.  POGS  Distribution:  1991,  Days  021  -  030,  F  Component  Residuals 
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Chart  4.  POGS  Distribution;  1991,  Days  041  -  050,  F  Component  Residuals 
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Chart  5.  POGS  Distribution:  1991,  Days  051  -  060,  F  Component  Residuals 
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Chart  6.  POGS  Distribution:  1991,  Days  121  - 130,  F  Component  Residuals 
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Chart  7.  POGS  Distribution:  1991,  Days  131  -  140,  F  Component  Residuals 


•  •  ♦  • 

mm  mm . 

CN)rs)  o  m  m  m  m  m 
—  —  on  00  fs  CD  m 

A 


m  m  m  m  •  •  m  m  m 

^cocM-^mm^ojtn 

I  1  I  I 


. m  m  mm 

mmmmmmo-^  rgrsj 
^mtDrvooon^^ 

I  I  I  I  I  I  I  I  I  I 

V 


I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 


o 

o 


n 

•O 

d  ^ 

00  o 


m 

u 

£ 

m 

2  = 

^  O 


11 


^  o  cvj 


£  H 

^  H  S  ^ 
«  ®  H  S 


O 

a 

II  % 

0)  *u 
o 

£  2 


II 


O  *3  f- 

§  §« 
X  .-S  § 

•  VI  g  j 

a  «  u  e 
^  o  03  u 


31 


Chart  8.  POGS  Distribution;  1991,  Days  171  -  180,  F  Component  Residuals 
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Chart  9.  POGS  Distribution;  1991,  Days  181  -  190,  F  Component  Residuals 
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Chart  10.  POGS  Distribution:  1991,  Days  201  -  210,  F  Component  Residuals 
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Chart  11.  POGS  Distribution:  1991,  Days  21 1  -  220,  F  Component  Residuals 
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Chart  12.  POGS  Distribution:  1991,  Days  231  -  240,  F  Component  Residuals 
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Chart  13.  POGS  Distribution:  1991,  Days  241  -  250,  F  Component  Residuals 
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Chart  14.  POGS  Distribution;  1991,  Days  251  -  260,  F  Component  Residuals 
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Chart  15.  POGS  Distribution:  1991,  Days  261  -  270,  F  Component  Residuals 


ldlo  m  ld  •*••**♦•  * 

rsjcM  —  OLOLninLOiniDLOinm 
^^CTiGor^cDLn'^torsj— i 


. in  in  inin 

•  •inininiOLninLOinLno-^  ojrvj 
mLD'-cMro^LniDt'^ooai^^ 


^  c/3 

2< 

Hg 

^  OQ 

pq  u 

(i:  >. 

cn  00 

2 

Q  u 

cn  z 


Jv%€ 


N  \ 


y  •rw' 


■■-<  )  I 


z' 

v'  • 

SC*^  “— " 

09 

V  ' 

•  ' 

:  /h 


W  q  -  H 

..  o  O  *3 
II  «  o  g 


«  VI  c  j 

a  «  V  e 
Q  m  u 


Chart  16.  POGS  Distribution:  1992,  Days  001  -  010,  F  Component  Residuals 
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Chart  17.  POGS  Distribution:  1992,  Days  01 1  -  020,  F  Component  Residuals 
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Chart  18.  POGS  Distribution:  1992,  Days  021  -  030,  F  Component  Residixals 
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Chart  19.  POGS  Distribution:  1992,  Days  041  -  050,  F  Component  Residuals 
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Chart  20.  POGS  Distribution:  1992,  Days  121  -  130,  F  Component  Residuals 
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Chart  22.  POGS  Distribution:  1992,  Days  141  - 150,  F  Component  Residuals 
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Chart  24.  POGS  Distribution:  1992,  Days  161  -  170,  F  Component  Residuals 
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Chart  25.  POGS  Distribution:  1992,  Days  171  -  180,  F  Component  Residuals 
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Chart  26.  POGS  Distribution:  1992,  Days  181  - 190,  F  Component  Residuals 
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Chart  27,  POGS  Distribution:  1992,  Days  191  -  200,  F  Component  Residuals 
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Chart  28.  POGS  Distribution:  1992,  Days  201  -  210,  F  Component  Residuals 
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Chart  29.  POGS  Distribution;  1992,  Days  211  -  220,  F  Component  Residuals 
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Chart  30.  POGS  Distribution:  1992,  Days  221  -  230,  F  Component  Residuals 
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Chart  31.  POGS  Distribution:  1992,  Days  231  -  240,  F  Component  Residuals 
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Chart  32.  POGS  Distribution:  1992,  Days  241  -  250,  F  Component  Residuals 
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Chart  33.  POGS  Distribution:  1993,  Days  021  -  030,  F  Component  Residuals 
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Chart  34.  POGS  Distribution:  1993,  Days  031  -  040,  F  Component  Residuals 
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Chart  35.  POGS  Distribution;  1993,  Days  041  -  050,  F  Component  Residuals 
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Chart  36.  POGS  Distribution:  1993,  Days  121  -  130,  F  Component  Residuals 
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Chart  37.  POGS  Distribution:  1993,  Days  131  - 140,  F  Component  Residuals 


lolo  m  LO 

fVJCM  —  O 


inminLDLnmLninLD 

ODoorvCDLnvforsj-^ 


in 


A 


•  9  •  « 

•  «••••«•*  lo  to 

•  lOLOioiOLOLniniOLOO^  ojoo 
LO— '(Nro^mcDc^oocn^*-^ 

I  I  I  I  I  I  I  I  1  I  I  I  I  I 

V 


I  i 


^  o 

^  O 

§2 

-  J 

^  Q) 

g’S 

a  » 
II  'S 
0) 

^  <£ 

b.  S 


g  ^ 

CVI  o 

H  S 

m  VI 


n 

h 

*  » 

»  ® 
o 

II  « 
«  II 
£  v 
H  £ 

•  ^ 

""  § 

e  u 

M-O 

a>  e 

Ca  bd 


61 


Chart  38.  POGS  Distribution:  1993,  Days  141  -  150,  F  Component  Residuals 
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Chart  39.  POGS  Distribution:  1993,  Days  171  - 180,  F  Component  Residuals 
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Chart  40.  POGS  Distribution:  1993,  Days  181  -  190,  F  Component  Residuals 
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Chart  41.  POGS  Distribution:  1993,  Days  191  -  200,  F  Component  Residuals 


POGS  DISTRIBUTION 

(MAGNETIC  TOTAL  INTENSITY  RESIDUALS,  UNITS:  nT) 


LOLO  IT)  LO  ••••••••  • 

rsJCM  OlDlDLOLniDinLOU^Ul  •  • 


.  .  •  •  *  •  •  *  *1/110  UllO 

•  inioioiOLOLOL/iLnioo-^  ojrs3 


rMto<rLntDN.cDO>^^  ^  — 

I  I  I  I  I  I  I . 

V 


•  iw^  ■= 


i 

V  *  O 

•:  M  * 

‘A-j 


>  / 


A  /  /  •© 


Chart  42.  POGS  Distribution:  1993,  Days  201  -  210,  F  Component  Residuals 


contain  the  value  of  M.  Sugium's  Dst(O)  index  divided,  rather  than  multiplied,  by  the  cosine  of 
the  geomagnetic  latitude.  This  is  not  a  serious  problem  since  Sugiura's  original  equatorial  Dst 
index  can  be  recovered  by  multiplying  the  POGS  Dst  values  by  the  cosine  of  the  geomagnetic 
latitude.  One  simply  needs  to  be  aware  of  how  the  data  has  been  treated.  The  geomagnetic 
coordinates  used  for  this  purpose  were  from  the  1990  Epoch.  The  Dst  index  has  a  range  of 
±  1000  nT  and  perhaps  higher  during  extremely  active  periods. 

The  WMM-95  modeling  effort  actually  used  the  POGS  Dst  values  for  data  selection. 
Consequently,  the  data  selection  criteria  became  more  stringent  approaching  the  geomagnetic 
poles  than  originally  intended.  This  was  actually  beneficial  since  less  geomagnetically  active 
Auroral  Zone  data  were  selected,  while,  as  charts  1  through  42  indicate,  no  noticeable 
deterioration  in  the  selected  POGS  data  distribution  was  observed  near  the  magnetic  poles.  The 
color  scheme  on  the  POGS  data  distribution  charts  corresponds  to  magnetic  Total-Intensity 
residuals,  computed  with  respect  to  initial  POGS-based  models  which,  as  input  data,  used  the 
POGS  file  indicated  in  the  chart  legend  and  other  files  nearby  in  time.  As  indicated  in  each 
chart  legend,  the  POGS  data  are  organized  into  10-day  files.  For  example,  the  file: 
pogs2121130.dst  corresponds  to  the  year  1992  and  covers  days  121  through  130.  The  extension 
"dst"  indicates  that  the  equatorial  value  of  the  Dst  index  {divided  by  the  cosine  of  the 
geomagnetic  latitude)  and  the  Kp  index  are  both  attached  to  each  data  record  within  the  10-day 
data  file. 

On  the  other  hand,  Project  MAGNET  data  are  not  considered  temporally  coherent,  and  by 
design,  these  data  are  primarily  confined  to  latitude  bands  that  straddle  the  geomagnetic  equator 
although  excursions  to  the  North  Magnetic  Pole  and  to  the  South  Magnetic  Pole  were  made,  as  is 
indicated  in  Aitoff  equal-area  charts  43a-d  and  44a-d.  In  the  equatorial  band  the  data 
distribution  is  considered  reasonably  uniform  and  of  sufficient  density  for  degree  12  spherical- 
harmonic  modeling  when  combined  with  POGS  data.  Charts  43a-d  display  Project  MAGNET 
data  coverage  for  the  X,  Y,  Z,  and  F  magnetic  components  of  the  C32  projects  indicated  in  the 
chart  legends.  These  surveys  span  the  time  interval  from  the  fourth  quarter  of  calendar  year 
1988  through  the  third  quarter  of  calendar  year  1990.  Charts  44a-d  display  Project  MAGNET 
data  coverage  for  the  X,  Y,  Z,  and  F  magnetic  field  components  of  the  D32  projects  indicated  in 
the  chart  legends.  These  surveys  span  the  time  interval  from  the  first  quarter  through  the  fourth 
quarter  of  calendar  year  1993.  An  additional  set  of  Project  MAGNET  data  was  collected  during 
the  first  calendar  quarter  of  1994  but  was  not  processed  in  time  to  use  in  the  1995  Epoch  model. 
The  project  numbers  on  these  charts  correspond  to  Fiscal  Years  (FY),  which  start  in  October. 
Thus,  Project  352  corresponds  to  the  second  quarter  of  FY93.  The  "5"  is  an  aircraft  designator. 
It  corresponds  to  the  fifth  aircraft  used  in  the  Project  MAGNET  program.  The  color  code  on 
these  charts  corresponds  to  the  magnetic  field  residuals  with  respect  to  a  modified  1990  Epoch 
model  generated  in  May  1994.  It  will  be  discussed  shortly.  Note  that  the  Y-component  appears 
to  be  biased  for  many  flights.  This  appears  to  have  been  the  result  of  a  bias  in  the  inertial 
attitude  device  on  the  aircraft.  Consequently,  this  component  received  less  weight  in  the 
modeling  process.  In  contrast  to  the  POGS  data,  the  aeromagnetic  data  were  not  selected  on  the 
basis  of  either  the  Kp  or  Dst  index.  Instead,  the  aircraft  is  flown  at  night,  and  as  far  as  it  was 
possible  to  do  so,  data  were  collected  during  periods  of  low  solar  activity. 
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Chart  43a.  C32  Project  MAGNET  Distribution:  X  -  Component  Residuals 
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Chart  43c.  C32  Project  MAGNET  Distribution:  Z  -  Component  Residuals 
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Chart  44a.  D32  Project  MAGNET  Distribution:  X  -  Component  Residuals 


PROJECT  MAGNET  DISTRIBUTION 
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Chart  44b.  D32  Project  MAGNET  Distribution:  Y  -  Component  Residuals 


Chart  44c.  D32  Project  MAGNET  Distribution:  Z  -  Component  Residuals 
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Chart  44d.  D32  Project  MAGNET  Distribution:  F  -  Component  Residuals 


The  POGS  and  Project  MAGNET  data  are  available  to  the  general  public  through  the  NGDC  in 
Boulder,  Colorado.  The  combined  data  sets  are  useful  for  both  MF  and  SV  modeling  and  are 
intended  to  compliment  each  other  in  the  sense  that  vector  data  are  required  as  input  to  the 
model  along  the  magnetic  equator  in  order  to  avoid  Backus  effect  modeling  errors  (Stem  and 
Bredekamp  [1975]  and  Stem  et  al.  [1980])  that  would  odierwise  arise  if  only  scalar  Total 
Intensity  data  were  available,  as  would  be  the  case  if  POGS  data  alone  were  used. 

Observatory  annual-means  data  have  been  traditionally  used  to  compute  the  SV  portion  of  the 
WMM.  The  primary  deficiency  of  this  data  set  is  its  poor  spatial  distribution.  Even  when 
supplemented  with  Project  MAGNET  intersect  data,  the  resulting  SV  models  are  limited  to 
spherical-harmonic  degree  8.  Observatory  data  are  not  directly  used  for  MF  modeling  due  to 
the  unknown  cmstal  biases  that  exist  under  each  observatory.  If  the  Earth  is  divided  into  equal 
area  cells  on  the  order  of  10°  X  10"  at  the  equator,  then  as  seen  in  chart  45  from  the  BGS, 
roughly  half  of  these  cells  contain  no  observatories.  In  this  chart,  cells  with  no  observatories 
contain  an  asterisk,  while  those  with  observatories  have  no  asterisk.  A  great  paucity  of 
observatory  data  can  be  seen  especially  in  the  southern  hemisphere  and  all  ocean  areas  in 
general.  The  status  of  currently  operating  geomagnetic  observatories,  of  which  there  are  only 
approximately  180,  is  discussed  by  McLean  et  al.  (1994).  As  indicated  in  chart  46  (from 
McLean  et  al.  [1994])  there  is  a  significant  time  lag  on  the  order  of  1  to  4  years  in  observatory 
data  reporting.  This  chart  indicates  the  actual  distribution  of  operating  observatories  as  of 
March  1994. 

The  modeling  objective  is  to  create  two  spherical-harmonic  models.  One,  the  MF  model, 
characterizes  the  Earth's  core-generated  magnetic  field  at  the  1995.0  epoch.  The  other,  the  SV 
model,  characterizes  the  Earth's  core-generated  slow-temporal-change  field  for  the  time  span 
from  1995.0  to  2000.0.  The  SV  model  is  thus  centered  at  1997.5  and  is  referred  to  as  the  1997.5 
Epoch  SV  model.  It  is  a  predictive  model  since  it  is  computed  prior  to  1995.0.  The  combined 
1995.0  MF  model  and  the  1997.5  SV  model  are  referred  to  collectively  as  the  1995.0  Epoch 
World  Magnetic  Model.  This  model  will  expire  at  the  end  of  December  1999,  when  another 
model  based  on  new  data  will  replace  it. 

Given  the  objective,  and  the  available  data,  the  following  procedure  was  adopted: 

a.  Use  the  Observatory  annual-magnetic-means  data  and  other  data  as  available  (e.g.. 
Project  MAGNET  intersect  data  and  repeat  station  data)  to  create  two  Definitive  SV  models,  the 
first  covering  the  5-year  interval  1985.0  to  1990.0  and  the  second  covering  the  5-year  interval 
1990.0  to  1995.0.  These  are  then  referred  to  as  the  BGS  Definitive  1987.5  and  1992.5  Epoch  SV 
models,  respectively.  These  were  provided  by  the  BGS  using  annual  means  data  available 
through  May  1994. 

b.  Use  the  BGS  Definitive  1987.5  SV  model  to  adjust  the  WC-85  (revised)  model  MF 
coefficients  (Quinn  et  al.  [1991])  to  the  1990.0  Epoch.  The  resulting  MF  model,  when  merged 
with  the  BGS  Definitive  1992.5  Epoch  SV  model  is  designated  as  the  WMM-90  (modified) 
model.  This  model  is  considered  to  be  the  best  a  priori  model  estimate  for  the  1990  to  1995 
time  span. 
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Chart  46.  NGDC  Geomagnetic  Observatory  Distribution  &  Reporting  History 


c.  Select  POGS  data  on  the  basis  of  the  Kp  and  Dst  indices  and  according  to  local  time 
in  accordance  with  the  criteria  stated  above  for  each  1 0-day  file. 

d.  Locate  consecutive  groups  of  10-day  POGS  select  files  from  which  the  total  data 
selected  consisted  of  approximately  18,000  uniformly  distributed  records  and  for  which  a 
maximum  of  six  10-day  files  was  used.  A  total  of  14  time  groups  satisfying  these  criteria  were 
identified.  Some,  corresponding  to  magnetically  very  quiet  times,  contained  just  two  10-day 
files,  while  others,  corresponding  to  less  quiet  times,  contained  the  maximum  of  six  10-day 
files.  Some  groups  overlapped,  white  other  groups  were  separated  significantly  in  time. 

e.  Make  preliminary  MF  models  for  each  of  the  14  POGS  defined  time  groups  at  their 
respective  mean  epochs  and  use  these  models  to  isolate  the  majority  of  the  ionospheric  magnetic 
field  and  to  devise  a  Magnetic  Latitude  dependent  Ionospheric  Correction  for  each  10-day 
POGS  select  file. 

f.  Select  from  each  ionospherically  corrected  10-day  file  for  all  14  time  groups,  that 
portion  of  the  data  which  passes  over  Europe,  where  the  WMM-90  (modified)  model  is 
expected  to  be  the  most  accurate  due  to  the  large  concentration  of  geomagnetic  observatories  in 
that  region.  This  select  European  data  subset  covers  the  approximate  time  from  1991  through 
1993.7.  Then  use  the  POGS  scalar  magnetic  residuals,  computed  with  respect  to  the  WMM-90 
(modified)  model,  to  compute  the  POGS  magnetometer  Total  Intensity  Bias  and  Drift-Rate  via 
least  squares,  assuming  a  linear  drift. 

g.  Apply  the  Bias  and  Drift-Rate  corrections  and  the  ionospheric  corrections  to  the 
10-day  POGS  select  files  in  each  of  the  14  time  groups. 

h.  Adjust  the  entire  Project  MAGNET  data  set  to  the  mean  epoch  of  each  of  the  14 
POGS-defined  time  groups,  using  the  WMM-90  (modified)  model. 

i.  Make  MF  models  at  each  of  the  14  group  epochs. 

j.  Perform  a  least-squares  fit  over  the  14  epochs  for  each  MF  Gauss  coefficient  using 
1992.5  as  a  reference  time.  A  1992.5  Epoch  MF  model  results,  along  with  a  1992.5  SV  model 
designated  as  SV  1992.5  (A).  Combined,  two  coefficient  sets  are  referred  to  as  the  WMM-92.5 
(optimum)  model.  Both  the  MF  and  SV  portions  of  this  model  are  completely  defined 
through  degree  and  order  12.  An  alternate  procedure,  independent  of  the  POGS  data  and 
involving  observatory  data  and  Project  MAGNET  intersect  data,  generated  another  1992.5  SV 
model  designated  as  SV  1992.5  (B).  It  is  defined  only  through  degree  and  order  8  due  to  the 
poor  distribution  of  these  two  data  sets. 

k.  Use  the  SV  portion  of  the  WMM-92.5  {optimum)  model  to  adjust  the  MF  d  portion  of 
the  WMM-92.5  {optimum)  model  backward  2.5  years  to  1990.0.  Then  merge  the  SV  portion  of 
the  WMM-92.5  {optimum)  model  with  this  new  adjusted  1990.0  Epoch  MF  model  to  form  the 
WMM-90  {revise^  model.  Subsequently  adjust  the  WMM-90  {revised)  MF  model  forward  to 
the  1995.0  Epoch  using  the  average  of  the  SV  1992.5  (A)  and  SV  1992.5  (B)  models  truncated 
to  degree  8.  The  result  is  the  WMM-95  MF  model. 
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l.  Use  Magnetic  Observatory  and  Project  MAGNET  survey-track  intersect  data,  filling 
in  remaining  spatial  gaps  with  SV  values  generated  from  the  WMM-90  (revised)  model  of  step 
"j"  where  necessary,  to  generate  a  1997.5  Predictive  SV  model. 

m.  Combine  the  1995.0  Epoch  MF  model  generated  in  step  "j"  with  the  1997.5  Epoch 
Predictive  S  V  model  of  step  "1"  to  obtain  the  WMM-95  model. 

The  WMM-90  (revised)  model  is  a  by-product  of  the  WMM-95  modeling  process.  When 
truncated  to  degree  10,  it  becomes  a  candidate  model  for  the  1990.0  Epoch  Definitive 
Geomagnetic  Reference  Field  (DGRF-90)  model  of  the  lAGA,  while  the  WMM-95  model,  when 
truncated  to  degree  10,  becomes  a  candidate  model  for  the  International  Geomagnetic  Reference 
Field  (IGRF-95)  model  of  lAGA. 

The  four  SV  models  discussed  above  are  listed  in  table  4.  Three  of  these  SV  models  (i.e., 
1987.5,  1992.5(B),  and  1997.5)  were  generated  by  BGS,  using  the  techniques  described  by 
Macmillan  (1994).  The  Gauss  coefficients  for  all  of  these  models  are  based  on  Observatory 
annual  means  data.  Repeat  Station  data,  and  Project  MAGNET  intersect  data  and  so  are 
necessarily  truncated  to  degree  and  order  8. 


2.1  Secular- Variation  Data  Analysis 

Traditionally,  the  most  accurate  data  available  for  SV  modeling  have  been  Quiet-Day 
observatory  magnetic-annual-means,  the  first  time-derivative  of  which  provides  information 
concerning  the  slow  rates  of  change  (greater  than  one  year)  of  various  components  of  the 
Earth's  Main  magnetic  field  at  various  locations  (roughly  180)  scattered  around  the  world.  Due 
to  the  sparsity  and  spatial  nonuniformity  of  this  data  set,  it  is  possible  to  generate  only  a  degree 
and  order  8  spherical-harmonic  SV  model,  even  when  these  data  are  supplemented  with  other 
data  such  as  Project  MAGNET  intersect  difference  data  and  Repeat  Station  data,  since  these 
supplementary  data  are  also  sparse,  nonuniformly  distributed  and  generally  less  accurate  than 
the  annual  means  data.  Furthermore,  the  Predictive  SV  model  is  necessarily  based  on 
extrapolations  of  observatory  annual  means  at  each  observatory  site.  The  BGS  has  been 
primarily  responsible  for  generating  the  Predictive  SV  portion  of  the  WMM  series  of  magnetic 
models  as  was  the  case  for  this  1995.0  Epoch. 

Examples  of  observatory  annual  means  data  from  a  few  selected  sites  (Honolulu,  Huancayo, 
Pilar,  and  Rude  Skov)  for  the  X-north,  Y-east,  and  Z-vertically  down  components  of  the  Earth's 
magnetic  field,  as  well  as  for  the  D-Declination,  I-Inclination,  and  F-Total  Intensity  components 
of  the  Earth's  magnetic  field,  are  given  in  figures  la,  lb,  Ic,  Id,  le,  and  If  through  4a,  4b,  4c, 
4d,  4e,  and  4f  The  discontinuities  in  the  field  components  at  Honolulu  are  due  to  the 

repositioning  of  the  observatory  at  two  different  instances.  The  geographic  location  and 

elevation  of  these  observatories  are  given  in  the  legend  in  the  upper  left-hand  comer  of  these 
figures.  Notice  that  the  rate  of  change  of  one  or  more  of  the  field  components  at  several 

observatory  sites  has  reversed  direction  over  time  intervals  as  short  as  2  or  3  years.  This 

behavior  is  an  illustration  of  the  sudden,  unpredictable  nature  of  the  Earth's  core  magnetic  field 
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Figure  lb.  East  Y  Component  at  Honolulu  (HON) 
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Figure  Ic.  Vertical  Z  Component  at  Honolulu  (HON) 
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Figure  Id.  Declination  D  Component  at  Honolulu  (HON) 
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Figure  le.  Inclination  I  Component  at  Honolulu  (HON) 
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Figure  If.  Total  Intensity  F  Component  at  Honolulu  (HON) 
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Figure  2a.  North  X  Component  at  Huancayo  (HU A) 
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Figure  2b.  East  Y  Component  at  Huancayo  (HU A) 
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Figure  2d.  Declination  D  Component  at  Huancayo  (HU A) 


OBSERVATORY  ANNUAL  MEANS 


Figure  2e.  Inclination  I  Component  at  Huancayo  (HU A) 
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Figure  2f.  Total  Intensity  F  Component  at  Huancayo  (HU A) 
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Component  at  Pilar  (PIL) 
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Figure  3b.  East  Y  Component  at  Pilar  (PIL) 
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Vertical  Z  Component  at  Pilar  (PIL) 
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Figure  Sd.  Declination  D  Component  at  Pilar  (PIL) 
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Inclination  I  Component 
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Figure  3f.  Total  Intensity  F  Component  at  Pilar  (PIL) 
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Component  at  Rude  Skov  (RSV) 
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Figure  4b.  East  Y  Component  at  Rude  Skov  (RSV) 
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Figure  4c.  Vertical  Z  Component  at  Rude  Skov  (RSV) 
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Figure  4d.  Declination  D  Component  at  Rude  Skov  (RSV) 
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Figure  4e.  Inclination  I  Component  at  Rude  Skov  (RSV) 
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Figure  4f.  Total  Intensity  F  Component  at  Rude  Skov  (RSV) 


changes  and  explains  why  it  is  necessary  to  make  new  spherical-harmonic  geomagnetic-field 
models  at  5-year  intervals. 

Quiet-Day  annual-means  data  are  used  for  SV  modeling  because  magnetic  fluctuations  generated 
by  external  current  sources,  having  periods  of  1  year  or  less,  tend  to  average  out  to  zero  over 
the  course  of  1  year.  These  fluctuating  fields  originate  from  a  variety  of  current  systems  located 
external  to  the  Earth's  surface  in  the  ionosphere  (e.g.,  the  Solar  quiet  (Sq)  current  system, 
Sabaka  and  Baldwin  [1993])  and  in  the  magnetosphere  (e.g.,  the  Ring  current  system).  They 
also  originate  in  current  systems  which  couple  the  magnetosphere  to  the  ionosphere  (e.g.,  the 
Field- Aligned  current  systems).  These  external  current  systems  are  driven  by  the  solar  wind. 
Their  associated  fluctuating  magnetic  fields  induce  currents  internal  to  the  Earth's  surface  iu  the 
Earth's  crust  and  mantle.  These  induced  currents  in  turn  generate  corresponding  induction  fields 
which  emanate  from  Earth's  crust  and  mantle.  This  inductive  process  takes  place  as  a 
consequence  of  the  finite  conductivities  associated  with  both  the  crust  and  mantle.  Magnetic 
fluctuations,  from  both  external  and  induced  internal  origins  having  periods  greater  than  1  year, 
still  contaminate  the  annual-means  data.  The  most  significant  of  these  magnetic  fluctuations  are 
related  to  the  11 -year  and  22-year  solar  cycles.  The  long-term  magnetic  fluctuations  in  the 
annual-means  data  set  are  considered  to  be  small  and  are  generally  ignored  because  many 
observatories  do  not  have  a  sufficiently  long  history  for  a  detailed  analysis  of  their  annual 
means  data  and  because  the  external  fields  are  not  sufficiently  well  understood  to  generate 
accurate  external-field  correction  models.  The  penalty  for  not  removing  the  long-period  external 
and  induced  fields  is  an  increased  uncertainty  in  the  SV  Gauss  coefficients. 


2.2  Main  Field  Data  Analysis 
2.2.1  Observatory  Data 

The  Observatory  magnetic  annual-means  data  were  not  directly  used  in  the  MF  modeling 
because  those  data  contain,  in  addition  to  small  external  field  contributions,  some  rather  large 
local  and  regional  magnetic  biases  of  crustal  origin.  A  detailed  regional  survey  around  each 
observatory  site  would  be  necessary  to  isolate  and  remove  these  biases.  Alternatively,  these 
biases  could  be  computed  as  part  of  the  MF  modeling  process.  However,  this  approach  is 
somewhat  controversial  since  modeling  at  different  epochs,  not  far  apart  in  time,  will  yield 
substantially  different  observatory  magnetic  biases  at  one  site.  This  bias  variability  is  partly  due 
to  the  local  induction  field  that  exists  at  each  site.  This  field  depends  on  the  local  conductivity 
distribution,  which  is  generally  unknown,  and  on  the  time-dependent  external  driving  fields  from 
the  ionosphere  and  magnetosphere.  The  annual-means  data  were  used  indirectly  in  the  MF 
modeling  through  the  two  definitive  SV  models  referred  to  in  the  previous  section.  These  SV 
models  are  derived  from  the  first  differences  of  the  annual  means  data,  which  are  independent 
of  the  observatory  biases.  They  were  used  to  adjust  the  Project  MAGNET  vector-aeromagnetic 
data  forward  and  backward  in  time  to  match  the  mean  epochs  of  14  POGS  data  subsets  which 
covered  the  time  span  from  1991 .0  to  1993.7.  Or,  put  another  way,  the  entire  Project  MAGNET 
data  set  was  copied  14  times  with  each  copy  adjusted  to  one  of  the  14  epochs  of  a  particular 
POGS  data  subset.  These  adjustments  were  accomplished  using  the  WMM-90  (modified)  SV 
model. 
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2.2.2  Project  MAGNET  Data 


The  Project  MAGNET  vector-aeromagnetic  data  consists  of  two  data  subsets.  The  first  data  set 
consists  of  78  high-altitude  (>  15,000  ft.)  survey  flights  which  were  flown  from  October  1988 
through  September  1990.  This  data  set  corresponds  to  projects  with  the  C32  designation  and 
was  collected  prior  to  the  structural  remodeling  of  the  Project  MAGNET  aircraft  which  took 
place  during  1991  and  1992.  The  second  data  set  consists  of  56  high-altitude  survey  flights 
which  were  flown  from  January  through  December  1993.  This  data  set  corresponds  to  projects 
with  the  D32  designation.  The  processed  data  consists  of  one  vector  sample  every  2  seconds 
(i.e.,  0.5  Hz  sample  rate)  yielding  1,216,686  vector  records  for  the  C32  project  flights  and 
1,050,462  vector  records  for  the  D32  project  flights.  Approximately  every  one-hundredth  vector 
record  was  selected  for  modeling  purposes,  yielding  10,658  vector  records  for  the  C32  flights 
and  9,360  vector  records  of  the  D32  flights  for  a  total  of  20,018  vector  records.  All  2,267,148 
vector  records  were  employed  to  determine  the  statistics  for  weighting  the  data  selected  from 
each  flight. 

The  Project  MAGNET  data  calibration  and  data  reduction  procedures  are  described  by  Coleman 
(1992).  The  statistics  for  each  magnetic  component  of  each  flight  for  the  C32  projects  are 
given  in  table  5,  while  statistics  for  the  D32  projects  are  given  in  table  6.  The  statistics  for  those 
data  collected  prior  to  1990.0  were  computed  with  respect  to  the  WC-85  {updated)  model  listed 
in  table  7.  It  is  composed  of  the  WC-85  {revised)  model  listed  in  NAVOCEANO  Technical 
Report  No.  304  (Quinn  et  al.  [1991])  and  the  revised  BGS  definitive  1987.5  SV  model  listed  in 
table  4  of  this  report.  The  statistics  for  those  data  collected  post  1990.0  were  computed  using 
the  WMM-90  (modified)  model,  the  coefficients  for  which  are  listed  in  table  8.  This  model  is 
composed  of  the  WC-85  {revised)  MF  Gauss  coefficients  adjusted  to  the  1990.0  Epoch  using 
the  revised  BGS  definitive  1987.5  SV  model  and  combining  the  resulting  MF  coefficients  with 
the  BGS  definitive  1992.5  SV  model  listed  in  table  4.  The  RMS  values  thus  computed  for  each 
flight  were  used  as  weight  factors  in  the  modeling  process.  These  statistics  indicate  that  the 
Y-component  magnetic  data  has  a  strong  bias  probably  associated  with  a  Y-component  attitude 
bias  error  in  the  ring-laser  gyro.  Nevertheless,  the  Y-component  magnetic  field  data  was  used 
in  the  modeling  but  was  given  less  weight  in  accordance  with  the  statistics  displayed  in  table  5 
and  table  6. 

The  Project  MAGNET  aircraft  is  a  Lockheed  RP3D  Orion.  Prior  to  1990  its  instrument  suite 
consisted  of  a  Honeywell  fluxgate  vector  magnetometer;  a  Texas  Instruments  ASQ-81 
absolute-scalar,  metastable-helium  magnetometer;  a  GPS  receiver;  a  precision  barometric 
altimeter;  and  an  electrostatic  gyro  for  attitude  determination.  Post  1990  the  electrostatic  gyro 
was  replaced  with  a  ring-laser  gyro,  and  a  radar  altimeter  was  added  to  the  instrument  suite.  For 
post- 1990  modeling  only  the  radar  altimeter  data  were  used  for  height  determinations  since  the 
barometric  data  yielded  slightly  larger  RMS  magnetic  field  statistics  with  respect  to  a  given  field 
model  than  did  the  edited  radar  altimeter  data.  The  accuracies  of  the  instruments  composing 
this  Project  MAGNET  instrument  suite  are  given  by  Coleman  (1992).  Procedures  used  for  both 
calibrating  the  vector  magnetometer  and  compensating  the  vector  magnetometer  data  for  the 
magnetic  fields  generated  by  the  aircraft  itself  are  also  given  by  Coleman  (1992).  Data  from 
these  instruments  was  time  synced  to  within  10  milliseconds.  Surveys  were  flown  at  night  (local 
time)  to  minimize  solar-driven  external  field  effects  that  contribute  to  the  magnetic  Daily 


106 


Table  5.  Project  MAGNET  Magnetic  Field  Flight  Statistics 
wrt.  WMM-90  (modified),  Year  >  1990 
wrt.  WC-85  (updated).  Year  <  1990 
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0.079 

0.181 

0.197 

0.019 

0.117 

0.119 

-0.019 

0.134 

0.135 

0.144 

0.125 

0.190 

0.067 

0.132 

0.148 

0.025 

0.136 

0.139 

-0.022 

0.155 

0.157 

-0.003 

0.133 

0.133 

0.165 

0.234 

0.286 

-0.055 

0.134 

0.145 

0.142 

0.049 

0.150 

0.215 

0.081 

0.230 

0.109 

0.054 

0.122 

0.012 

0.090 

0.090 

-0.044 

0.082 

0.093 

-0.062 

0.082 

0.103 

0.017 

0.089 

0.090 

0.362 

0.539 

0.649 

-0.211 

0.753 

0.782 

-0.035 

0.459 

0.460 

-0.224 

0.593 

0.633 

-0.116 

0.632 

0.643 

0.116 

0.054 

0.128 

-0.207 

0.170 

0.268 

-0.404 

0.239 

0.469 

-0.229 

0.162 

0.281 

-26.5 

61.8 

67.3 

-14.1 

53.2 

55.1 

-28.9 

55.0 

62.2 

-0.9 

42.8 

42.8 

-64.5 

85.6 

107.2 

-20.6 

73.3 

76.1 

-22.1 

78.1 

81.2 

-18.2 

97.5 

99.2 

3.6 

106.4 

106.5 

-43.1 

47.7 

64.3 

16.0 

61.2 

63.2 

-33.6 

78.0 

84.9 

19.8 
59.6 

62.8 

-26.2 

60.8 

66.2 

-23.6 

52.3 

57.4 

-18.3 

56.2 

59.1 

17.8 
48.6 

51.8 

-68.9 

84.7 

109.2 

-18.2 

74.3 

76.5 

0*08 

Z'LL 

L’OZ- 

-13.3 

98.4 

99.3 

13,1 

106.9 

107.6 

-48.7 

42.8 

64.8 

28.6 

60.7 

67.1 

-30.1 

78.9 

84.4 

18.9 

60.6 

63.4 

35.2 

106.9 

112.5 

19.9 

74.1 

76.7 

-9.3 

80.7 

81.3 

88.0 

68.3 

111.4 

-51.4 

95.2 

108.2 

-12.7 

81.4 

82.3 

-24.1 

83.2 

86.7 

-15.5 

107.6 

108.7 

63.1 

125.3 

140.3 

-11.4 

65.9 

66.9 

87.0 

63.0 

107.4 

43.4 

85.5 
95.8 

59.2 

67.5 

89.8 

-57.5 

62.7 

85.1 

-21.6 

51.0 

55.4 

-17.1 

55,0 

57.6 

10.2 

38.9 

40.2 

-81.8 

116.6 

142.4 

-29.1 

117.3 

120.8 

14.0 

134.3 

135.0 

-6.6 

110.7 

110.9 

-124.4 

195.5 

231.7 

-51.2 

48.8 

70.7 

-79.4 

35.5 

87.0 

-189.0 

65.2 

200.0 

-81.2 

38.0 

89.7 

1.4 

50.2 

50.2 

-29.4 

50.1 

58.1 

-37.6 

47.6 

60.7 

10.6 

48.8 

50.0 

95.0 

115.8 

149.7 

-42.0 

168.3 

173.5 

-29.4 

112.2 

116.0 

-22.3 

134.0 

135.9 

18.1 

159.7 

160.7 

68.8 

29.7 

75.0 

-75.5 

70.4 

103.3 

-70.8 

66.6 

97.2 

-66.7 

63.2 

91.8 

-57.8 

62.4 

85.1 

Cn.  CO  CO 

d  CO 
*-  u)  in 

-11.0 

53.8 

54.9 

8.5 

39.0 

39.9 

-27.8 

95.3 

99.3 

-42.6 

101.8 

110.3 

1.6 

133.2 

133.2 

-15.7 

113.8 

114.9 

-124.0 

192.0 

228.6 

-32.0 

45.3 

55.5 

-100.5 

34.8 

106.4 

-222.8 

67.1 

232.7 

-100.2 

31.7 

105.1 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

Mean 

Std 

Rms 

16545 

16229 

17470 

18325 

00 

10 

12679 

18140 

16469 

16737 

17301 

13055 

17895 

CM 

00 

CM 

166 

169 

175 

179 

207 

212 

218 

223 

228 

267 

270 

234 

244 

166 

169 

175 

179 

207 

212 

00 

CM 

223 

00 

CM 

CM 

267 

270 

CO 

CM 

244 

1989 

1989 

1989 

1989 

1989 

1989 

1989 

1989 

1989 

1989 

1989 

1989 

6861 

3182 

3183 

4137 

4138 

1184 

1186 

00 

00 

1189 

1190 

1195 

1196 

1198 

1199 

CO 

in 

O) 

I 

CM 

CO 

u 


CO 

LO 

lO 

0) 

O) 

CM 

CM 

CO 

CO 

O 

O 

U) 

CT) 

CM 

CO 

O 


109 


-0.200 

0.120 

0.234 

-0.272 

0.120 

0.210 

-0.063 

0.094 

0.113 

-0.184 

0.159 

0.243 

-0.003 

0.147 

0.147 

-0.159 

0.097 

0.186 

-0.218 

0.067 

0.228 

0.234 

0.134 

0.269 

0.339 

0.085 

0.350 

-0.172 

0.146 

0.226 

-0.009 

0.112 

0.112 

-0.188 

0.103 

0.215 

-0.044 

0.110 

0.119 

0.150 

0.129 

0.198 

0.094 

0.189 

0.211 

-0.015 

0.142 

0.142 

0.194 

0.216 

0.291 

-0.014 

0.078 

0.079 

-0.091 

0.252 

0.268 

0.054 

0.092 

0.107 

-0.019 

0.093 

0.095 

0.257 

0.151 

0.298 

0.109 

0.070 

0.129 

0.054 

0.130 

0.141 

0.079 

0.097 

0.125 

-0.004 

0.124 

0.124 

-23.8 

38.7 

45.4 

-85.7 

92.0 

125.7 

31.6 

57.8 

65.9 

-4.1 

59.9 

60.0 

-80.7 

51.6 

95.8 

-1.5 

89.4 

89.4 

11.6 

38.2 

39.9 

-0.1 

69.9 

69.9 

-4.2 

58.7 

59.0 

-21.5 

51.5 

55.8 

-69.0 

84.2 

108.8 

-50.7 

40.4 

64.8 

-41.0 

72,2 

83.0 

-28.2 

49.0 

56.6 

-83.8 

91.8 

124.3 

44.2 

57.8 

72.8 

-3.1 

58.6 

58.7 

-70.9 

52.4 

88.2 

-3.9 

85.2 

85.3 

-0.5 

36.5 

36.5 

17.9 

72.5 

74,7 

0.4 

69.3 

69.3 

-12.1 

59.1 

60.3 

-58.7 

86.5 

104.5 

-35.3 

41.0 

54.1 

-44.5 

72.2 

84.9 

-89.0 

58.1 

106.2 

3.7 

87.1 

87.2 

-54.8 

61.4 

82.3 

-54.6 

80.9 

97.6 

42.7 

69.0 

81.1 

-56.8 

91.6 

107.8 

-77.5 

43.8 

89.0 

134.7 

81.8 

157.5 

159.2 

81.3 

178.8 

-95.0 

79.0 

123.6 

29.7 
94.2 

98.8 

-94.4 

63.1 

113.5 

7.2 

90.9 

91.2 

-16.1 

86.0 

87.5 

-111.0 

61.3 

126.8 

1.6 

51.8 

51.9 

-51.9 

39.0 

64.9 

-68.4 

51.9 

85.9 

-109.3 

76.5 

133.4 

-111.0 

35.1 

116.4 

-46.8 

68.4 

82.6 

-204.6 

51.5 

211.0 

-46.9 

25.4 

53.4 

-52.1 

56.3 

76.7 

-70.4 

39.9 

80.9 

-56.8 

61.0 

83.3 

76.3 

77.0 

108.4 

67.8 

51.8 
85.3 

-5.8 

47.9 

48.3 

72.1 

64.3 

96.7 

11.1 

33.3 

35.1 

-86.0 

106.2 

136.6 

6.5 

34.7 

35.3 

-17.6 

50.4 

53.4 

70.3 

72.8 

101.2 

45.1 
33.5 

56.2 

16.3 

67.1 

69.0 

29.5 

53.6 
61.2 

-16.5 

66.3 

68.3 

6.3 

73.8 

74.1 

-92.4 

65.2 

113.1 

1.3 

52.5 

52.5 

-30.2 

40.7 

50.8 

-68.1 

51.8 

85.6 

-67.5 

68.7 

96.3 

111.3 

36.8 

117.2 

-45.8 

67.8 

81.8 

>26.4 

45.5 

>31.0 

-55.1 

24.3 

60.3 

■57.2 

56.4 

80.3 

•77.5 

40.3 

87.3 

54.5 
60.7 

81.6 

0)  "2  P 

2  55  ir 


c 

®  -P,  2 

®  2  P 
5  W  (T 


©  2  P 
^  (n  QC 


©■DC 

^  w  tr 


©■DC 

^  (O  cc 


c 

w  -p,  2 

©  2  P 

5 


c 

®  -p,  2 

®  :2  E 

S  (/)  DC 


c 

©  _  W 

®  :2  E 

^  CO  DC 


I  'B  £ 
:  CO  DC 


i  m 

:  B  £ 

CO  DC 


c 

2  -o 


:2  £ 
(O  DC 


©  CO 

©  2  P 

2  W  IT 


U) 

in 

in 

o 

00 

00 

in 

CM 

00 

o 

CO 

rv 

in 

o 

in 

(0 

00 

T— 

in 

(0 

CM 

CM 

CM 

CM 

o 

CM 


CO 

(O 

CM 


CO 

a> 

in 

0) 

T— 

CM 

00 

CO 

00 

o 

CO 

CM 

CO 

in 

CM 

00 

CM 

in 

in 

in 

CM 

CM 

CO 

CO 

in 

<n 

CO 

CM 

O 

IV 

in 

CO 

in 

t— 

O 

CO 

T- 

CO 

cn 

O) 

CM 

CO 

O 

O 

O 

o 

CO 

CO 

CO 

CO 

CO 

|v 

00 

CO 

in 

in 

CO 

CO 

in 

CM 

CM 

CM 

CM 

CM 

CM 

0) 

o 

CJ) 

(5> 

CJ) 

o 

00 

00 

00 

00 

00 

00 

cyj 

a 

0) 

O 

(J> 

cn 

r- 

r— 

r- 

Oi 

CM 

in 

CO 

|V 

cn 

CD 

in 

in 

in 

in 

CO 

r— 

O 

O 

O 

O 

CM 

CM 

CM 

CM 

CM 

00 

CM 

CM 

CO 

o 

in 

'Cf 

in 

o 

o 

o 

CM 

CO 

CO 

CO 

CO 

CO 

CO 

cn 

0) 

CD 

cn 

cn 

O) 

CJ) 

00 

00 

00 

00 

00 

00 

00 

(j) 

<n 

<n 

a 

0) 

r— 

r— 

CM 

in 

O 

CM 

CO 

r— 

00 

00 

o 

fv 

CO 

CO 

5 

IT) 

<J) 

CM 

CO 


O 


in 

in 

Of) 

o 

CM 

CM 

CO 

CO 

o 

o 

ir> 

9 


CM 

CO 

o 


no 


CO  «-  CO 

CM  <5>  CD 

q  q  p 

00  0)  CO  CD 

CO  lo  CD  CO  r** 

o  o  q  q  o  q 

O  «-  CD 

CO  If)  If) 
p  p  p 

CO  ^ 

CO  CD  CD 

r-  q 

00  00  |v 

IV  ^  CD 

q  q  q 

OO  CO  CO 
O  CD  CD 

q  q 

d  6  CD 

d  d  d  d  d  d 

1 

CD  d  d 

d  d  d 

1 

d  d  d 

d  d  d 

1 

0)0  0) 

CM  0)  CO 

CM  CM  q 

CO  00  CD  0)  CO  00 
tf)  tv  f-  If)  CO  00 

q  q  q  q  q 

|V  CO  O 

CM  O  ^ 
r-  q  q 

^  ID  00 

CM  ^ 

q  q 

^  O  CD 
o  CO  CO 

O  CO  CM 
Tt  CO  CD 
CM  q  q 

CD  d  d 

d  d  d  d  d  d 

odd 

d  d  d 

odd 

d  d  d 

CD  q  u) 

00  q  q  '5*; 

q  q  q 

q  q 

1-;  q  q 

r-  q  Tj- 

00  d 
t-  o  o 

d  oo’ d  CO  d 

to  ^  |v  CO  CO  |v 

d  cd  CD 

CM  iv 

1 

1—  «—  CM 
r-  O  O 

'  r—  T“ 

d  od  00* 

r-  rv 

d  »- 

f-  CO  CO 

r*«»  f—  ^ 

q  q  q  q  q  q 

q  q  q 

q  q  f-; 

q  q  q 

q  q  q 

d  d 
*  o  o 

T—  r- 

If)  cd  00*  ^ 

CO  ^  If)  If)  ''t  fv 

00*  od  d 

r-  |V  00 

CD  CM  ^ 

O  O 
'  ^  ^ 

O*  CD  cd 
■  <D  CO 

cd  cd  ^ 

r-  CO  CO 

q  CM  q 

q  q  q  q  q  q 

q  q  q 

q  q 

«-  q  q 

q  q  ^ 

cd  CO  ^ 

«-  O  O 

id  CD  If)  CM  00  00* 

CM  ^  If)  CD  ^  CD 

od  od  CD 
|v 

CM  |v  CM 

CO  O  *“ 

'  f-  r- 

^*  0)  00* 

^  0)  O 

T— 

CM  od  *- 
q  00  0) 

"(t  q 

(0  q  q  q  q  q 

t-;  q  q 

q  q  q 

Tt  q 

q  q  q 

U)  CD 
«-  00  00 

id  id  (D  If)  iv 
^  <0  |v  If)  CO  CO 

cd  d 
q  If)  CD 

CM  d 
r-  00  ^ 

CM  cd  ^ 
ID  U) 

d  ^ 

q  CO  '4' 

q  q  q 

q  q  q  q  q  p 

■«t  q  q 

q  q  q 

q  q  q 

q  q  q 

If)  o>  d 
't  00  o 

1  ^ 

od  00  id  cd  CM 

CM  ^  if)  ^  ^  CO 

od  id  od 
T-  ID  ID 

id  id  id 

'  IV  1^ 

0) 

q  ID  iv 

'It  CM  id 
f-  |v  <0 

q  q  q 

'?t  q  q  q  ^ 

q  q  q 

q  q  q 

q  q  q 

q  q  q 

CM  0) 
CO 

CM  ^  rv  id  iv  IV 
^  <D  |V  ID  CO  CD 

cd  id  d 

^  CD  00 

CM  od  •- 

CM  00  ID 

od  cd  cd 
ID  ID 

d  id  od 

CO  CO  CO 

1  ^  ^ 

-0.014 

0.061 

0.062 

ID  |v  CO 

O  CO  CD 

odd 

0.109 

0.106 

0.152 

-0.051 

0.081 

0.096 

-0.029 

0.122 

0.126 

0.076 

0.092 

0.120 

-0.004 

0.074 

0.074 

-0.029 

0.103 

0.107 

0.072 

0.206 

0.219 

0.102 

0.073 

0.126 

0.379 

0.166 

0.414 

0.378 

0.177 

0.418 

0.384 

0.150 

0.413 

0.361 

0.105 

0.376 

43.0 

50.0 

70.0 

54.5 

80.3 

97.1 

-27.6 

56.4 

62.8 

-0.8 

103.6 

103.6 

-26.3 

108.4 

111.5 

12.2 

89.5 

90.3 

q  q  q 
CM  od  od 

'  ID  LO 

52.1 
49.8 

72.1 

66.8 

80.3 

104.5 

-25.0 

66.8 

71.3 

-4.0 

102.9 

103.0 

-23.3 

108.5 

111.0 

13.4 

89.9 

90.9 

-2.5 

55.4 

55.5 

37.5 

58.9 

69.8 

-88.3 

100.2 

133.6 

75.6 

61.2 

97.3 

-20.9 

105.7 

107.7 

3.9 

122.1 

122.1 

32.9 

104.0 

109.1 

0.0 

58.6 

58.6 

38.2 

38.9 

54.5 

14.6 
77.2 

78.6 

-50.5 

41.5 

65.4 

-43.8 

74.6 

86.5 

-34.1 

97.6 

103.4 

q  q  ^ 

CD  CM  cd 
ID  CD  00 

-10.5 

64.2 

65.0 

-10.4 

47.2 

48.3 

34.0 

92.1 

98.0 

45.5 

44.1 

63.4 

178.5 

78.9 

195.2 

194.2 

95.5 

216,4 

218.5 

83.1 

233.8 

201.7 

55.3 

209.1 

40.9 

38.0 

55.8 

8*08 

ro8 

e'l7 

-58.9 

38.0 

70.1 

-57.5 

76.7 

95.9 

-62.3 

95.1 

113.7 

43.0 

64.5 

77.5 

-24.5 

63.2 

67.8 

If) 

9 

CN 

CO 

o 


CM 

ID 

ID 

9 

9 

CM 

CM 

CO 

CO 

O 

o 

0^4 

U) 

O 

04 

CO 

U 


111 


-0.056 

0.050 

0.075 

<0  0^  U)  (0  CM  CM 

00  CM  If)  in  in 

o  d  d  d  d  d 

0.084 

0.091 

0.124 

0.020 

0.174 

0.175 

-0.126 

0.084 

0.151 

0.025 

0.123 

0.126 

-0.061 

0.117 

0.132 

0.092 

0.109 

0.142 

0.104 

0.121 

0.159 

0.076 

0.099 

0.125 

0.008 

0.094 

0.095 

-0.060 

0.084 

0.103 

0.791 

0.077 

0.795 

o  in  o  o  1^ 

o  (?)  o  rv 

^  CM  in  q  q  q 

d  d  d  d  d  d 

1 

0.327 

0.183 

0.375 

0.224 

0.310 

0.382 

00  CM  O 
in  00 
q  r-  q 

d  d  d 

0.293 

0.245 

0.382 

0.475 

0.252 

0.538 

0.344 

0.201 

0.399 

0.315 

0.094 

0.328 

0.310 

0.078 

0.320 

0.361 

0.120 

0.381 

-0.323 

0.222 

0.392 

-76.3 

46.5 

89.4 

q  q  q  ^  CM  CM 

in  hs' d  d  iW 

CM  00  cn  1-  in  in 

1 

-40.4 

73.9 

84.2 

43.8 

103.1 

112.0 

-3.3 

57.3 

57.4 

25.7 

126.1 

128.7 

2.4 

101.4 

101.4 

17.4 

95.3 

96.9 

-31.8 

40.2 

51.2 

-50.4 

42.9 

66.2 

-6.7 

65.1 

65.4 

6.7 

99.2 

99.4 

-78.5 

46.7 

91.3 

r-*CMr^q^r- 
d  CM  f-‘  (d 
CM  00  00  CM  in  in 

-43.2 

73.2 

85.0 

41.9 

102.3 

1  10.5 

-2.0 

58.0 

58.1 

25.9 

126.3 

129.0 

VZOi 

Z'ZOl 

9*9 

23.4 

96.5 
99.3 

-28.2 

42.0 

50.6 

-54.6 

54.1 

76.9 

-5.1 

66.9 

67.1 

5.9 

99.6 

99.8 

46.0 

51.3 

68.9 

^  q  q  ^  q  q 

CO  cd  CO  d  rv  ^ 

'(t  0)  O  CO  00  O) 

*  ' 

76.7 

96.7 
123.4 

-35.3 

107.0 

112.7 

-56.9 

66.7 

87.7 

-4.1 

134.5 

134.6 

-33.9 

117.8 

122.5 

55.7 

107.9 

121.4 

48.1 

65.8 

81.5 

48.4 

54.3 

72.7 

8>Z 

9*0 

-13.8 

97.8 

98.8 

-84.1 

45.2 

95.4 

q  q  q  q  q  CM 

in  cd  d  d  d 

CO  to  ^  in  <0 

49.5 

64.3 

81.1 

25.9 

80.1 

84.2 

-69.0 

48.1 

84.1 

31.5 

106.0 

110.6 

-45.8 

94.9 

105.3 

-52.7 

94.3 

108.0 

CO  CM  f- 

in  d  CM 

CO  ^ 

-50.8 

41.6 

65.6 

-11.3 

66.6 

67.5 

51.2 

78.2 
93.5 

347.6 

34.8 

349.3 

115.9 

90.6 

147.1 

-3.2 

49.0 

49.1 

142.0 

77.2 

161.6 

42.4 

81.4 
91.8 

185.5 
78.9 

201.6 

151.5 

114.7 

190.0 

222.1 

140.1 

262.6 

142.0 

78.8 

162.4 

146.8 
38.9 

151.9 

139.9 
44.8 

146.9 

170.5 

52.4 

178.4 

-95.0 

74.0 

120.4 

-69.6 

44.9 

82.8 

q  ^  q  q  q  q 
d  d  CM  d  d 
^  <D  1^  ^  in  CO 

55.1 

64.4 

84.8 

42.9 

74.8 

86.2 

-34.6 

49.6 

60.5 

22.3 

105.2 

107.5 

-79.5 

90.9 

120.7 

-75.1 

87.6 

115.4 

-33.3 

38.4 

50.9 

-63.0 

41.4 

75.4 

-1  1.3 

64.9 

65.9 

71.0 

79.6 

106.6 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

CD 

W 

CD 

(/i 

CD 

W 

CD 

Vi 

(D 

Vi 

CD 

—  « 

CD 

Vi 

IS 

CO 

© 

CO 

© 

CO 

© 

(0 

© 

Vi 

© 

Vi 

03 

^  E 

(D 

^  E 

(D 

^  £ 

03 

'B  £ 

© 

E 

0) 

B  £ 

© 

:2  £ 

© 

B  £ 

© 

B  £ 

© 

”2 

£ 

© 

B  £ 

© 

-O 

£ 

© 

B  £ 

2 

W  DC 

cn  DC 

<n  DC 

<n  DC 

:e 

55 

DC 

to  DC 

:e 

to  DC 

(0  DC 

2 

(0  DC 

cc 

<0  DC 

S) 

DC 

2 

<0  DC 

00 

CO 

in 

CO 

CO 

CM 

CO 

00 

o 

CO 

t- 

00 

in 

CO 

03 

O 

CO 

03 

0) 

03 

o 

CO 

o 

00 

CO 

T— 

0) 

O 

03 

00 

in 

in 

0) 

00 

CO 

00 

CO 

CO 

03 

CO 

CO 

r- 

t— 

"T- 

»— 

r— 

r- 

r— 

0) 

CM 

in 

00 

O 

in 

CM 

in 

o 

o 

in 

in 

<o 

CO 

in 

in 

in 

in 

CO 

CO 

CO 

r— 

CM 

<J3 

CM 

»— 

in 

00 

o 

in 

CM 

in 

o 

o 

in 

in 

CO 

CO 

in 

in 

in 

in 

CO 

<0 

CO 

r— 

»— 

CM 

o 

o 

O 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

<J3 

(73 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

a 

(73 

03 

03 

03 

<73 

03 

03 

03 

03 

03 

03 

03 

r— 

«“ 

T— 

«— 

r— 

T— 

1^ 

00 

03 

f— 

CM 

CO 

in 

in 

Tt 

00 

C73 

03 

03 

03 

o 

o 

o 

o 

o 

o 

in 

in 

in 

03 

O 

o 

o 

CM 

in 

in 

in 

in 

in 

in 

in 

in 

CO 

CM 

CM 

CO 

CO 

in 

U3 

in 

in 

in 

o 

o 

o 

o 

O 

CM 

CM 

CM 

CM 

CM 

CO 

CO 

CO 

CO 

CO 

O 

o 

o 

o 

o 

112 


Table  6.  Project  MAGNET  Magnetic  Field  Flight  Statistics 
wrt.  WMM-90  (modified) 
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D32-451  4168  1993  295  295  18565  Mean  -23.3  -18.5  -69.1  10.9  -20.1  -17.4  2.597  -0.064 

Magnetic  Pole  Flight  Std  104.9  72.9  86.8  99.9  97.9  96.9  32.200  0.074 

Rms  107.4  75.2  110.9  100.4  99.9  98.5  32.300  0.098 
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Table  7.  WC-85  (updated)  Model  Coefficients  (Con.) 
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10.8 

-0.1 
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Table  8.  WMM-90  (modified)  Model  Coefficients 
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2 

1,690.2 

-375.1 

-0.1 

-9.5 

3 

0 

1,313.7 

0.0 

1.8 

0.0 

3 

1 

-2,243.6 

-285.4 

-7.0 

4.1 

3 

2 

1,244.8 

293.7 

-0.6 

1.9 

3 

3 

806.9 

-353.7 

-7.9 

-12.0 

4 

0 

935.5 

0.0 

0.8 

0.0 

4 

1 

784.0 

248.0 

0.7 

2.0 

4 

2 

324.0 

-237.6 

-6.7 

1.4 

4 

3 

-422.2 

89.7 

0.6 

2.8 

4 

4 

136.5 

-299.5 

-5.5 

-1.3 

5 

0 

-211.8 

0.0 

0.6 

0.0 

5 

1 

353.0 

44.2 

0.0 

0.2 

5 

2 

247.7 

152.7 

-1.1 

0.9 

5 

3 

-112.1 

-156.1 

-2.4 

0.6 

5 

4 

-162.3 

-68.7 

0.3 

1.9 

5 

5 

-38.2 

99.0 

2.5 

1.0 

6 

0 

62.5 

0.0 

1.2 

0.0 

6 

1 

63.7 

-14.7 

0.0 

0.3 

6 

2 

59.0 

82.6 

0.9 

-1.4 

6 

3 

-178.9 

70.5 

1.8 

-0.2 

6 

4 

2.3 

-53.3 

-0.3 

-0.6 

6 

5 

15.9 

1.1 

-0.2 

0.8 

6 

6 

-92.3 

22.2 

0.4 

2.0 

7 

0 

75.8 

0.0 

0.1 

0.0 

7 

1 

-63.2 

-80.5 

-0.9 

0.5 

121 


Table  8.  WMM-90  (modified)  Model  Coefficients  (Con.) 


n 

m 

jylTl 

fen 

um 

•  m 

gn 

•  m 

h„ 

7 

2 

0.8 

-25.7 

-0.6 

0.3 

7 

3 

29.6 

1.6 

0.6 

0.5 

7 

4 

2.2 

19.4 

0.9 

-0.4 

7 

5 

7.9 

16.4 

0.5 

0.0 

7 

6 

10.1 

-20.5 

0.1 

-0.1 

7 

7 

-1.3 

-8.3 

-0.7 

-0.7 

8 

0 

21.7 

0.0 

0.0 

0.0 

8 

1 

3.8 

10.7 

-0.4 

0.5 

8 

2 

-0.9 

-19.3 

-0.1 

0.0 

8 

3 

-10.2 

5.2 

0.2 

0.2 

8 

4 

-14.0 

-20.7 

-0.8 

0.5 

8 

5 

1.7 

11.3 

0.0 

-0.1 

8 

6 

2.6 

9.5 

0.0 

-1.3 

8 

7 

0.5 

-18.4 

-0.7 

-0.2 

8 

8 

-5.2 

-5.6 

-0.3 

-0.7 

9 

0 

3.6 

0.0 

0.0 

0.0 

9 

1 

9.5 

-21.9 

0.0 

0.0 

9 

2 

-0.9 

14.3 

0.0 

0.0 

9 

3 

-10.7 

9.5 

0.0 

0.0 

9 

4 

10.7 

-6.7 

0.0 

0.0 

9 

5 

-3.2 

-6.4 

0.0 

0.0 

9 

6 

-1.4 

9.1 

0.0 

0.0 

9 

7 

6.3 

8.9 

0.0 

0.0 

9 

8 

0.8 

-8.0 

0.0 

0.0 

9 

9 

-5.5 

2.1 

0.0 

0.0 

10 

0 

-3.3 

0.0 

0.0 

0.0 

10 

1 

-2.6 

2.6 

0.0 

0.0 

10 

2 

4.5 

1.2 

0.0 

0.0 

10 

3 

-5.6 

2.6 

0.0 

0.0 

10 

4 

-3.6 

5.7 

0.0 

0.0 

10 

5 

3.9 

-4.0 

0.0 

0.0 

10 

6 

3.2 

-0.4 

0.0 

0.0 

10 

7 

1.7 

-1.7 

0.0 

0.0 

10 

8 

3.0 

3.8 

0.0 

0.0 

10 

9 

3.7 

-0.8 

0.0 

0.0 

10 

10 

0.7 

-6.5 

0.0 

0.0 
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Table  8.  WMM-90  (modified)  Model  Coefficients  (Con.) 


n 

m 

„m 

fen 

h!r 

•  m 

gn 

•  m 

K 

11 

0 

1.3 

0.0 

0.0 

0.0 

11 

1 

-1.4 

0.0 

0.0 

0.0 

11 

2 

-2.5 

1.0 

0.0 

0.0 

11 

3 

3.2 

-1.6 

0.0 

0.0 

11 

4 

0.2 

-2.2 

0.0 

0.0 

11 

5 

-1.1 

1.1 

0.0 

0.0 

11 

6 

0.3 

-0.7 

0.0 

0.0 

11 

7 

-0.3 

-1.7 

0.0 

0.0 

11 

8 

0.9 

-1.5 

0.0 

0.0 

11 

9 

-1.1 

-1.3 

0.0 

0.0 

11 

10 

2.4 

-1.1 

0.0 

0.0 

11 

11 

3.0 

0.6 

0.0 

0.0 

12 

0 

-1.3 

0.0 

0.0 

0.0 

12 

1 

0.1 

0.7 

0.0 

0.0 

12 

2 

0.5 

0.7 

0.0 

0.0 

12 

3 

0.7 

1.3 

0.0 

0.0 

12 

4 

0.4 

-1.5 

0.0 

0.0 

12 

5' 

-0.2 

0.3 

0.0 

0.0 

12 

6 

-1.1 

0.2 

0.0 

0.0 

12 

7 

0.9 

-1.1 

0.0 

0.0 

12 

8 

-0.6 

1.2 

0.0 

0.0 

12 

9 

0.8 

-0.2 

0.0 

0.0 

12 

10 

0.2 

-1.3 

0.0 

0.0 

12 

11 

0.4 

0.6 

0.0 

0.0 

12 

12 

0.2 

0.6 

0.0 

0.0 
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Variation  (DV).  Otherwise  no  DV  corrections  were  made.  The  flights  are  generally  of  long 
range  over  remote  ocean  areas,  which  precludes  the  monitoring  of  DV.  Each  flight  typically 
lasts  10  to  12  hours  and  is  flown  at  an  average  speed  of  about  440  km/hr.  A  low-pass  filter  with 
a  cutoff  wavelength  of  7  km  is  routinely  applied  to  Project  MAGNET  survey  data  collected  at 
altitudes  greater  than  15,000  feet.  The  Honeywell  vector  magnetometer  is  calibrated  at  least 
once  a  year  at  the  National  Aeronautics  and  Space  Administration’s  (NASA's)  Coil  Room 
Facility  in  Greenbelt,  Maryland. 

These  Project  MAGNET  survey  data  sets  were  collected  primarily  in  ocean  areas  straddling  the 
geomagnetic  equator  with  the  intent  to  provide  absolute  vector-magnetic  measurements  in  the 
equatorial  region  to  counter  the  Backus  effect  (also  known  as  the  Perpendicular  Error  effect). 
This  effect  is  a  spherical-harmonic  modeling  error  that  generates  spurious  magnetic  anomalies 
along  the  geomagnetic  equator  if  only  scalar  Total-Magnetic-Intensity  measurements  are 
available  for  the  modeling.  Flights  were  also  made  to  the  North  and  South  Magnetic  Poles  in 
order  to  supply  data  to  pin  down  the  geomagnetic  pole  positions. 


2.2.3  POGSData 

The  POGS  satellite  was  launched  April  11,  1990.  The  salient  design  characteristics,  processing 
procedures,  and  calibrations  applied  are  given  by  Quinn  et  al.  (1993).  The  satellite  was  inserted 
into  a  near  circular,  polar  orbit  ranging  in  altitude  between  700  km  and  750  km.  The  orbit  was 
not  sun-synchronous  and  so,  covered  all  local  times.  The  fluxgate  vector  magnetometer  was 
located  at  the  tip  of  an  8-ft.  non-magnetic  Copper-Beryllium  boom  and  had  a  resolution  of  2  nT 
per  vector  axis.  Although  the  satellite  carried  a  vector  magnetometer,  there  was  essentially  no 
attitude  determination  system  on  the  satellite  to  orient  the  magnetometer.  We  say  essentially  no 
attitude  determination  system  because  a  coarse  attitude  device  was  constructed  from  9  of  the  52 
solar  panels,  which  provided  some  attitude  for  the  satellite  (but  not  the  boom-mounted 
magnetometer  relative  to  the  satellite)  with  respect  to  the  sun-line.  However,  the  accuracy  of 
this  system  was  at  best  only  ±  3”  with  degrading  accuracy  as  the  satellite  approached  the 
dawn-dusk  meridian.  So,  the  solar  panels  were  used  as  a  coarse  sun-sensor  which  provided  no 
data  on  the  night  side  of  the  Earth.  Consequently,  the  POGS  mission  was  essentially  a  scalar 
magnetic  survey  of  the  Earth.  The  attitude  data  were  nevertheless  made  available,  since 
applications  other  than  world  magnetic  modeling  may  find  it  useful. 

The  POGS  satellite  mission  carried  a  fluxgate  vector  magnetometer.  These  magnetometers  are 
not  absolute  instruments  and  so,  have  a  tendency  to  drift.  Extra  thermal  mass  was  added  to  the 
magnetometer  to  minimize  this  drift.  Even  so,  it  was  estimated  at  the  time  of  launch  that  the 
magnetometer  drift  could  be  as  much  as  50  nT/yr.  Additionally,  the  quartz  clock  on  the 
satellite,  to  which  the  magnetics  data  are  time  tagged  prior  to  being  telemetered  back  to  a  ground 
station,  also  has  a  tendency  to  drift.  If  not  properly  taken  into  account,  both  of  these  drifts  can 
be  misinterpreted  as  SV. 

The  POGS  quartz  clock  bias  and  drift  were  determined  by  periodically  monitoring  the  POGS 
clock  via  scheduled  timing  pulses  sent  from  the  satellite  to  the  Master  Ground  station  located  at 
NAVOCEANO,  which  was  equipped  with  a  cesium  standard  clock  calibrated  at  the  Naval 
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Observatory.  The  pulse  was  of  known  duration  (20  seconds).  The  difference  between  the 
scheduled  end  of  the  pulse,  which  was  known  a  priori,  and  the  actual  end  of  the  pulse, 
accounting  for  the  time  taken  to  travel  the  distance  between  the  satellite  and  the  ground  station, 
yielded  the  POGS  quartz-clock  timing  error,  which  could  be  measured  to  within  10  milliseconds. 
Repeating  this  procedure  at  roughly  10-day  interval  during  the  1991.0  to  1993.7  time  frame 
yielded  the  POGS  quartz  clock  bias  and  drift  rates  indicated  by  Quinn  et  al.  (1993).  The  time- 
corrected  magnetics  data  could  then  be  properly  merged  with  the  ephemeris  data  supplied  by 
DMA.  The  ephemeris  data  were  generated  from  Doppler-tracking  data  collected  from  DMA  s 
Tracking  Network  (TRANET),  which  monitored  the  emissions  from  two  beacons  of  different 
frequencies,  which  were  located  onboard  the  POGS  satellite.  The  satellite  ephemeris  was 
reconstructed  using  the  GEM-lOB  gravity  model.  The  resulting  satellite  positioning  Spherical 
Error  Probable  (SEP)  was  estimated  to  be  less  than  75  meters  by  DMA.  The  TRANET  stations 
were  dismantled  in  October  1993.  This  event  effectively  ended  the  POGS  mission  of  magnetic 
data  collection  for  world  magnetic  modeling  purposes. 

The  POGS  magnetometer  bias  and  drift  were  determined  by  selecting  the  Solar  quiet-time  POGS 
Total  Intensity  data  over  the  European  area  bounded  by  40**  N  and  60“  N  and  by  5  W  and  35  E, 
correcting  for  external  field  effects,  subsequently  removing  the  Total  Intensity  as  computed 
from  the  WMM-90  (modified)  model,  and  minimizing  the  residual  with  respect  to  the  bias  and 
drift  coefficients  of  a  linear  drift  model  via  least  squares.  The  main  assumption  in  this  procedure 
is  that  the  WMM-90  (modified)  model,  which  is  based  on  the  most  recent  available  magnetic 
observatory  annual-means  data  as  of  May  1994,  is  extremely  accurata  over  the  European  area 
due  to  the  substantial  clustering  of  Magnetic  Observatories  in  the  region.  During  the  course  of 
estimating  the  magnetometer  bias  and  drift  coefficients,  it  was  determined  that  for  a  short 
segment  of  data  between  1991.0  and  1991.4  the  timing  pulse  corrections  were  incorrectly 
applied.  So,  the  method  used  to  compute  the  scalar  magnetic  bias  and  drift  coefficients  for  the 
magnetometer  was  also  used  to  generate  an  additional  magnetic  correction  to  compensate  for  the 
timing-pulse  problem.  The  linear  drift  corrections  thus  determined  take  the  following  forms: 


a.  Magnetic  Compensation  for  Clock  Drift  Errors: 


ABc  =  be  +  be  (t  -  Xc) 
where  the  reference  time  in  years  is: 

Xc  =  1991.175394115 
and  where  the  bias  and  drift  rates  are: 

be  =  -18.783  nT 

be  =  144.338  nT/yr 


1991.0  <  t  <  1991.4  (29) 

(30) 

(31a) 

(31b) 
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b.  Magnetic  Correction  for  Magnetometer  Drift: 


ABm  -  bjvi  +  b(t  -  Tm) 

1991.0  <  t  <  1993.7 

(32) 

where  the  reference  time  in  years  is: 

Tm  =  1992.506748293 

(33) 

and  where  the  bias  and  drift  rates  are: 

bM  =  1 1.620  nT 

(34a) 

bM  =  -  16.167  nT/yr 

(34b) 

These  two  corrections  are  illustrated  in  figure  5.  When  these  corrections  are  applied  to  the 
POGS  data  (i.e.,  subtracted  from  the  POGS  Total  Intensity  data),  then  the  POGS  data  are 
considered  as  having  been  Ground  Truthed  and  are  thus  considered  to  be  absolute 
measurements.  It  should  be  noted  here  that  the  POGS  data  available  from  NGDC  do  not  include 
these  corrections.  The  POGS  10-day  file  statistics  with  these  drift  corrections  applied  are  given 
in  table  9.  They  have  been  computed  with  respect  to  the  WMM-90  (modified)  model. 

The  ionospheric/magnetospheric  field  contributions  to  the  POGS  data  were  considerably  larger 
than  the  magnetometer  drift  corrections  even  though  the  quiet-time  selection  criteria  previously 
given  in  the  Overview  was  employed.  Given  the  altitude  of  the  satellite  (about  725  km),  the 
ionospheric  portion  of  the  magnetic  field  measured  by  POGS  was  tentatively  attributed  to 
equatorial  and  mid-latitude  Spread-F  effects  (Kelley  [1989]).  This  is  a  phenomenon  occurring 
in  the  F-layer  of  the  ionosphere.  In  extremely  quiet  times,  this  field  was  limited  to  the 
equatorial  region.  During  these  times  the  phenomenon  is  referred  to  as  the  equatorial  Spread-F 
effect.  Frequently,  the  Spread-F  effect  was  prominent  both  at  mid-latitudes  and  at  high-latitudes 
and  exhibits  an  interesting  latitudinal  banding  effect.  The  number  and  magnetic  intensity  of 
these  bands  depends  on  the  amount  of  energy,  as  measured  by  the  Kp  index,  being  dumped  into 
the  ionosphere  as  a  consequence  of  solar  activity. 

Detailed  analyses  of  the  POGS  data  with  respect  to  ionospheric  and  magnetospheric  effects 
could  take  years.  For  the  WMM-95  modeling  effort  it  was  necessary  only  to  identify  and 
remove  fields  of  external  origin.  To  isolate  that  portion  of  the  observed  field  due  to  this 
presumed  Spread-F  effect,  magnetic  field  models  were  generated  from  the  POGS  data  at  the  14 
epochs  corresponding  to  the  14  groups  of  quiet-time  POGS  data  selected  for  modeling  use. 
The  groups  of  10-day  POGS  files  selected  for  each  model  are  listed  in  table  10.  Vector  data 
from  the  WMM-90  (modified)  model  were  inserted  in  an  equatorial  band  straddling  +  20“  about 
the  geomagnetic  equator  to  control  the  Backus  effect,  while  POGS  data  were  used  exclusively 
outside  of  this  band.  Except  for  equal  area  weighting,  no  weights  were  applied.  Internal  Gauss 
coefficients  up  to  degree  and  order  12  and  external  Gauss  coefficients  up  to  degree  and  order  5 
were  computed  at  each  of  the  14  epochs.  That  is,  14  models  were  generated  at  14  distinct 
epochs.  Scalar  Total  Intensity  field  values  computed  from  these  models,  using  both  internal  and 
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Table  9.  POGS  10-DAY  FILE  STATISTICS 
wrt.  WMM-90  (modified) 


|Dst| 

<  50 

Kp  <  2  + 

7  pm  < 

Local  Time  <  5  am 

File 

Year 

Day 

Begin 

Day 

End 

Records 

Bfv 

(nT) 

Used  in 
Model  #s 

pogslOI  1020.dst 

1991 

11 

20 

10210 

Mean 

Std 

15.0 

42.6 

1 

Rms 

45.1 

pogs1021  OSO.dst 

1991 

21 

30 

8491 

Mean 

Std 

19.2 

51.4 

1  &  2 

Rms 

54.8 

pogs1031 040.dst 

1991 

31 

40 

4262 

Mean 

Std 

14.0 

50.1 

2 

Rms 

52.0 

pogs1041050.dst 

1991 

41 

50 

5701 

Mean 

Std 

22.0 

49.1 

2  &  3 

Rms 

53.8 

pogs1051060.dst 

1991 

51 

60 

11914 

Mean 

Std 

20.8 

42.5 

3 

Rms 

47.3 

pogsl061 070.dst 

1991 

61 

70 

3572 

Mean 

Std 

11.3 

60.5 

Rms 

61.6 

pogsl071  OSO.dst 

1991 

71 

80 

5181 

Mean 

Std 

23.2 

59.3 

Rms 

63.7 

pogs1081090.dst 

1991 

81 

90 

2183 

Mean 

Std 

5.9 

64.2 

Rms 

64.4 

pogs1091  lOO.dst 

1991 

91 

100 

790 

Mean 

Std 

20.0 

55.8 

Rms 

59.2 

pogs1 1011  lO.dst 

1991 

101 

110 

56 

Mean 

Std 

-7.6 

18.5 

Rms 

19.9 

pogsl  1111  20.dst 

1991 

111 

120 

17 

Mean 

Std 

-2.3 

22.4 

Rms 

21.9 

pogsl  1211 30.dst 

1991 

121 

130 

3669 

Mean 

Std 

37.4 

37.6 

4 

Rms 

53.0 

pogsl  1311 40.dst 

1991 

131 

140 

8686 

Mean 

Std 

30.5 

37.5 

4 

Rms 

48.3 

pogsl  1411  50.dst 

1991 

141 

150 

1634 

Mean 

Std 

10.7 

41.5 

4 

Rms 

42.9 

pogsl  1511  eo.dst 

1991 

151 

160 

365 

Mean 

Std 

-27.5 

41.9 

4 

Rms 

50.1 

pogsl  1611 70.dst 

1991 

161 

170 

1977 

Mean 

Std 

12.1 

41.5 

4&5 

Rms 

43.2 

129 


pogsl  171 180. dst 

1991 

171 

180 

3938 

Mean 

24.5 

4  &  5 

Std 

43.0 

Rms 

49.5 

pogsl  181 190. dst 

1991 

181 

190 

3157 

Mean 

28.3 

5 

Std 

39.4 

Rms 

48.5 

pogsl  1 91 200  .dst 

1991 

191 

200 

0 

Mean 

0.0 

5 

Std 

0.0 

Rms 

0.0 

pogsl  201  210. dst 

1991 

201 

210 

3651 

Mean 

25.7 

5 

Std 

38.5 

Rms 

46.3 

pogsl  21 1  220.dst 

1991 

211 

220 

3966 

Mean 

33.9 

5 

Std 

31.0 

Rms 

45.9 

pogsl  221  230. dst 

1991 

221 

230 

1243 

Mean 

11.6 

5 

Std 

35.0 

Rms 

36.9 

pogsl  231 240.dst 

1991 

231 

240 

3079 

Mean 

24.3 

6 

Std 

37.0 

Rms 

44.2 

pogsl  241  250. dst 

1991 

241 

250 

2773 

Mean 

19.2 

6 

Std 

41.5 

Rms 

45.7 

pogsl  251  260. dst 

1991 

251 

260 

3483 

Mean 

29.6 

6 

Std 

36.7 

Rms 

47.2 

pogsl  261  270. dst 

1991 

261 

270 

8956 

Mean 

29.7 

6 

Std 

35.0 

Rms 

45.9 

pogsl  271  280.dst 

1991 

271 

280 

0 

Mean 

0.0 

Std 

0.0 

Rms 

0.0 

pogsl  281  290.dst 

1991 

281 

290 

5 

Mean 

60.9 

Std 

6.6 

Rms 

61.2 

pogsl  291  300. dst 

1991 

291 

300 

10 

Mean 

66.2 

Std 

18.9 

Rms 

68.5 

pogsl  301  31 0. dst 

1991 

301 

310 

296 

Mean 

8.5 

Std 

39.3 

Rms 

40.1 

pogsl  31 1 320.dst 

1991 

311 

320 

1208 

Mean 

39.4 

Std 

41.6 

Rms 

57.3 

pogsl  321 330.dst 

1991 

321 

330 

1630 

Mean 

13.6 
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43.4 

7 

Rms 

43.7 
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Std 
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Table  10.  WORLD  MAGNETIC  MODELS  BASED  ON  POGS  DATA 
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external  parts,  were  then  subtracted  from  the  Total  Intensity  data  of  the  10-day  POGS  files, 
creating  Total  Intensity  Residual  files  for  each  10-day  period.  Each  10-day  Residual  file  was 
averaged  over  all  Geomagnetic  Longitudes  for  each  1®  Geomagnetic  Latitude  band,  yielding  an 
estimated  external  field  correction  profile  for  each  10-day  file  as  a  function  of  Geomagnetic 
Latitude. 

Figures  6  through  19  exhibit  sample  external  field  correction  profiles  corresponding  to  specific 
10-day  POGS  data  files.  One  file  has  been  selected  from  each  of  the  14  modeling  groups  listed 
in  table  10.  They  span  the  1991.0  to  1993.7  time  frame.  Rotating  these  profiles  dirough  360®  of 
longitude  and  transforming  back  to  geodetic  coordinates  yields  charts  (47)  through  (60),  which 
display  a  global,  10-day  average  picture  of  the  mainly  Ionospheric  corrections.  This  procedure 
may  also  remove  residual  magnetospheric  effects  as  well.  These  figures  and  charts  exhibit  the 
high  degree  of  temporal  variability  that  is  clearly  external  in  origin. 

Magnetic  fields  of  crastal  origin  are  not  resolved  in  the  POGS  data  primarily  because  of  the 
satellite's  altitude  which  is  nearly  twice  that  of  the  1979-1980  MAGSAT  satellite  mission.  In 
unusually  quite  periods,  as  exemplified  by  charts  47, 48,  49,  54,  59,  and  60,  substantial  magnetic 
activity  associated  with  Field-Aligned  currents  can  still  be  seen  in  the  auroral  zones,  while 
pronounced  equatorial  Spread-F  effects  are  clearly  defined  along  the  geomagnetic  equator.  The 
POGS  satellite's  altitude  (about  725  km)  also  precludes  the  possibility  of  interpreting  the 
observed  external  equatorial  field  as  being  due  to  the  Equatorial  Electrojet  (EE).  This  is  in 
contrast  to  the  MAGSAT  mission  (about  400  km),  which  did  encounter  EE  effects  (Langel  et  al., 
1993b). 

Charts  47  and  54  make  an  interesting  contrast.  The  external  fields  are  generally  quiet  for  both, 
but  the  activity  along  the  equatorial  zone  has  moved  from  the  southern  geomagnetic  hemisphere 
to  the  northern  geomagnetic  hemisphere.  This  effect  is  apparently  associated  with  a  change  in 
direction  of  the  Interplanetary  Magnetic  Field  (IMF)  from  its  usual  southward  direction  to  a 
northerly  one.  Some  more  disturbed  10-day  periods  exemplified  by  charts  50  and  51  exhibit 
substantial  magnetic  activity  at  both  mid-latitudes  and  high-latitudes.  Charts  52  and  53  and 
charts  55  through  58  correspond  to  somewhat  less  disturbed  10-day  periods  and  prominently 
illustrate  varying  degrees  of  low-latitude  and  mid-latitude  banding. 


2.3  Weight  Factors 

The  modeling  procedure  employed  was  a  modification  of  that  formulated  by  Cain  et  al.  (1967). 
The  objective  was  to  minimize  the  chi-square  (x^)  function: 

=  Xr  +  Xe  +  Xv  + 

with  respect  to  168  MF  (i.e.,  internal)  Gauss  coefficients  corresponding  to  a  degree  and  order  12 
spherical-harmonic  model,  35  external-field  Gauss  coefficients  corresponding  to  a  degree  and 
order  5  spherical  harmonic  model  and  3  additional  External-Gauss  coefficients  corresponding  to 
a  degree  and  order  1  spherical-harmonic  model,  which  accounts  for  residual  Ring-current 
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IONOSPHERIC  CORRECTION 


Figure  6.  POGS  Total  Intensity  Ionospheric-Field  Correction:  1991,  Days  Oil  -  020 
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Ionospheric-Field  Correction:  1991,  Days  031  -  040 
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Figure  8.  POGS  Total  Intensity  Ionospheric-Field  Correction:  1991,  Days  051  -  060 
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Figure  9.  POGS  Total  Intensity  Ionospheric-Field  Correction:  1991,  Days  131  -  140 
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Ionospheric-Field  Correction;  1991,  Days  211  -  220 
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Figure  11.  POGS  Total  Intensity  Ionospheric-Field  Correction;  1991,  Days  261  -  270 
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Figure  12.  POGS  Total  Intensity  Ionospheric-Field  Correction:  1991,  Days  331  -  340 
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Figure  13.  POGS  Total  Intensity  Ionospheric-Field  Correction:  1992,  Days  121  -  130 
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Ionospheric-Field  Correction:  1992,  Days  171  - 180 
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Figure  15.  POGS  Total  Intensity  Ionospheric-Field  Correction:  1992,  Days  181  -  190 
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Intensity  Ionospheric-Field  Correction:  1992,  Days  241  -  250 
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Figure  17.  POGS  Total  Intensity  Ionospheric-Field  Correction:  1993,  Days  021  -  030 
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Intensity  Ionospheric-Field  Correction:  1993,  Days  141  - 150 
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Figure  19.  POGS  Total  Intensity  Ionospheric-Field  Con-ection:  1993,  Days  191  -  200 
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Chart  47.  Ionospheric  Magnetic  Field:  1991,  Days  01 1  -  020 
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Chart  48.  Ionospheric  Magnetic  Field:  1991,  Days  031  -  040 
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Chart  49.  Ionospheric  Magnetic  Field:  1991,  Days  051  -  060 
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Chart  50.  Ionospheric  Magnetic  Field:  1991,  Days  131-140 
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Chart  52.  Ionospheric  Magnetic  Field;  1991,  Days  261  -  270 
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Chart  53.  Ionospheric  Magnetic  Field;  1991,  Days  331  -  340 
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Chart  54.  Ionospheric  Magnetic  Field:  1992,  Days  121  -  130 
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Chart  55.  Ionospheric  Magnetic  Field:  1992,  Days  171  - 180 
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Chart  56.  Ionospheric  Magnetic  Field:  1992,  Days  181  -  190 
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Chart  58.  Ionospheric  Magnetic  Field:  1993,  Days  021  -  030 
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Chart  59.  Ionospheric  Magnetic  Field:  1993,  Days  141  - 150 
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Chart  60.  Ionospheric  Magnetic  Field:  1993,  Days  191  -  200 


effects,  as  monitored  by  the  Dst  index,  not  already  removed  via  the  subtraction  of  the  external 
fields  discussed  in  the  previous  subsection,  where: 

tr  =  (36a) 

1=1 

Xe  =  2  W0<(^0'-  “  ^0')^ 

1=1 

h 

1=1 

Xf  =  2  (^^''  “ 

f=l 

Uppercase  B's  refer  to  the  model  values  of  their  respective  magnetic  components,  while  the 
lower  case  b's  refer  to  the  observed  (measured)  values  of  their  respective  magnetic  components. 
The  subscript  i  refers  to  a  particular  data  point.  The  total  number  /  of  data  points  may  be 
different  for  each  magnetic  component.  Each  data  point  has  a  weight  w  which  depends  on 
several  factors: 


a.  Data  Type: 

Wt(POGS)  =  1  m=l 

w^(Project  MAGNET)  =  1  m  =  2 

Project  MAGNET  observatory  airswing  calibrations  yield  RMS  errors  on  the  order  of  35  nT, 
while  POGS  RMS  residual  errors  for  each  10-day  file  with  respect  to  the  WMM-90  (modified) 
model  average  approximately  45  nT.  These  two  numbers  are  indicators  of  the  relative  quality 
of  the  two  data  sets.  The  ratio  of  these  two  numbers  is  approximately  equal  to  one. 
Consequently,  the  two  data  sets  are  assigned  unit  weight  relative  to  each  other. 


b.  Equal  Area  Weighting: 

On  a  sphere  satellite  and  aircraft  magnetic  data  tend  to  be  more  concentrated  at  the  geographic 
poles  than  at  the  equator.  Since  there  is  no  reason  to  give  one  area  of  the  world  more  weight 
than  any  other,  a  colatitude  (6)  dependent  weight  factor  is  applied  to  each  data  point  such  that: 
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w^(0)  =  sin0 


(37) 


c.  Geomagnetic  Latitude  Weighting; 

If  only  scalar  Total  Intensity  data  are  used  for  modeling  purposes,  then  the  resulting  model  will 
generate  large  spurious  magnetic  anomalies  along  the  geomagnetic  equator  which  can  extend  to 
mid-latitudes.  On  the  other  hand,  it  is  desirable  to  use  as  much  scalar  Total  Intensity  data  and  as 
little  vector-magnetic  data  as  possible  in  the  modeling  process  because  the  vector  data  has 
additional  attitude  determination  errors  associated  with  it  which  tend  to  make  vector  data 
generally  less  accurate  than  the  scalar  Total  Intensity  data  computed  from  it  (Lowes  and  Martin 
[1987]).  Consequently,  in  order  to  maximize  the  accuracy  of  the  resulting  model.  Project 
MAGNET  vector-aeromagnetic  data  were  used  in  a  band  straddling  ±  20®  about  the  geomagnetic 
equator  to  counter  the  Backus  effect,  while  outside  of  this  band,  both  POGS  and  Project 
MAGNET  scalar  Total  Intensity  data  were  used.  However,  by  design.  Project  MAGNET 
surveys  were  concentrated  in  the  equatorial  band  where  vector  data  was  critically  needed. 
Therefore,  very  little  Project  MAGNET  data  were  collected  outside  of  this  band,  the  exception 
being  occasional  excursions  to  the  geomagnetic  poles.  Taking  0^  as  the  geomagnetic  latitude, 
this  weight  factor  takes  the  following  form: 


I0A/I  <  20® 


|0mI  >  20® 


>Va/(0m)  =  1 

for  k  =  l,2,3(r,0,(p)  1  ^  _  j 

(38a) 

=  0 

for  ^  =  4(F)  J  ’ 

=  0 

for  k  =  l,2,3(r,0,(p)  1.^  =  2 

(38b) 

=  1 

for  k  =  4(F)  J 

d.  Relative  Weighting  Among  Project  MAGNET  Flights  and  POGS  10-day  Files: 

RMS  statistics  were  computed  for  the  r,  6,  (p,  and  F  magnetic  components  of  Project  MAGNET 
data  on  a  flight  by  flight  basis.  Statistics  were  also  computed  for  the  F  magnetic  component  of 
the  POGS  data  on  a  10-day  file  by  10-day  file  basis.  These  statistics  were  computed  relative  to 
the  WMM-90  (modified)  model.  Relative  weights  were  assigned  for  each  magnetic  component 
using  the  squared  ratio  of  the  average  RMS  error  for  a  given  magnetic  component  for  all  flights 
(10-day  files)  used  in  a  model  to  that  of  a  particular  flight  (10-day  file)  for  the  same  magnetic 
component.  If  k  is  taken  as  the  index  for  a  particular  flight  (10-day  file)  and  if  I  is  taken  as  the 
index  representing  a  particular  magnetic  component  (r,0,(p,  F),  then  : 

x-A  =  (39) 
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e.  Outlier  Weighting: 

Editing  die  quiet-time  POGS  data  was  not  a  straightforward  matter  due  to  the  fact  that  the 
Spread-F  effect  was  broadly  present  at  essentially  all  latitudes  and  tended  to  blend  in  with  the 
Field-Aligned  current  effects  in  the  auroral  regions.  Consequently,  there  were  no  localized 
features  indicative  of  solar-induced  ionospheric/magnetospheric  activity  that  could  be  easily 
removed  from  the  data  through  interactive  computer  graphic  editing  or  through  statistical 
winnowing  without  a  substantial  loss  of  data.  Consequently,  the  average  Spread-F  and 
Field-Aligned  current  effects  were  removed  via  preliminary  modeling  of  the  POGS  data  as 
previously  indicated.  However,  some  of  the  external  field  activity  could  deviate  substantially 
from  the  10-day  average.  Therefore,  these  outliers,  which  would  primarily  occur  in  the  auroral 
zones,  were  down-weighted  via  the  following  weight  factor: 


where  AB  is  the  Total  Intensity  residual  of  the  data  point,  and  a  is  the  RMS  error  of  its 
corresponding  10-day  file  with  respect  to  the  WMM-90  (modified)  model. 

Project  MAGNET  vector  and  scalar  data,  having  been  collected  much  nearer  to  the  Earth's 
smface,  were  subject  to  magnetic  contamination  from  crustal  sources  as  well  as  from 
Ionospheric  sources  such  as  the  Sq  currents  and  the  EE  currents,  both  of  which  are  located 
between  100-km  and  200-km  altitude.  Outliers  in  Project  MAGNET  survey  data  due  to  these 
sources  were  down-weighted  in  the  same  manner  as  POGS,  data  using  eq.  (40),  with  AB 
representing  the  vector  or  scalar  residuals  with  respect  to  the  WMM-90  (modified)  model  and  a 
representing  the  RMS  statistics  for  a  particular  survey  flight  and  magnetic  component.  This 
weighting  is  in  lieu  of  applying  along-track,  low-pass  filters  with  750-km  wavelength  cutoffs  to 
the  aeromagnetic  data.  Such  filtering  also  results  in  a  substantial  loss  of  data. 

Since  the  Project  MAGNET  data  and  the  POGS  data  were  essentially  collected  concurrently  in 
time,  there  was  no  age-dependent  weight  factor  applied  to  the  data  in  contrast  to  the  1990  Epoch 
modeling  procedures.  The  total  weight  applied  to  a  specified  data  point  is  the  product  of  the 
foregoing  weight  factors,  so  that: 


Wiklmn 


WTm'WMni^Mi) 


where  the  indices  correspond  to  the  following: 


ith  -  data  point 

kth  -  aircraft  flight  or  satellite  10-day  file 

Ith  -  magnetic  component  (r,0,(p,F) 

mth  -  data  type  (Project  MAGNET  or  POGS) 

nth  -  geomagnetic  latitude  band  (equatorial  or  nonequatorial) 


(41) 
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2.4  Inverse  Modeling  Mathematical  Details 


Although  the  modeling  objective  is  to  determine  the  spherical-harmonic  Gauss  coefficients 
associated  with  the  magnetic  field  generated  in  the  Earth's  core  as  indicated  in  eqs.  (7a),  (7b), 
and  (7c),  the  satellite  magnetic  field  measurements  consist  of  magnetic  fields  of  core, 
ionosphere,  and  magnetosphere  origins.  Since,  unlike  the  aircraft  data,  these  data  are  temporally 
coherent,  it  is  possible  to  spherically  harmonic  model  the  external  fields  measured  by  POGS  as 
well  as  the  internal  fields.  In  so  doing,  we  minimize  the  possibility  of  contaminating  the 
internal  Gauss  coefficients  with  noise  of  external  origin.  Thus,  the  aircraft  data  are  used  to 
compute  only  the  168  internal  Gauss  coefficients  to  degree  and  order  12,  while  the  POGS  data 
are  used  to  compute  the  168  internal  Gauss  coefficients  to  degree  and  order  12,  the  35  external 
Gauss  coefficients  to  degree  and  order  5,  and  (following  Langel  and  Estes  [1985],  Quinn  et  al. 
[1986],  and  Langel  [1993])  the  3  external  Gauss  coefficients,  associated  with  the  so-called  Ring 
current  and  monitored  by  the  Dst  Index,  to  degree  and  order  1.  The  internal  Gauss  coefficients 
of  degree  n  and  order  m  at  time  t  are  denoted  SlS  g„Jt)  while  the  external  Gauss 

coefficients  of  degree  n  and  order  m  at  time  f  are  denoted  as  g„Jt)  snds„Jt).  The  external 
Gauss  coefficients  corresponding  to  the  Ring  current  for  degree  n  and  order  m  at  time  t  are 
denoted  as  a„Jt)  and  Using  this  notation,  the  potential  function  that  includes  all  three 

contributions  to  the  observed  field  has  the  following  form: 


T(r,0,(p,r)  =  a 


Nj  n 

XXIf 

n-l  m=0 


«+l 


[gn«(0  cos(w(p)  +  h„m{t)  sin(w9)] P”(0) 


-‘-EE  U  lqnm(t)  cos(/«cp)  +  sin(w(p)] P“(0)  (42) 

11  =  1  m=0  ^  ^ 

+  aDst(Q,t)  EE  U  [a„m(0  cos(w(p)  +  P„mSin(w(p)]P”(0) 

M  =  1  m=o'^  ' 


where  Nj  =  72,  N^.  =  5,  and  Nj^  —  1. 

Now,  defining  the  coefficient  set  {CJ  for  /=  1,2, ...,L  such  that; 


gnn; 

for  n(n+ 1)/2 

+ 

m 

< 

1  {n,  m) 

< 

Lz 

hnmy 

for  n(n— 1)/2 

+ 

m 

-1- 

L, 

< 

1  in,  m) 

< 

Lh 

Qnmy 

for  «(«  + 1)/2 

+ 

m 

+  Lh 

< 

1  (n,  m) 

< 

Lq 

Snm\ 

for  n (n-l)  12 

+ 

m 

+ 

Lq 

< 

l(n,m) 

< 

Ls 

^nm\ 

for  n{n  +  \)ll 

+ 

m 

Ls 

< 

1  (n,  m) 

< 

La 

for  n{n-  l)/2 

+ 

m 

La 

< 

1  (n,  m) 

< 

Ep 

(43) 


where: 
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A7(W’/+3) 

~  2 

(44a) 

Lh  =  Lg  + 

NiiNi+l) 

2 

(44b) 

Lq  =  Zz/i  4- 

Ne(Ne  +  V 

2 

(44c) 

Z/5  =  Lq  + 

NE{NE  +  i) 

2 

(44d) 

La  ~  Ls  + 

Nr{Ne+3) 

2 

(44e) 

Lp  =  La  + 

Nr(Nr  +  \) 

2 

(44f) 

it  is  possible  to  put  eq.  (42)  into  the  following  form: 

L 

V(r,  0,  (p,  t)  =  2  0, 9. 0  (45) 

1=1 

where  L-  L^  =  206. 

The  corresponding  geometric  functions  are: 

"  «  +  l 

a  j  cos(w  <p)  (0) ;  /  (n,  m)  =  +  m 

sin(/M(p)P2'(0) ;  l(n,m)  =  +  m  +  Lg 

a{^  cos(w(p)PS'{0) ; /(«,»i)  =  +  m  +  Lj, 

Qviir,Q,<?,t)  =  \  ^  [(46) 

af-jj  sin(w9)P“(0) ;  /  (n,m)  =  +  m  +  Lq 

aDsti^J)  cos(/n(p)PJ'(0) ;  /(n,w)  =  +  m  +  Ls 

aDstiQ,  t)  j  sin(7M  (p)P^(0) ;  / (n,  m)  =  +  m  +  Pa 
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Taking  the  negative  gradient  of  eq.  (45)  yields  the  magnetic  field  vector  5(r,0,(p,t),  with 
components: 


L 


Brir,  0,  9,  f)  =  Z 
/=1 

(47a) 

L 

5e(r,  0, 9, 0  =  Z 
/=1 

(47b) 

L 

B^(r,  0, 9, 0  =  Z  2q)/(^>  6. 9. 0 

1=1 

(47c) 

where: 

^  /  /X  ..  dQvi(r-A<?,t) 

Qr](r,d,(p,t)  -- 

(48a) 

Qdi(r,Q,(?,t)  -  -  r  3Q 

(48b) 

^  ^  ^  1  30  9.9.0 

e,p/(r,0,9,^)  -  -  ,sin9  atp 

(48c) 

The  Total  Intensity,  given  in  terms  of  these  three  components,  is: 

5f(r,0,9,t)  =  jBr(r,e,<p,ty  +50(r,  0,9,0^  +  ^<p(r,0,9,t)^ 

'  (49) 

It  is  clearly  nonlinear  in  terms  of  the  Gauss  coefficients  C,  .  Under  the  assumption  that  the 
Gauss  coefficients  vary  slowly  in  time,  eq.  (49)  may  be  linearized  via  Taylor  expansion. 

yielding: 


L 

hBpir,  0,  (p,  0  =  Z  Qpiir,  0,  (p,  t)  (50) 

/=i 

where  the  5  indicates  that  a  small  variation  of  the  indicated  quantity  is  to  be  taken,  and  where 
the  functions  Qp^r,d,(p,i),  defined  as: 


0F/(^0, 9.0 


g,.,(/-,e,(p,0  +5e(r,e,(p,0  g9,(r,e,<p,0  +£(p(r,e,c(),0  S,p/(r,9,(p,l) 


(51) 
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are  evaluated  at  a  fixed  epoch  t  =  t.  Thus,  eq.  (50)  may  be  integrated,  assuming  Qp(j',9,<p,t) 
varies  minimally  with  respect  to  small  changes  in  the  Gauss  coefficients,  to  yield: 

L 

B  F{f,  6,  (p,  0  =  ^^Ci(t)  Qpiif ,  0, 


which  is  a  useful  form  for  setting  up  the  least-squares  problem  to  determine  the  Gauss 
coefficients. 

Minimization  of  the  chi-square  function  of  eq.  (35)  requires  that: 


=  EII-sq  =  0 


which  is  satisfied  when: 


j  =  1,2, . .  .,206 


Therefore,  we  must  have: 


^  ^  ^ 

dCj  dCj  dCj  dCj 


j  =  1,2, . .  .,206 


which  is  a  nonlinear  system  of  206  equations  for  the  206  unknown  coefficients  in  the  set  {C,}. 
This  system  of  equations  is  nonlinear  because  depends  on  Bp  ,  which  in  turn  depends 
nonlinearly  on  the  Gauss  coefficients  through  eq.  (49)  which  depends  on  eq.  (51).  The 
individual  terms  in  eq.  (55)  are  determined  from  eqs.  (37a)  through  (37d),  eqs.  (47a)  through 
eqs.  (47c),  and  eq.  (52)  and  by  taking  the  indicated  derivatives  with  respect  to  the  Gauss 
coefficients,  so  that: 

^  2  L  Ir  Ir 

Cl  ^ri  ^ri  bri  Qrjifif9i,(pi,X)  (56a) 

^  1=1  i=l  i=l 


2  ^  /o 

=  X  X  0.-,  9^ OGjin,  0/, 9m t)  -  X  0/,  9m 'c)  (56b) 

^  /=1  Z-l  i=l 


^2  ^ 

^  M  QfS  P/? W(pi b(pi  01)  (p/, x)  (56c) 


1=1  i=l 
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L  Ip 


=  S  C/  X  ^FiiQriin,  Qi,  9/,  'C)  Qrjin,  e,-,  (?,,  X)  +  Qmiru  0/,  (?i,  'C)  Gejin,  0/,  9/,  T) 


1=1  i=l 


(56d) 


+  Q<fi{n,  0/,  9m  'c)  0M  9m  t:)}  -  X 

/=1 

where  eq.  (51)  has  been  used  to  obtain  eq.  (56d). 

Inserting  these  expressions  into  eq.  (55)  yields  the  following  system  of  equations: 

L 

=  j  =  1,2, . .  .,206 


where: 


J/y  =  ^yVriQrtMi,<^u'^)Qrj{ru^h^u<)  +  Wqi  Q9iiri,ei,(?i,x)  QBjiri,ei,(?i,'Z) 


j  9  ^ 

X^<p^eq./(^M0M9M'c)ecpAM0M9M'c)  +  T^WFi{Qri(ri,ei,(?i,x)Qrj(ri,ei,(?i,x)  (58) 


=  X  bri  Qrjiru  0i,  9^  'C)  +  X  >^9/^0/  Q^jin,  0M  9m  'c)  + 


-« <p  -r 

X  w<p'-  ^9'-  S<P7(''m  0m  9/.  t:)  +  X  '') 

1=1  i=l 

This  system  of  206  linear  algebraic  equations  can  be  written  in  matrix  form  as: 

cs  =  m 

which  may  be  inverted  to  yield: 

C  = 
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which  is  an  estimate  of  the  Gauss  coefficients  that  best  fit  the  observed  data  in  least-squares 
sense. 

Since  the  fourth  term  on  the  right-hand  side  of  eq.  (59)  depends  on  the  Gauss  coefficients,  it  is 
necessary  to  solve  eq.  (60)  iteratively.  If  p  is  the  iteration  index,  then  we  may  write  the  p'th 
iterative  solution  for  the  Gauss  coefficients  as; 

C(P)  =  p  =  1,2,  . . .,  pmax  (62) 

where  p„^  is  the  maximum  number  of  iterations,  and  where  the  WMM-90  (modified)  model 
was  used  as  the  initial  guess.  Note  that  the  2  matrix  in  eq.  (58)  is  independent  of  the  Gauss 
coefficients  and  so  must  be  computed  only  once.  One  set  of  Gauss  coefficients  was  computed 
via  the  above  technique  for  each  of  the  14  model  epochs  specified  in  table  10. 
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3.0  The  WMM-92.5  (optimum)  Model 

As  table  10  indicates,  14  World  Magnetic  Models  were  computed  at  distinct  epochs  covering  the 
time  frame  from  1991.0  to  about  1993.7.  These  models  are  given  the  designations  WMD-01, 
WMD-02  .  .  .  WMD-14.  The  internal  coefficients  are  derived  from  both  Project  MAGNET  and 
POGS  data,  while  the  external  coefficients  refer  only  to  the  POGS  data.  The  internal  Gauss 
coefficients  g^^^  corresponding  to  these  models  are  listed  in  table  1 1  and  are  plotted  as  functions 
of  time  in  figure  20.  The  internal  Gauss  coefficients  h^  corresponding  to  these  models  are 
given  in  table  12  and  are  plotted  as  a  functions  of  time  in  figure  21.  Also,  the  external  Gauss 
coefficients  and  s^  corresponding  to  these  models  are  given  in  tables  13  and  14, 
respectively,  while  the  external  Gauss  coefficients  and  (5^  ,  which  were  intended  to  extract 
any  residual  Ring-current  effects  not  already  accounted  for  by  the  direct  removal  of  the  10-day 
mean  ionospheric  correction  field,  a  procedure  that  was  primarily  intended  to  eliminate  Spread-F 
effects  and  Field- Aligned  current  effects,  are  listed  in  tables  15  and  16,  respectively.  These 
latter  sets  of  coefficients  are  quite  small,  indicating  that  most  of  the  magnetosphenc 
Ring-current  contribution  was  eliminated  as  part  of  the  mean  ionospheric  correction.  The  set  of 
external  Gauss  coefficients  q^  and  s^  is  asymmetric  in  the  sense  that  the  former  set  of 
coefficients  is  considerably  larger  than  the  latter.  It  is  not  clear  why  this  is  so,  although  it  is 
possible  that  the  portion  of  the  external  field  that  normally  contributes  to  the  s^  set  of 
coefficients  may  also  have  been  removed  by  the  10-day  mean  ionospheric  correction.  In 
addition,  the  q^^^  coefficient  seem  to  exhibit  a  quasi-annual  variation  as  one  would  expect. 

The  internal  Gauss  coefficients,  as  seen  in  figures  20  and  21,  clearly  exhibit  essentially  linear 
secular  behavior.  A  few  coefficients  exhibit  some  small  but  noticeable  external  field 
contamination  which  also  displays  a  quasi-annual  variation.  This  contamination  is  perhaps  most 
noticeable  in  the  gjo  coefficient,  although  in  general  the  contamination  is  quite  minimal.  This 
quasi-annual  external  field  contamination  is  detectable  primarily  as  a  consequence  of  the  POGS 
data  set  longevity  of  nearly  3  years.  The  effects  of  this  contamination  can  be  eliminated  via 
a  weighted,  linear  least-squares  fit  of  each  internal  Gauss  coefficient.  The  weights  come  from 
the  RMS  statistics  of  each  model  with  respect  to  all  POGS  and  Project  MAGNET  data  used  to 
produce  that  model.  Thus,  we  fit  the  14  values  of  a  single  Gauss  coefficient  to  a  linear  model  of 
the  form: 


C,(t)  =  C,(T)  -I-  C,  (t  -  T)  (63) 

where  the  reference  time  T  is  chosen  to  be  1992.5.  The  coefficients  CfT^  are  then  the  MF  values 

at  1992.5,  while  the  coefficients  Ci  represent  the  SV  coefficients  at  1992.5.  Performing  this 
least  squares  fit  on  all  of  the  internal  Gauss  coefficients  yields  the  WMM-92.5  (optimum)  model, 
the  coefficients  of  which  are  listed  in  table  17.  Using  the  SV  portion  of  the  WMM-92.5  model 
to  adjust  the  MF  coefficients  of  the  same  model  backward  2.5  years  yields  a  set  of  MF 
coefficients  at  the  1990.0  epoch,  which,  when  merged  with  the  WMM-92.5  (optimum)  SV 
coefficients,  yields  the  WMM-90  (revised)  model  listed  in  table  18. 
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Table  11.  Schmidt  Normalized  Gauss  Coefficients,  for  g,  Modeled  at  14  Epochs 
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Table  15.  Dst  Correction  Coefficients,  for  a,  Modeled  at  14  Epochs 
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FIGURE  20. 

GEOMAGNETIC  SPHERICAL-HARMONIC  COEFFICIENTS  g„„  VERSUS  TIME 

(units:  nT) 

(Pages  187-206) 
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FIGURE  21. 

GEOMAGNETIC  SPHERICAL-HARMONIC  COEFFICIENTS  h„^  VERSUS  TIME 

(units:  nT) 


(Pages  207-224) 
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Table  17.  WMM-92.5  (optimum)  Schmidt  Normalized  Gauss 

Coefficients 
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Table  17.  WMM-92.5  (optimum)  Schmidt  Normalized  Gauss 
Coefficients  (Con.) 
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Table  18.  WMM-90  (revised)  Schmidt  Normalized  Gauss  Coefficients  (Con.) 
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3.1  The  WMM-95  Model 


Truncating  the  SV  portion  of  the  WMM-92.5  (optimum)  model  to  degree  and  order  8  yields  the 
WMM-92.5(A)  SV  model  in  table  4.  This  SV  model  was  averaged  with  another  degree  and 
order  8  SV  model  (WMM-92.5(B)  SV),  which  was  derived  independently  by  the  BGS  and 
which  is  also  listed  in  table  4.  The  resulting  degree  and  order  8  SV  model  was  then  used  to 
adjust  the  full  degree  and  order  12  WMM-92.5  (optimum)  MF  model  coefficients  forward  2.5 
years  to  yield  MF  Gauss  coefficients  at  the  1995.0  epoch.  When  merged  with  the  Predictive 
1997.5  SV  model,  which  was  also  supplied  by  the  BGS,  the  desired  1995  Epoch  World 
Magnetic  Model  listed  in  table  3  results.  The  BGS  method  for  determining  the  WMM-92.5(B) 
SV  model  and  the  Predictive  1997.5  SV  model  is  discussed  by  Macmillan  (1994). 

An  indication  of  the  erratic  nature  of  the  geomagnetic  field  is  provided  by  the  wanderings  of 
the  North  Geomagnetic  Pole  and  the  South  Geomagnetic  Pole.  The  movements  of  these  poles, 
which  are  sometimes  referred  to  as  Dip  poles,  since  1945  are  illustrated  in  charts  61  and  62. 
These  charts  are  derived  from  the  Definitive  International  Geomagnetic  Reference  Field  (DGRF) 
models,  the  WMM-90  (revised)  model  and  the  WMM-95  model.  The  pole  movements  illustrate 
a  poorly  understood  phenomena  known  as  the  geomagnetic  jerk  which  occurred  around  1970.  It 
is  more  pronounced  for  the  South  Magnetic  Pole  movement,  which  exhibits  a  sudden  change  in 
direction  at  about  that  time.  These  jerks  occur  only  a  few  times  per  century  and  are  thought  to 
be  due  to  the  sudden  release  of  energy  built  up  over  a  long  period  at  the  core-mantle  boundary. 
This  build-up  is  presumed  to  arise  as  a  consequence  of  electromagnetic  coupling  between  the 
top  of  the  fluid  core  and  the  bottom  of  the  lower  mantle,  both  of  which  have  substantial 
electrical  conductivities.  These  jerks  have  been  correlated  to  irregular  changes  in  the 
Length-of-Day.  The  numerical  values  of  the  pole  positions  at  1-year  intervals  for  both  poles  are 
listed  in  table  19. 

It  should  be  noted  that  roughly  90  percent  of  the  Earth's  magnetic  field  is  contained  in  the 
degree-1  spherical-harmonic  coefficients:  g,®,  h,®,  and  h/.  These  three  coefficients  characterize 
the  Earth's  magnetic  dipole  field  and  form  the  basis  of  the  geomagnetic  coordinate  system, 
which  for  the  1995  epoch  is  illustrated  in  chart  63.  The  axis  of  the  geomagnetic  coordinate 
system  pierces  the  Earth's  surface  at  the  poles  of  the  Earth's  magnetic  dipole.  These  poles  are 
different  from  the  geomagnetic  dip  poles  which  are  derived  from  the  full  degree  and  order  12 
model.  The  location  of  the  dipole  poles  is  determined  when  the  horizontal  (H)  component, 
computed  from  just  the  degree- 1  spherical-harmonic  coefficients,  is  equal  to  zero.  The  dip  poles 
on  the  other  hand  are  determined  when  the  H  component  of  the  magnetic  field,  computed  from 
the  full  set  of  degree- 12  spherical-harmonic  coefficients,  is  equal  to  zero.  For  the  WMM-95 
model  at  1995.0,  the  North-magnetic  dipole  pole-position  is  located  at  79.30  degrees  latitude  and 
-71.46  degrees  longitude,  while  the  South-magnetic  dipole  pole-position  is  located  at  -79.30 
degrees  latitude  and  +108.54  degrees  longitude.  The  displacement  vector  for  the  eccentric 
(off-center)  dipole  for  1995.0  in  the  usual  Earth-fixed  spherical  coordinate  system  (i.e.,  Z-axis 
is  the  rotation  axis  (positive  North),  X-axis  points  to  the  prime  meridian,  and  the  Y-axis  is 
orthogonal  to  the  other  two,  thereby  creating  a  right-handed  system)  is  527.20  km  radially 
outward  from  the  Earth's  center,  with  a  colatitude  of  21.43  degrees  and  a  longitude  of  144.77 
degrees.  The  off-center  dipole  pole-positions  are  computed  from  the  degree  and  order  1  and  2 
spherical-harmonic  Gauss  coefficients. 
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Magnetic  Pole  Movement:  1945  -  2000 


62.  South  Magnetic  Pole  Movement:  1945  -  2000 
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Chart  63.  Geomagnetic  Coordinate  System:  1995.0  Epoch 


A  5-degree  grid  of  MF  and  SV  values  for  seven  components  (X,  Y,  Z,  H,  F,  D,  I)  of  the 
Earth's  magnetic  field  derived  from  WMM-95  are  tabulated  in  table  20.  Contours  of  five  of 
these  components  (Z,  H,  F,  D,  I)  for  the  MF  are  illustrated  in  charts  64  through  68,  while 
contours  of  the  SV  for  these  five  components  are  illustrated  in  charts  69  through  73.  These 
charts  are  plotted  on  a  corrected  Mercator  projection.  North  polar  stereographic  projections  of 
contours  for  these  same  five  components  plus  the  Grid  Variation  are  illustrated  in  charts  74 
through  79,  while  corresponding  SV  contours  of  these  components  for  the  north  pole  are 
illustrated  in  charts  80  through  85.  Similarly,  for  the  south  polar  region,  the  main  field  contours 
are  given  in  charts  86  through  91,  while  the  corresponding  Secular  Variation  contours  are 
illustrated  in  charts  92  through  97.  The  magnetic  field  components  contoured  on  both  the 
Mercator  and  the  polar  stereographic  charts  were  generated  with  respect  to  the  1984  World 
Geodetic  System  (WGS-84)  ellipsoid  at  the  1995.0  Epoch  for  the  MF  and  at  the  1997.5  Epoch 
for  the  SV  using  gridded  data  derived  from  the  WMM-95  model  and  the  GEOMAG  algorithm. 
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TABLE  20. 

\VMM-95  MAIN  FIELD  AND  ANNUAL  CHANGE  GRID  VALUES 
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NORTH  COMPONENT  (X)  WMM-95 
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EAST  COMPONENT  (Y)  WMM-95 
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9251.9  7295.7  4954.5  2313.7  -502.5  -3332.9  -6004.6  -8371.7  -10348.0  -11914.7  -13105.6  -13979.6 

-136.0  -161.2  -177.7  -184.3  -181.9  -172.3  -157.7  -140.0  -120.6  -100.6  -80.8  -62.5 
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16415.2  16415.2  16415.2  16415,2  16415.2  16415.2  16415.2  16415.2  16415.2  16415.2  16415,2  16415.2 

19.4  19.4  19.4  19.4  19.4  19.4  19.4  19,4  19.4  19.4  19.4  19.4 
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TOTAL  INTENSITY  (F)  WMM-95 


41388.5  40943.0  40444.1  39889.4  39282.9  38638.6  37976.9  37318.1  36679.0  36072.6  35508.6  34992.9 

-5.5  -3.3  -2.6  -3.5  -6.1  -9.7  -13.7  -17.6  -21.2  -24.2  -26.8  -28.7 


TOTAL  INTENSITY  (F)  WMM-95 


TOTAL  INTENSITY  (F)  WMM-95 


TOTAL  INTENSITY  (F)  WlVIM-95 


TOTAL  INTENSITY  (F)  WMM-95 


TOTAL  INTENSITY  (F)  WMM-95 


TOTAL  INTENSITY  (F)  WMM-95 


TOTAL  INTENSITY  (F)  \VMM-95 


297 


TOTAL  INTENSITY  (F)  \VMM-95 


TOTAL  INTENSITY  (F)  WMM-95 


DECLINATION  (D)  \VMM-95 


E.  Long 


DECLINATION  (D)  WMM-95 


100  105  110  115  E.  Long 


DECLINATION  (D) 


DECLINATION  (D)  WMM-95 


DECLINATION  (D)  WMM-95 


V©  ©S 

iri  ov 
<s 


VO  o\ 

V)  ov 


VO  rH  mo 


»  fS  ^  ^ 

Tf  VO  o  ri  It 


VO  OS 
lO  o\ 
lO  fS 


m  »H  m  fh 


VO  Ov 
d  ON 


VO  ON 
d  Ov 
On  r< 


VO  X 
ON  TT  »0  M 

fs,  ^  <s  fH 


«0  VO  X  ON 

m  vj  ON  m  VO  M 


m  ON 


^  ^  On 

no’  m  'V 

fr>  rH  m  ^ 


X  *<T 
X 

VJ  (S 


ON  V) 
VO  <s 


NO  fS 
TT 


'Nf 

On  d 


VO  On 

d  m 

VO  n 


m  m 

vi  d 

VO  n 


m  0\  1-1  t 


r»  *1  On  o 


NO  r'  ^  ®i 

On  m  d  1-N  On 


X  ©  ^ 

ri  ON  VO  d 
m  1-H 


t-.  VO 
X 

m  ^ 


»-i  r>  NO  VI 
d  On 


VO 

ON 


On  ^ 
ON  l> 


00  (S  ^  ^ 

00  d  i>  d 


X 


305 


DECLINATION  (D)  WMM-95 


E.  Long 


DECLINATION  (D)  WMM-95 


DECLINATION  (D)  WMM-95 


309 


52.5  147.5  142.5  137.5  132.5  127.5  122.5  117.5  112.5  107.5  102.5  97.5 

-0.5  -0.5  -0.5  -0.5  -0.5  -0.5  -0.5  -0.5  -0.5  -0.5  -0.5  -0.5 


DECLINATION  (D)  \VMM-95 


DECLINATION  (D)  WMM-95 


u 


O  O  >-11 


O  tT 


t-  VO  »H  1-H 


aovo  oo  vjo 
rif*i  fO>-?  vjo  |jn> 


00  On 


c*^  '•t 

r»  ' 


>1-1  On  TT 


O  »H  “IT 


VO  00 


i/i  fH  ®® 


1  VO  i-H  n 


o\  o 

00  o  ri  o 


tt  M  ov  t~:  fi 

ovd  00'^  voro  vifs  Mjjo 


I  1-H  o 


w  >-i  < 


ft  < 


'^1-3  1-N  o^  o 


nJ 


311 


INCLINATION  (I)  WMM-95 


nJ 

u 


00  in 

K  d 


d  d 


t--  m 

d 


in 

Qo  d  in  d  d  d 


ri  v© 
M  d 


o  i-H  < 


o  oo 


mm  isoo  yoTt  ooo 


O  1-1  i-H  r-1 


f'l  O  1-1  < 


312 


E.  Long 


INCLINATION  (I)  WMM-95 


nJ 

U 


vi  Ti  fo  ^ 


rn  iH  sJ  t  ’ 


1-5  ri  t  ’ 


>1-1  1-1 


VO 

3^ 


VO  fS 

3^ 


t^o 

Tfi*  vo'^  vdit  vo«r) 


r^o 

oorn  r4rn  iS*^  iS'^  i2*^ 


o«r)  if»ov  »-i''T 

ovri  vorr  votj-  vjin 


00  ®v 

3  vi  ^ 


3  - 


VO  00 

3® 


m  r< 

3® 


o  o 

3® 


IT)  n 
00  lo 


•rr  1-t  ID  m 


t~-  i-H  OV  1-1 


nJ 

U 


313 


100  105  110  115  E.  Long 


INCLINATION  (I)  WMM-95 


(  n  VO 


S  So  ^ 


vovo  1-Ho  f«^v«  mo  \OTf  r'Tr 

•^ri  or4  ^ori  i-!ri  'Ori  fhi-h  vdi-5  cSfH  mr4  i/ri  vom  t^’«T  ri-^  <s*^ 


m  <s 

3  ^ 


On  VI 
VO  m 


o  m 


in 

m  m 


o  m 

VO  m 


VO  m 
in  m 


hJ 

M 


314 


INCLINATION  (I)  WMM-9S 


INCLINATION  (I)  \VMM-95 


E.  Long  300  305  310  315  320  325  330  335  340  345  350  355  E.  Long 


INCLINATION  (I)  WMM-95 


Q\  tn 
wi 


V© 

3 


■Tt  so 

^  r4 


O  1-4  11  l-( 


o  vq 
On  ih 


7 


VO  VO 

O  fS 


vq 


i/j 

s  ^ 


t^vq 

mr4  irjrn  t'f*? 

VO  \0  \£j  \0 


•  H  n  *S 


3  ^ 


fn 

3  ^ 


ir>  ffi 

3  ^ 


On  On 

ri  o 


n 

3^‘ 


fS  iq 

«s 
tt  • 


o>  <s 


fM  V© 


n  On 


I©  IS 

3  ^  3  ^ 


iq  rq 

so 

^  ■ 


•<t  ^  ( 


IS  fS 

so  1/^ 

T  ■ 


3"? 


'o  fn 
vq  . 


■"tt  ON 
O  '«T 


ov  m 

3  ^ 


IS 

3  ^ 


t-  VO 

T  ‘ 


SO 


t^  ON 

tn  K 
vq  . 


IS  1^ 

5g 


l^  rs 
3  ^ 


t' 

3  ^ 


3  T 


318 


E.  Long  0  5  10  15  20  25  30  35  40  45  50  55 


INCLINATION  (I)  WMM-95 


320 


E.  Long  120  125  130  135  140 


INCLINATION  (I)  \VMM-95 


INCLINATION  (I)  WMM-95 


INCLINATION  a)  WMM-95 


o  o 
^  fs 


V]  l/J 
VV  ON 

T  ' 


V)  ^  »D  1 


Tf 

7  7 


^  'O 

7  7 


VO  <s 

o\  o 

M  M 


1-?  d 

n 


1-1  < 


■  ^  VO  r» 


00 

1/1  ff?  00  ri 

VI  F-t  V)  FN 


<S  1/1 

ON  vi 

7  7 


00 

vi  ve 

7  7 


F-(  r> 

7  7 


fH  o 
vd  vd 
7  7 


r» 
ri  vd 

7  7 


r;  w 

fS  00 
fF)  vX 


VO  *«r 

7  7 


ri  r^i 

7 


-Tf  o\ 

vd  d 

M  i-< 


VO  Tf 
fFl  00 


VO  VI  ON  VO 

d  fh  1-!  on 

VO  fH  VO  • 


O  VO  VO  o 

SS 


■  *-(  VI  1-H  VI  1 


7  7  ’ 


o  O  O  00  o  < 

■FT  vi  On  ri  fFi  < 

fn  — (  m  F-v  FT  1 


fF)  oo 

d  F# 

VO  I 


fF}  © 

0\  r4 

7  7 


00  ON 

7  ■ 


O  Ov 

7  ’ 


fF}  00 

7  • 


VI  M 

FH 

7  ' 


VO 

7  ■ 


VI  fFj 

vi  vd 

tF>  • 


© 

On  vi 

7  ’ 


00  n  On  ON 

FT  vd  FT 

7  .  FT  . 


fFl  V) 

7  ■ 


00  FH  fFl  00 

ON  ft  1^  rn 


Ov  FT  FT  fFl 


<5  FH  (^  © 


VI 

fFl  ^ 

7 


VI  VI 

d  d 

7 


VO 

ON  ri 


fFl  VO 


VI  fFl 


r«  N 

7 


323 


E.Long  300  305  310  315  320  325  330  335  340  345  350  355  E.  Long 


Vertical  Component  (Z) 


SOI  3  ,S9  3  „se  M  „Sl  M  „SS  M  „S6  M  „S£t  M  „S£I 


326 


Chart  65.  Horizontal  Intensity  (H) 


Total  Intensity  (F) 


1995.0  at  surface  of  WGS-84  reference  ellipsoid.  units:  degrees 


Chart  68.  Inclination  (I) 


units:  nT/y: 


Chart  70.  Horizontal  Intensity  (H) 
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Chart  79.  Grid  Variation  (GV) 
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Chart  90.  Inclination  (I) 


351 


180  W 


0  £ 

1995.0  at  surface  of  WGS-84  reference  ellipsoid.  Units:  degrees 


Chart  9 1 .  Grid  Variation  (GV) 


352 


reference  ellipsoid. 


Chart  92.  Vertical  Component  (Z) 
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Chart  93.  Horizontal  Intensity  (H) 
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Chart  96.  Inclination  (I) 
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INTERNAL  WMM-95  MODEL  COEFFICIENTS 
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o  o  n 


Q*********************************************************************** 


SUBROUTINE  GEOMAG  (GEOMAGNETIC  FIELD  COMPUTATION) 

C 

C 

Q*********************************************************************** 

C 

c 

C  WMM-95  is  a  Defense  Mapping  Agency  (DMA)  standard  product.  It  is 
C  covered  under  DMA  Military  Specification:  MIL‘W-89500  (1993). 

C  For  information  on  the  use  and  applicability  of  this  product  contact 
C 

C  DIRECTOR 

C  DEFENSE  MAPPING  AGENCY/HEADQUARTERS 

C  ATTN:  CODE  PR 

C  8613  LEE  HIGHWAY 

C  FAIRFAX,  VA  22031-2137 

C 
C 

Q*********************************************************************** 

C 

C  GEOMAG  PROGRAMMED  BY; 

C 

C  JOHN  M.  QUINN  7/19/90 

C  FLEET  PRODUCTS  DIVISION,  CODE  N342 

C  NAVAL  OCEANOGRAPHIC  OFFICE  (NAVOCEANO) 

C  STENNIS  SPACE  CENTER  (SSC),  MS  39522-5001 

C  PHONE:  COM:  (601)  688-5828 

C  AV:  485-5828 

C  FAX:  (601)  688-5521 


C*********************************************************************** 

C 

C 

C  PURPOSE:  THIS  ROUTINE  COMPUTES  THE  DECLINATION  (DEC), 

C  INCLINATION  (DIP),  TOTAL  INTENSITY  (TI)  AND 

C  GRID  VARIATION  (GV  -  POLAR  REGIONS  ONLY,  REFERENCED 

C  TO  GRID  NORTH  OF  POLAR  STEREOGRAPHIC  PROJECTION)  OF 

C  THE  EARTH'S  MAGNETIC  FIELD  IN  GEOTETIC  COORDINATES 

C  FROM  THE  COEFFICIENTS  OF  THE  CURRENT  OFFICIAL 

C  DEPARTMENT  OF  DEFENSE  (DOD)  SPHERICAL  HARMONIC  UORLD 

C  MAGNETIC  MODEL  (UHM-95).  THE  UMM  SERIES  OF  MODELS  IS 

C  UPDATED  EVERY  5  YEARS  ON  JANUARY  1'ST  OF  THOSE  YEARS 

C  WHICH  ARE  DIVISIBLE  BY  5  (I.E.  1980,  1985,  1990  ETC.) 

C  BY  THE  NAVAL  OCEANOGRAPHIC  OFFICE  IN  COOPERATION 

C  WITH  THE  BRITISH  GEOLOGICAL  SURVEY  (BGS).  THE  MODEL 

C  IS  BASED  ON  GEOMAGNETIC  SURVEY  MEASUREMENTS  FROM 

C  AIRCRAFT,  SATELLITE  AND  GEOMAGNETIC  OBSERVATORIES. 


Q*******************************ifr******Klf****************i»r*****r********** 

C 

C 

C  MODEL:  THE  WMM  SERIES  GEOMAGNETIC  MODELS  ARE  COMPOSED 

C  OF  TWO  PARTS;  THE  MAIN  FIELD  MODEL,  WHICH  IS 

C  VALID  AT  THE  BASE  EPOCH  OF  THE  CURRENT  MODEL  AND 

C  A  SECULAR  VARIATION  MODEL,  WHICH  ACCOUNTS  FOR  SLOW 

C  TEMPORAL  VARIATIONS  IN  THE  MAIN  GEOMAGNETIC  FIELD 

C  FROM  THE  BASE  EPOCH  TO  A  MAXIMUM  OF  5; YEARS  BEYOND 

C  THE  BASE  EPOCH.  FOR  EXAMPLE,  THE  BASE  EPOCH  OF 

C  THE  WHM-95  MODEL  IS  1995.0.  THIS  MODEL  IS  THEREFORE 

C  CONSIDERED  VALID  BETWEEN  1995.0  AND  2000.0.  THE 

C  COMPUTED  MAGNETIC  PARAMETERS  ARE  REFERENCED  TO  THE 

C  WGS-84  ELLIPSOID. 


Q*********************************************************************** 

C 

C 
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c 


ACCURACY:  IN  OCEAN  AREAS  AT  THE  EARTH'S  SURFACE  OVER  THE 
ENTIRE  5  YEAR  LIFE  OF  A  DEGREE  AND  ORDER  12 
SPHERICAL  HARMONIC  MODEL  SUCH  AS  WMM-95,  THE  ESTIMATED 
RMS  ERRORS  FOR  THE  VARIOUS  MAGENTIC  COMPONENTS  ARE: 

DEC  -  0.5  Degrees 

DIP  -  0.5  Degrees 

TI  -  280.0  nanoTeslas  (nT) 

GV  -  0.5  Degrees 

OTHER  MAGNETIC  COMPONENTS  THAT  CAN  BE  DERIVED  FROM 
THESE  FOUR  BY  SIMPLE  TRIGONOMETRIC  RELATIONS  WILL 
HAVE  THE  FOLLOWING  APPROXIMATE  ERRORS  OVER  OCEAN  AREAS: 

X  -  140  nT  (North) 

Y  -  140  nT  (East) 

Z  -  200  nT  (Vertical)  Positive  is  down 

H  -  200  nT  (Horizontal) 

OVER  LAND  THE  RMS  ERRORS  ARE  EXPECTED  TO  BE  SOMEWHAT 
HIGHER,  ALTHOUGH  THE  RMS  ERRORS  FOR  DEC,  DIP  AND  GV 
ARE  STILL  ESTIMATED  TO  BE  LESS  THAN  1.0  DEGREE,  FOR 
THE  ENTIRE  5-YEAR  LIFE  OF  THE  MODEL  AT  THE  EARTH'S 
SURFACE.  THE  OTHER  COMPONENT  ERRORS  OVER  LAND  ARE 
MORE  DIFFICULT  TO  ESTIMATE  AND  SO  ARE  NOT  GIVEN. 

THE  ACCURACY  AT  ANY  GIVEN  TIME  OF  ALL  FOUR 
GEOMAGNETIC  PARAMETERS  DEPENDS  ON  THE  GEOMAGNETIC 
LATITUDE.  THE  ERRORS  ARE  LEAST  AT  THE  ECaJATOR  AND 
GREATEST  AT  THE  MAGNETIC  POLES. 

IT  IS  VERY  IMPORTANT  TO  NOTE  THAT  A  DEGREE  AND 
ORDER  12  MODEL,  SUCH  AS  WMM-95,  DESCRIBES  ONLY 
THE  LONG  WAVELENGTH  SPATIAL  MAGNETIC  FLUCTUATIONS 
DUE  TO  EARTH'S  CORE.  NOT  INCLUDED  IN  THE  SERIES 
MODELS  ARE  INTERMEDIATE  AND  SHORT  WAVELENGTH 
SPATIAL  FLUCTUATIONS  OF  THE  GEOMAGNETIC  FIELD 
WHICH  ORIGINATE  IN  THE  EARTH'S  MANTLE  AND  CRUST. 
CONSEQUENTLY,  ISOLATED  ANGULAR  ERRORS  AT  VARIOUS 
POSITIONS  ON  THE  SURFACE  (PRIMARILY  OVER  LAND,  IN 
CONTINENTAL  MARGINS  AND  OVER  OCEANIC  SEAMOUNTS, 

RIDGES  AND  TRENCHES)  OF  SEVERAL  DEGREES  MAY  BE 
EXPECTED.  ALSO  NOT  INCLUDED  IN  THE  MODEL  ARE 
NONSECULAR  TEMPORAL  FLUCTUATIONS  OF  THE  GEOMAGNETIC 
FIELD  OF  MAGNETOSPHERIC  AND  IONOSPHERIC  ORIGIN. 

DURING  MAGNETIC  STORMS,  TEMPORAL  FLUCTUATIONS  CAN 
CAUSE  SUBSTANTIAL  DEVIATIONS  OF  THE  GEOMAGNETIC 
FIELD  FROM  MODEL  VALUES.  IN  ARCTIC  AND  ANTARCTIC 
REGIONS,  AS  WELL  AS  IN  EQUATORIAL  REGIONS,  DEVIATIONS 
FROM  MODEL  VALUES  ARE  BOTH  FREQUENT  AND  PERSISTENT. 

IF  THE  REQUIRED  DECLINATION  ACCURACY  IS  MORE 
STRINGENT  THAN  THE  WMM  SERIES  OF  MODELS  PROVIDE,  THEN 
THE  USER  IS  ADVISED  TO  REQUEST  SPECIAL  (REGIONAL  OR 
LOCAL)  SURVEYS  BE  PERFORMED  AND  MODELS  PREPARED  BY 
NAVOCEANO,  WHICH  OPERATES  THE  PROJECT  MAGNET 
AIRCRAFT  AND  THE  POLAR  ORBITING  GEOMAGNETIC  SURVEY 
(POGS)  SATELLITE.  REQUESTS  OF  THIS  NATURE  SHOULD 
BE  MADE  THROUGH  DMA  AT  THE  ADDRESS  ABOVE. 


Q*********************************************************************** 


C 

C 

C  USAGE:  THIS  ROUTINE  IS  BROKEN  UP  INTO  TUO  PARTS: 

C 

C  A)  AN  INITIALIZATION  MODULE,  WHICH  IS  CALLED  ONLY 

C  ONCE  AT  THE  BEGINNING  OF  THE  MAIN  (CALLING) 

C  PROGRAM 

C  B)  A  PROCESSING  MODULE,  WHICH  COMPUTES  THE  MAGNETIC 

C  FIELD  PARAMETERS  FOR  EACH  SPECIFIED  GEODETIC 
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POSITION  (ALTITUDE,  LATITUDE,  LONGITUDE)  AND  TIME 


INITIALIZATION  IS  HADE  VIA  A  SINGLE  CALL  TO  THE  MAIN 
ENTRY  POINT  (GEOMAG),  WHILE  SUBSEQUENT  PROCESSING 
CALLS  ARE  HADE  THROUGH  THE  SECOND  ENTRY  POINT  (GEOMG1). 
ONE  CALL  TO  THE  PROCESSING  MODULE  IS  REQUIRED  FOR  EACH 
POSITION  AND  TIME. 

THE  VARIABLE  MAXDEG  IN  THE  INITIALIZATION  CALL  IS  THE 
MAXIMUM  DEGREE  TO  WHICH  THE  SPHERICAL  HARMONIC  MODEL 
IS  TO  BE  COMPUTED.  IT  MUST  BE  SPECIFIED  BY  THE  USER 
IN  THE  CALLING  ROUTINE.  NORMALLY  IT  IS  12  BUT  IT  MAY 
BE  SET  LESS  THAN  12  TO  INCREASE  COMPUTATIONAL  SPEED  AT 
THE  EXPENSE  OF  REDUCED  ACCURACY. 

THE  PC  VERSION  OF  THIS  SUBROUTINE  MUST  BE  COMPILED 
WITH  A  FORTRAN  77  COMPATIBLE  COMPILER  SUCH  AS  THE 
MICROSOFT  OPTIMIZING  FORTRAN  COMPILER  VERSION  4.1 
OR  LATER. 


Q********************************************************************** 


REFERENCES: 

JOHN  M.  QUINN,  DAVID  J.  KERRIDGE  AND  DAVID  R.  BARRACLOUGH, 
WORLD  MAGNETIC  CHARTS  FOR  1985  -  SPHERICAL  HARMONIC 
MODELS  OF  THE  GEOMAGNETIC  FIELD  AND  ITS  SECULAR 
VARIATION,  GEOPHYS.  J.  R.  ASTR.  SOC.  (1986)  87, 

PP  1143-1157 

DEFENSE  MAPPING  AGENCY  TECHNICAL  REPORT,  TR  8350.2: 

DEPARTMENT  OF  DEFENSE  WORLD  GEODETIC  SYSTEM  1984, 
SEPT.  30  (1987) 

JOSEPH  C.  CAIN,  ET  AL.;  A  PROPOSED  MODEL  FOR  THE 

INTERNATIONAL  GEOMAGNETIC  REFERENCE  FIELD  -  1965, 

J.  GEOMAG.  AND  GEOELECT.  VOL.  19,  NO.  4,  PP  335-355 
(1967)  (SEE  APPENDIX) 

ALFRED  J.  ZMUDA,  WORLD  MAGNETIC  SURVEY  1957-1969, 
INTERNATIONAL  ASSOCIATION  OF  GEOMAGNETISM  AND 
AERONOMY  (lAGA)  BULLETIN  #28,  PP  186-188  (1971) 

JOHN  M.  QUINN,  RACHEL  J.  COLEMAN,  MICHAEL  R.  PECK,  AND 
STEPHEN  E.  LAUBER;  THE  JOINT  US/UK  1990  EPOCH 
WORLD  MAGNETIC  MODEL,  TECHNICAL  REPORT  NO.  304, 

NAVAL  OCEANOGRAPHIC  OFFICE  (1991) 

JOHN  M.  QUINN,  RACHEL  J.  COLEMAN,  DONALD  L.  SHIEL,  AND 
JOHN  M.  NIGRO;  THE  JOINT  US/UK  1995  EPOCH  WORLD 
MAGNETIC  MODEL,  TECHNICAL  REPORT  NO.  314,  NAVAL 
OCEANOGRAPHIC  OFFICE  (1995) 


C*********************************************************************** 


PARAMETER  DESCRIPTIONS: 


SEMIMAJOR  AXIS  OF  WGS-84  ELLIPSOID  (KM) 

SEMIMINOR  AXIS  OF  WGS-84  ELLIPSOID  (KM) 

MEAN  RADIUS  OF  IAU-66  ELLIPSOID  (KM) 

SCHMIDT  NORMALIZATION  FACTORS 

GAUSS  COEFFICIENTS  OF  MAIN  GEOMAGNETIC  MODEL  (NT) 

GAUSS  COEFFICIENTS  OF  SECULAR  GEOMAGNETIC  MODEL  (NT/YR) 
TIME  ADJUSTED  GEOMAGNETIC  GAUSS  COEFFICIENTS  (NT) 

TIME  ON  PREVIOUS  CALL  TO  GEOMAG  (YRS) 

GEODETIC  ALTITUDE  ON  PREVIOUS  CALL  TO  GEOMAG  (YRS) 
GEODETIC  LATITUDE  ON  PREVIOUS  CALL  TO  GEOMAG  (DEG.) 


C  OLON  - 

C  TIME 

C 

C  ALT 

C  GLAT  - 

C  GLON 

C  EPOCH  - 

C  DTR 

C  SP(M)  - 

C  CP(M)  - 

C  ST 

C  CT 

C  R 

C  CA 

C  SA 

C  BR 

C  BT 

C  BP 

C  P(N,M)  - 

C  PP(N)  - 

C  DP(N,M)- 

C  BX 

C  BY 

C  BZ 

C  BH 

C  DEC 

C 

c 

C  DIP 

C 
C 

C  TI 

C  GV 

C 
C 
C 
C 

C  MAXDEG  - 

C  MOXORD  - 

C 
C 

Q*********************************************************************** 

C 

c 

C  NOTE:  THIS  VERSION  OF  GEOMAG  USES  THE  WMM-95  GEOMAGNETIC 
C  MODEL  REFERENCED  TO  THE  WGS-84  GRAVITY  MODEL  ELLIPSOID 

C 
C 

Q*********************************************************************** 

c 

C 

c 

c 

c 

c 

c 

c 

c 

Q*********************************************************************** 

c 

c 

C  INITIALIZATION  MODULE 

C 

C 

Q*********************************************************************** 

C 

C 

SUBROUTINE  GECIMAG(MAXDEG) 

C 

C 

DIMENSION  C(0:12,0:12),CD(0:12,0:12),TC(0:12,0:12) 


GEODETIC  LONGITUDE  ON  PREVIOUS  CALL  TO  GEOHAG  (DEG.) 


COMPUTATION  TIME  (YRS)  (INPUT) 
(EG,  1  JULY  1995  =  1995.500) 

GEODETIC  ALTITUDE  (KM)  (INPUT) 
GEODETIC  LATITUDE  (DEG.)  (INPUT) 
GEODETIC  LONGITUDE  (DEG.)  (INPUT) 


BASE  TIME  OF  GEOMAGNETIC  MODEL  (YRS) 

DEGREE  TO  RADIAN  CONVERSION 

SINE  OF  (M*SPHERICAL  COORD.  LONGITUDE) 

COSINE  OF  (M*SPHERICAL  COORD.  LONGITUDE) 

SINE  OF  (SPHERICAL  COORD.  LATITUDE) 

COSINE  OF  (SPHERICAL  COORD.  LATITUDE) 

SPHERICAL  COORDINATE  RADIAL  POSITION  (KM) 

COSINE  OF  SPHERICAL  TO  GEODETIC  VECTOR  ROTATION  ANGLE 
SINE  OF  SPHERICAL  TO  GEODETIC  VECTOR  ROTATION  ANGLE 
RADIAL  COMPONENT  OF  GEOMAGNETIC  FIELD  (NT) 

THETA  COMPONENT  OF  GEOMAGNETIC  FIELD  (NT) 

PHI  COMPONENT  OF  GEOMAGNETIC  FIELD  (NT) 

ASSOCIATED  LEGENDRE  POLYNOMIALS  (UNNORMALIZED) 
ASSOCIATED  LEGENDRE  POLYNOMIALS  FOR  M=1  (UNNORMALIZED) 
THETA  DERIVATIVE  OF  P(N,H)  (UN NORMAL I ZED) 

NORTH  GEOMAGNETIC  COMPONENT  (NT) 

EAST  GEOMAGNETIC  COMPONENT  (NT) 

VERTICALLY  DOWN  GEOMAGNETIC  COMPONENT  (NT) 

HORIZONTAL  GEOMAGNETIC  COMPONENT  (NT) 

GEOMAGNETIC  DECLINATION  (DEG.) 

EAST-POSITIVE  ANGLES 

UEST=NEGATIVE  ANGLES 
GEOMAGNETIC  INCLINATION  (DEG.) 

DOWN=POSITIVE  ANGLES 
UP=NEGATIVE  ANGLES 
GEOMAGNETIC  TOTAL  INTENSITY  (NT) 

GEOMAGNETIC  GRID  VARIATION  (DEG.) 

REFERENCED  TO  GRID  NORTH 
GRID  NORTH  REFERENCED  TO  0  MERIDIAN 
OF  A  POLAR  STEREOGRAPHIC  PROJECTION 
(ARCTIC/ ANTARCTIC  ONLY) 

MAXIMUM  DEGREE  OF  SPHERICAL  HARMONIC  MODEL  (INPUT) 
MAXIMUM  ORDER  OF  SPHERICAL  HARMONIC  MODEL 


(OUTPUT) 

(OUTPUT) 

(OUTPUT) 

(OUTPUT) 
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DIMENSION  P(0:12,0:12),DP(0:12,0:12),SNORM(0:12,0:12) 
DIMENSION  SP<0:12),CP(0:12),FN(0:12),FM(0:12),PP(0:12) 
REAL  K(0:12,0:12) 

EQUIVALENCE  (SNORM,P) 

C 

C 


DATA  EPOCH/1995.0/ 


DATA  C/ 

0.0,- 

29682.1,  -2194.7, 

1318.8, 

940.0, 

-209.5, 

* 

68.5, 

78.0,  24.7, 

2.9, 

-2.9, 

1.7, 

* 

■1.8, 

5315.6,  -1782.2, 

3078.6, 

-2273.6, 

782.9, 

* 

354.0, 

65.6,  -68.1, 

3.4, 

7.5, 

-3.3, 

* 

-1.6, 

0.9,  -2359.1, 

-418.6, 

1685.7, 

1246.9, 

* 

290.9, 

238.2,  64.1, 

0.1, 

-1.5, 

0.4, 

* 

2.8, 

-3.6,  -0.1, 

-261.1, 

301.0, 

-416.5, 

* 

766.3, 

-418.9,  -122.1, 

-169.1, 

29.6, 

-9.6, 

* 

-10.3, 

-4.3,  1.2, 

-0.5, 

259.4, 

-230.9, 

* 

99.8, 

-306.1,  113.8, 

-162.8, 

-0.5, 

6.0, 

* 

-16.5, 

9.7,  -3.1, 

-0.6, 

0.8, 

43.7, 

* 

157.6, 

-150.1,  -59.2, 

104.4, 

-23.3, 

16.5, 

* 

8.7, 

2.6,  -2.3, 

2.4, 

0.1, 

0.2, 

* 

-15.2, 

74.3,  69.4, 

-55.3, 

3.0, 

33.3, 

* 

-91.0, 

9.2,  3.6, 

-2.4, 

2.8, 

-0.7, 

* 

0.5, 

-76.1,  -24.5, 

1.6, 

20.0, 

16.5, 

* 

-23.6, 

-6.8,  -2.4, 

-4.9, 

6.8, 

0.7, 

* 

-0.8, 

0.4,  14.9, 

-19.5, 

6.3, 

-20.4, 

4r 

12.2, 

7.0,  -19.0, 

-8.8, 

-8.5, 

-0.5, 

* 

4.1, 

1.3,  -0.4, 

-19.8, 

14.6, 

10.9, 

* 

-7.5, 

-6.8,  9.3, 

7.7, 

-8.1, 

2.6, 

* 

-6.5, 

3.6,  -0.3, 

0.3, 

3.2, 

1.7, 

it 

2.9, 

5.6,  -3.4, 

-0.7, 

-2.9, 

2.3, 

it 

-1.6, 

-6.6,  0.6, 

2.2, 

0.2, 

0.3, 

it 

1.0, 

-3.6,  -1.4, 

1.9, 

0.2, 

-1.3, 

it 

-2.4, 

-0.6,  -2.2, 

1.3, 

4.2, 

0.4, 

it 

0.3, 

1.4,  0.8, 

-3.0, 

0.7, 

0.5, 

* 

-0.8, 

0.6,  0.1, 

-1.3, 

1 

o 

0.9, 

it 

0.6/ 

C 

C 


DATA  CD/ 

0.0, 

17.6,- 

13.7, 

0.8, 

1.2, 

0.9,  0.4,  -0.3, 

it 

0.1, 

0.0, 

0.0, 

0.0, 

0.0, 

-18.0,  13.2,  4.0, 

it 

-6.6, 

1.1. 

0.5, 

-0.3, 

-1.1, 

0.0,  0.0,  0.0, 

it 

0.0, 

0.0,- 

14.6, 

-7.2, 

-0.3, 

-0.5,  -6.8,  -1.4, 

it 

0.3, 

-0.5, 

0.4, 

0.0, 

0.0, 

0.0,  0.0,  4.0, 

it 

2.2, 

-12.6, 

-8.5, 

0.3, 

-1.7, 

2.1,  0.5,  0.3, 

it 

0.0, 

0.0, 

0.0, 

0.0, 

1.3, 

1.0,  2.5,  -1.2, 

* 

-4.5, 

0.0, 

0.0, 

1.3, 

-1.3, 

o 

o 

o 

o 

o 

o 

★ 

0.0, 

0.5, 

1.5, 

0.6, 

1.7, 

0.6,  2.1,  -0.4, 

* 

0.1, 

0.5, 

0.0, 

0.0, 

0.0, 

0.0,  0.7,  -1.5, 

it 

-0.5, 

-0.7, 

1.1, 

2.6, 

-0.4, 

0.0,  0.4,  0.0, 

* 

0.0, 

0.0, 

0.0, 

0.3, 

0.0, 

0.7,  -0.6,  0.1, 

* 

-0.6, 

-0.4, 

-0.9, 

-0.9, 

0.0, 

0.0,  0.0,  0.0, 

* 

0.4, 

-0.3, 

0.1, 

0.8, 

-0.1, 

-1.3,  -0.9,  -1.1, 

* 

0.1, 

56*0.0/ 

C 

C 

C 

C 

C 

C  INITIALIZE  CONSTANTS 

C 

C 

IF  (MAXOEG  .GT,  12)  MAXDEG=12 

MAXORD=MAXDEG 

PI=3. 14159265359 

DTR=PI/180.0 

SP(0)=0. 

CP(0)=1. 

P<0,0)=1. 

PP(0)=1. 
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DP(0,0)=0. 

A=6378.137 

8=6356.7523142 

RE=6371.2 

A2=A**2 

B2=B**2 

C2=A2-B2 

A4=A2**2 

B4=B2**2 

C4=A4-B4 

CONVERT  SCHMIDT  NORMALIZED  GAUSS  COEFFICIENTS  TO  UNNORMALIZED 

SNORM<0,0)=1. 

DO  20  N=1,MAXORD 

SNORM<N,0)=SNORM(M-1,0)*FLOAT(2*N-1)/FLOAT(N) 

J=2 

DO  10  M=0,N 

K(N,H)=FLOAT<(N-1)**2-M**2)/FLOAT((2*N-1)*(2*N-3)) 

IF  (M  .GT.  0)  THEN 
FLNMJ=FLOAT((N-M+1 )*J)/FLOAT(N+M) 
SNORM(N,M)=SNORM(N,M-1)*SQRT(FLNMJ) 

J=1 

C<M-1,N)=SNORM(M,M)*C(M-1,N) 

CD(M-1,N)=SNORM(N,M)*CD(M-1,N) 

END  IF 

C(N,M)=SNORM(N,M)*C(N.M) 

CD(M,M)=SNORM(N,M)*CD(N,M) 

10  CONTINUE 

FN{N)=FLOAT(N+1) 

FM(N)=FLOAT(N) 

20  CONTINUE 
IC(1,1)=0. 


OTIME=-1000. 

OALT=-1000. 

OLAT=-1000. 

OL0N=-1000. 


r*********************************************************************** 


PROCESSING  MODULE 


0*********************************************************************** 

C 

c 

ENTRY  GEOMG1(ALT,GLAT,GLON, TIME, DEC, DIP, TI,GV) 


DT=TIME-EPOCH 
IF  (OTIME  .LT.  0.  .AND. 
PRINT  *,  '  ■ 


.LT.  0.  .OR.  DT 


5.))  THEN 


PRINT  * 
PRINT  * 
PRINT  * 
PRINT  * 
PRINT  * 
PRINT  * 
PRINT  * 
PRINT  * 
PRINT  * 
PRINT  * 
PRINT  * 
END  IF 


WARNING  -  TIME  EXTENDS  BEYOND  MODEL  5-YEAR  LIFE  SPAN' 
CONTACT  DMA  FOR  PRODUCT  UPDATES:' 

I 

DEFENSE  MAPPING  AGENCY' 

ATTN:  Code  PRS' 

8613  LEE  HIGHWAY' 

FAIRFAX,  VA  22031-2137' 

(703)285-9197,  AUTOVON  356-9197' 

I 

EPOCH  =  ', EPOCH 
TIME  =  ',TIME 
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r>  o  o 


c 

c 

RLOM=GLON*DTR 

RLAT=GLAT*DTR 

SRLON=SIN(RLON} 

SRLAT=SIN(RLAT) 

CRLON=COS(RLON) 

CRLAT=COS(RLAT) 

SRLON2=SRLON**2 
SRLAT2=SRLAT**2 
CRLON2=CRLON**2 
CRLAT2=CRLAT**2 
SP(1)=SRL0N 
CP(1)=CRL0N 

CONVERT  FROM  GEODETIC  COORDS.  TO  SPHERICAL  COORDS. 

IF  (ALT  .NE.  OALT  .OR.  GLAT  .ME.  OLAT)  THEN 
Q=SQRT(A2-C2*SRLAT2) 

Q1=ALT*Q 

Q2=((Q1+A2)/(Q1+B2))**2 
CT=SRLAT/SQRT(Q2*CRLAT2+SRLAT2) 

R2=ALT**2+2.0*Q1+(A4-C4*SRLAT2)/Q**2 
R=SQRT(R2) 

D=SQRT(A2*CRLAT2+B2*SRLAT2) 

CA=(ALT+D)/R 
SA=C2*CRLAT*SRLAT/(R*0) 

END  IF 
C 
C 

IF  (GLON  .NE.  OLON)  THEN 
DO  40  M=2,HAXORD 

SP(M)=SP(1)*CP(H-1)+CP(1)*SP(M-1) 

CP(H)=CP(1)*CP(M-1>-SP(1)*SP(M-1) 

40  CONTINUE 
END  IF 
C 
C 

AOR=RE/R 

AR=AOR**2 

C 

C 

BR=0. 

BT=0. 

BP=0. 

BPP=0. 

C 

C 

DO  70  N=1,MAXORD 
AR=AR*AOR 
DO  60  H=0,N 
C 

C  COMPUTE  UMNORMALIZED  ASSOCIATED  LEGENDRE  POLYNOMIALS 

C  AND  DERIVATIVES  VIA  RECURSION  RELATIONS 

C 

IF  (ALT  .NE.  OALT  .OR.  GLAT  .NE.  OLAT)  THEM 
IF  (N  .EQ.  M)  THEN 
P(N,M)=ST*P(N-1,M-1) 
DP(N,M)=ST*DP(N-1,M-1)+CT*P(M-1,M-1) 

GO  TO  50 
END  IF 

IF  (N  .EQ.  1  .AND.  M  .EQ.  0)  THEN 
P(N,M)=CT*P(M-1,M) 

DP(N,M)=CT*DP(N-1,M)-ST*P(N-1,M) 

GO  TO  50 
END  IF 

IF  (N  .GT.  1  .AND.  N  .NE.  M)  THEN 
IF  (M  .GT.  M-2)  P(N-2,M)=0.0 
IF  (M  .GT.  M-2)  DP(N-2,M)=0.0 
P(N,M)=CT*P(M-1,M)-K(N,M)*P(N-2,M) 


371 


DP(N,M)=CT*DP(N-1,M)-ST*P(N-1,M)-K(N,M)*DP(N“2,M) 

END  IF 
END  IF 

50  CONTINUE 
C 

C  TIME  ADJUST  THE  GAUSS  COEFFICIENTS 

C 

IF  (TIME  .NE.  OTIME)  THEN 
TC<N,M)=C(N,M)+DT*CD(N,M) 

IF  (M  .NE.  0)  THEN 
TC(M-1,N)=C(M-1,N)+DT*CD(M-1,N) 

END  IF 
ENDIF 
C 

C  ACCUMULATE  TERMS  OF  THE  SPHERICAL  HARMONIC  EXPANSIONS 

C 

PAR=AR*P(N,M) 

IF  (M  .EQ.  0)  THEN 
TEMP1=TC(N,M)*CP<M) 

TEMP2=TC(N,M)*SP(M) 

ELSE 

TEMP1=TC(N,M)*CP(M)+TC(M-1,N)*SP(M) 

TEMP2*TC(N,M)*SP(M)-TC(M*1,N)*CP(M) 

ENDIF 

BT=BT-AR*TEMP1*DP(N,M) 

BP=BP+FM(M)*TEMP2*PAR 

BR=BR+FN(N)*TEMP1*PAR 

C 

C  SPECIAL  CASE:  NORTH/SOUTH  GEOGRAPHIC  POLES 

C 

IF  (ST  .EQ.  0.0  .AND.  M  .EQ.  1)  THEN 

IF  (N  .EQ.  1)  THEN 

PP(N)=PP(N-1) 

ELSE 

PP(N)=CT*PP(N-1)-IC(N,M)*PP(N-2) 

ENDIF 

PARP=AR*PP(N) 

BPP=BPP+FM(M)*TEMP2*PARP 

ENDIF 

C 

C 

60  CONTINUE 
70  CONTINUE 
C 
C 

IF  (ST  .EQ.  0.0)  THEN 

BP=BPP 

ELSE 

BP=8P/ST 

ENDIF 

C 

C  ROTATE  MAGNETIC  VECTOR  COMPONENTS  FROM  SPHERICAL  TO 

C  GEODETIC  COORDINATES 

C 

BX=-BT*CA-BR*SA 

BY=BP 

BZ=BT*SA-BR*CA 

C 

C  COMPUTE  DECLINATION  (DEC),  INCLINATION  (DIP)  AND 

C  TOTAL  INTENSITY  (TI) 

C 

BH=SQRT(BX**2+BY**2) 

TI=SQRT(BH**2+BZ**2) 

DEC=ATAN2(BY,BX)/DTR 

DIP=ATAN2(BZ,BH)/DTR 

C 

C  COMPUTE  MAGNETIC  GRID  VARIATION  IF  THE  CURRENT 

C  GEODETIC  POSITION  IS  IN  THE  ARCTIC  OR  ANTARCTIC 

C  (I.E.  GLAT  >  +55  DEGREES  OR  GLAT  <  -55  DEGREES) 

C 

C  OTHERWISE,  SET  MAGNETIC  GRID  VARIATION  TO  -999.0 
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c 

GV=-999.0 

IF  (ABS(GLAT)  -GE.  55.)  THEM 

IF  (GLAT  .GT.  0.  .AMD.  GLOM  .GE.  0.)  GV=DEC-GLON 

IF  (GLAT  .GT.  0.  .AND.  GLOM  .LT.  0.)  GV=DEC+ABS(GLOM) 

IF  (GLAT  .LT.  0.  .AND.  6LON  .GE.  0.)  GV=DEC+GLON 

IF  (GLAT  .LT.  0.  .AND.  6LON  .LT.  0.)  GV=DEC-ABS(GLON) 

IF  (GV  .GT.  +180.)  GV=GV-360. 

IF  (GV  .LT.  -180.)  GV=GV+360. 

END  IF 
C 
C 

OTIME=TIME 

OALT=ALT 

OLAT=GLAT 

OLON=GLON 

C 

C 

RETURN 

C 

C 

END 
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